The functional transfer of genes from the mitochondria to the nucleus: the effects of selection, mutation, population size and rate of self-fertilization

Viruses and Mitochondrial Dysfunction in Neurodegeneration and Cognition: An Evolutionary Perspective

Mitochondria, the cellular powerhouses with bacterial evolutionary origins, play a pivotal role in maintaining neuronal function and cognitive health. Several viruses have developed sophisticated mechanisms to target and disrupt mitochondrial function which contribute to cognitive decline and neurodegeneration. The interplay between viruses and mitochondria might be traced to their co-evolutionary history with bacteria and may reflect ancient interactions that have shaped modern mitochondrial biology.

The Spiral Model of Evolution: Stable Life Forms of Organisms and Unstable Life Forms of Cancers

If one must prioritize among the vast array of contributing factors to cancer evolution, environmental-stress-mediated chromosome instability (CIN) should easily surpass individual gene mutations. CIN leads to the emergence of genomically unstable life forms, enabling them to grow dominantly within the stable life form of the host. In contrast, stochastic gene mutations play a role in aiding the growth of the cancer population, with their importance depending on the initial emergence of the new system. Furthermore, many specific gene mutations among the many available can perform this function, decreasing the clinical value of any specific gene mutation. Since these unstable life forms can respond to treatment differently than stable ones, cancer often escapes from drug treatment by forming new systems, which leads to problems during the treatment for patients. To understand how diverse factors impact CIN-mediated macroevolution and genome integrity-ensured microevolution, the concept of two-phased cancer evolution is used to reconcile some major characteristics of cancer, such as bioenergetic, unicellular, and multicellular evolution. Specifically, the spiral of life function model is proposed, which integrates major historical evolutionary innovations and conservation with information management. Unlike normal organismal evolution in the microevolutionary phase, where a given species occupies a specific location within the spiral, cancer populations are highly heterogenous at multiple levels, including epigenetic levels. Individual cells occupy different levels and positions within the spiral, leading to supersystems of mixed cellular populations that exhibit both macro and microevolution. This analysis, utilizing karyotype to define the genetic networks of the cellular system and CIN to determine the instability of the system, as well as considering gene mutation and epigenetics as modifiers of the system for information amplification and usage, explores the high evolutionary potential of cancer. It provides a new, unified understanding of cancer as a supersystem, encouraging efforts to leverage the dynamics of CIN to develop improved treatment options. Moreover, it offers a historically contingent model for organismal evolution that reconciles the roles of both evolutionary innovation and conservation through macroevolution and microevolution, respectively.

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.

Ecological Predictors of Organelle Genome Evolution: Phylogenetic Correlations with Taxonomically Broad, Sparse, Unsystematized Data

Comparative analysis of variables across phylogenetically linked observations can reveal mechanisms and insights in evolutionary biology. As the taxonomic breadth of the sample of interest increases, challenges of data sparsity, poor phylogenetic resolution, and complicated evolutionary dynamics emerge. Here, we investigate a cross-eukaryotic question where all these problems exist: which organismal ecology features are correlated with gene retention in mitochondrial and chloroplast DNA (organelle DNA or oDNA). Through a wide palette of synthetic control studies, we first characterize the specificity and sensitivity of a collection of parametric and non-parametric phylogenetic comparative approaches to identify relationships in the face of such sparse and awkward datasets. This analysis is not directly focused on oDNA, and so provides generalizable insights into comparative approaches with challenging data. We then combine and curate ecological data coupled to oDNA genome information across eukaryotes, including a new semi-automated approach for gathering data on organismal traits from less systematized open-access resources including encyclopedia articles on species and taxa. The curation process also involved resolving several issues with existing datasets, including enforcing the clade-specificity of several ecological features and fixing incorrect annotations. Combining this unique dataset with our benchmarked comparative approaches, we confirm support for several known links between organismal ecology and organelle gene retention, identify several previously unidentified relationships constituting possible ecological contributors to oDNA genome evolution, and provide support for a recently hypothesized link between environmental demand and oDNA retention. We, with caution, discuss the implications of these findings for organelle evolution and of this pipeline for broad comparative analyses in other fields.

Paternal Inheritance of Mitochondrial DNA May Lead to Dioecy in Conifers

In angiosperms cytoplasmic DNA is typically passed on maternally through ovules. Genes in the mtDNA may cause male sterility. When male-sterile (female) cytotypes produce more seeds than cosexuals, they pass on more copies of their mtDNA and will co-occur with cosexuals with a neutral cytotype. Cytoplasmic gynodioecy is a well-known phenomenon in angiosperms, both in wild and crop plants. In some conifer families (e.g. Pinaceae) mitochondria are also maternally inherited. However in some other families (e.g. Taxaceae and Cupressaceae) mtDNA is paternally inherited through the pollen. With paternal mtDNA inheritance, male cytotypes that produce more pollen than cosexuals are expected to co-occur with cosexuals. This is uncharted territory. An ESS model shows that the presence of male cytotypes selects for more female allocation in the cosexual, i.e. for sexual specialisation. An allele that switches sex from male to female can then invade. This leads to rapid loss of the neutral cytotype of the cosexual, fixation of the male cytotype and dioecy with 50% males and 50% females. The models suggest that paternal inheritance of mtDNA facilitates the evolution dioecy. Consistent with this hypothesis the Pinaceae are 100% monoecious, while dioecy is common in the Taxaceae family and in the genus Juniperus (Cupressaceae). However, no reliable data are yet available on both mode of inheritance of mtDNA and gender variation of the same species. When cosexuals benefit from reproductive assurance (high selfing rate, low inbreeding depression, low fertilisation) they maintain themselves next to males and females. This predicted pattern with three sex types present in the same population is observed in conifers in nature.

Connecting Species-Specific Extents of Genome Reduction in Mitochondria and Plastids

Mitochondria and plastids have both dramatically reduced their genomes since the endosymbiotic events that created them. The similarities and differences in the evolution of the two organelle genome types have been the target of discussion and investigation for decades. Ongoing work has suggested that similar mechanisms may modulate the reductive evolution of the two organelles in a given species, but quantitative data and statistical analyses exploring this picture remain limited outside of some specific cases like parasitism. Here, we use cross-eukaryote organelle genome data to explore evidence for coevolution of mitochondrial and plastid genome reduction. Controlling for differences between clades and pseudoreplication due to relatedness, we find that extents of mtDNA and ptDNA gene retention are related to each other across taxa, in a generally positive correlation that appears to differ quantitatively across eukaryotes, for example, between algal and nonalgal species. We find limited evidence for coevolution of specific mtDNA and ptDNA gene pairs, suggesting that the similarities between the two organelle types may be due mainly to independent responses to consistent evolutionary drivers.

Collective mitochondrial dynamics resolve conflicting cellular tensions: From plants to general principles

Mitochondria play diverse and essential roles in eukaryotic cells, and plants are no exception. Plant mitochondria have several differences from their metazoan and fungal cousins: they often exist in a fragmented state, move rapidly on actin rather than microtubules, have many plant-specific metabolic features and roles, and usually contain only a subset of the complete mtDNA genome, which itself undergoes frequent recombination. This arrangement means that exchange and complementation is essential for plant mitochondria, and recent work has begun to reveal how their collective dynamics and resultant "social networks" of encounters support this exchange, connecting plant mitochondria in time rather than in space. This review will argue that this social network perspective can be extended to a "societal network", where mitochondrial dynamics are an essential part of the interacting cellular society of organelles and biomolecules. Evidence is emerging that mitochondrial dynamics allow optimal resolutions to competing cellular priorities; we will survey this evidence and review potential future research directions, highlighting that plant mitochondria can help reveal and test principles that apply across other kingdoms of life. In parallel with this fundamental cell biology, we also highlight the translational "One Health" importance of plant mitochondrial behaviour - which is exploited in the production of a vast amount of crops consumed worldwide - and the potential for multi-objective optimisation to understand and rationally re-engineer the evolved resolutions to these tensions.

Mitochondrial genomes revisited: why do different lineages retain different genes?

The mitochondria contain their own genome derived from an alphaproteobacterial endosymbiont. From thousands of protein-coding genes originally encoded by their ancestor, only between 1 and about 70 are encoded on extant mitochondrial genomes (mitogenomes). Thanks to a dramatically increasing number of sequenced and annotated mitogenomes a coherent picture of why some genes were lost, or relocated to the nucleus, is emerging. In this review, we describe the characteristics of mitochondria-to-nucleus gene transfer and the resulting varied content of mitogenomes across eukaryotes. We introduce a 'burst-upon-drift' model to best explain nuclear-mitochondrial population genetics with flares of transfer due to genetic drift.

The Mighty NUMT: Mitochondrial DNA Flexing Its Code in the Nuclear Genome

Nuclear-mitochondrial DNA segments (NUMTs) are mitochondrial DNA (mtDNA) fragments that have been inserted into the nuclear genome. Some NUMTs are common within the human population but most NUMTs are rare and specific to individuals. NUMTs range in size from 24 base pairs to encompassing nearly the entire mtDNA and are found throughout the nuclear genome. Emerging evidence suggests that the formation of NUMTs is an ongoing process in humans. NUMTs contaminate sequencing results of the mtDNA by introducing false positive variants, particularly heteroplasmic variants present at a low variant allele frequency (VAF). In our review, we discuss the prevalence of NUMTs in the human population, the potential mechanisms of de novo NUMT insertion via DNA repair mechanisms, and provide an overview of the existing approaches for minimizing NUMT contamination. Apart from filtering known NUMTs, both wet lab-based and computational methods can be used to minimize the contamination of NUMTs in analyses of human mtDNA. Current approaches include: (1) isolating mitochondria to enrich for mtDNA; (2) applying basic local alignment to identify NUMTs for subsequent filtering; (3) bioinformatic pipelines for NUMT detection; (4) k-mer-based NUMT detection; and (5) filtering candidate false positive variants by mtDNA copy number, VAF, or sequence quality score. Multiple approaches must be applied in order to effectively identify NUMTs in samples. Although next-generation sequencing is revolutionizing our understanding of heteroplasmic mtDNA, it also raises new challenges with the high prevalence and individual-specific NUMTs that need to be handled with care in studies of mitochondrial genetics.

Cellular and environmental dynamics influence species-specific extents of organelle gene retention

Mitochondria and plastids rely on many nuclear-encoded genes, but retain small subsets of the genes they need to function in their own organelle DNA (oDNA). Different species retain different numbers of oDNA genes, and the reasons for these differences are not completely understood. Here, we use a mathematical model to explore the hypothesis that the energetic demands imposed by an organism's changing environment influence how many oDNA genes it retains. The model couples the physical biology of cell processes of gene expression and transport to a supply-and-demand model for the environmental dynamics to which an organism is exposed. The trade-off between fulfilling metabolic and bioenergetic environmental demands, and retaining genetic integrity, is quantified for a generic gene encoded either in oDNA or in nuclear DNA. Species in environments with high-amplitude, intermediate-frequency oscillations are predicted to retain the most organelle genes, whereas those in less dynamic or noisy environments the fewest. We discuss support for, and insight from, these predictions with oDNA data across eukaryotic taxa, including high oDNA gene counts in sessile organisms exposed to day-night and intertidal oscillations (including plants and algae) and low counts in parasites and fungi.

Global mitochondrial protein import proteomics reveal distinct regulation by translation and translocation machinery

Most mitochondrial proteins are translated in the cytosol and imported into mitochondria. Mutations in the mitochondrial protein import machinery cause human pathologies. However, a lack of suitable tools to measure protein uptake across the mitochondrial proteome has prevented the identification of specific proteins affected by import perturbation. Here, we introduce mePROD^mt, a pulsed-SILAC based proteomics approach that includes a booster signal to increase the sensitivity for mitochondrial proteins selectively, enabling global dynamic analysis of endogenous mitochondrial protein uptake in cells. We applied mePROD^mt to determine protein uptake kinetics and examined how inhibitors of mitochondrial import machineries affect protein uptake. Monitoring changes in translation and uptake upon mitochondrial membrane depolarization revealed that protein uptake was extensively modulated by the import and translation machineries via activation of the integrated stress response. Strikingly, uptake changes were not uniform, with subsets of proteins being unaffected or decreased due to changes in translation or import capacity.

•

Proteomics approach to quantify protein uptake into mitochondria

•

Determination of protein uptake rates for >700 mitochondrial proteins

•

Characterization of differential protein uptake changes during mitochondrial stress

•

Protein translation and translocation integrate to change protein uptake

Rapid shifts in mitochondrial tRNA import in a plant lineage with extensive mitochondrial tRNA gene loss

In most eukaryotes, transfer RNAs (tRNAs) are one of the very few classes of genes remaining in the mitochondrial genome, but some mitochondria have lost these vestiges of their prokaryotic ancestry. Sequencing of mitogenomes from the flowering plant genus Silene previously revealed a large range in tRNA gene content, suggesting rapid and ongoing gene loss/replacement. Here, we use this system to test longstanding hypotheses about how mitochondrial tRNA genes are replaced by importing nuclear-encoded tRNAs. We traced the evolutionary history of these gene loss events by sequencing mitochondrial genomes from key outgroups (Agrostemma githago and Silene [=Lychnis] chalcedonica). We then performed the first global sequencing of purified plant mitochondrial tRNA populations to characterize the expression of mitochondrial-encoded tRNAs and the identity of imported nuclear-encoded tRNAs. We also confirmed the utility of high-throughput sequencing methods for the detection of tRNA import by sequencing mitochondrial tRNA populations in a species (Solanum tuberosum) with known tRNA trafficking patterns. Mitochondrial tRNA sequencing in Silene revealed substantial shifts in the abundance of some nuclear-encoded tRNAs in conjunction with their recent history of mt-tRNA gene loss and surprising cases where tRNAs with anticodons still encoded in the mitochondrial genome also appeared to be imported. These data suggest that nuclear-encoded counterparts are likely replacing mitochondrial tRNAs even in systems with recent mitochondrial tRNA gene loss, and the redundant import of a nuclear-encoded tRNA may provide a mechanism for functional replacement between translation systems separated by billions of years of evolutionary divergence.

Gene tree and species tree reconciliation with endosymbiotic gene transfer

MOTIVATION

It is largely established that all extant mitochondria originated from a unique endosymbiotic event integrating an α-proteobacterial genome into an eukaryotic cell. Subsequently, eukaryote evolution has been marked by episodes of gene transfer, mainly from the mitochondria to the nucleus, resulting in a significant reduction of the mitochondrial genome, eventually completely disappearing in some lineages. However, in other lineages such as in land plants, a high variability in gene repertoire distribution, including genes encoded in both the nuclear and mitochondrial genome, is an indication of an ongoing process of Endosymbiotic Gene Transfer (EGT). Understanding how both nuclear and mitochondrial genomes have been shaped by gene loss, duplication and transfer is expected to shed light on a number of open questions regarding the evolution of eukaryotes, including rooting of the eukaryotic tree.

RESULTS

We address the problem of inferring the evolution of a gene family through duplication, loss and EGT events, the latter considered as a special case of horizontal gene transfer occurring between the mitochondrial and nuclear genomes of the same species (in one direction or the other). We consider both EGT events resulting in maintaining (EGTcopy) or removing (EGTcut) the gene copy in the source genome. We present a linear-time algorithm for computing the DLE (Duplication, Loss and EGT) distance, as well as an optimal reconciled tree, for the unitary cost, and a dynamic programming algorithm allowing to output all optimal reconciliations for an arbitrary cost of operations. We illustrate the application of our EndoRex software and analyze different costs settings parameters on a plant dataset and discuss the resulting reconciled trees.

AVAILABILITY AND IMPLEMENTATION

EndoRex implementation and supporting data are available on the GitHub repository via https://github.com/AEVO-lab/EndoRex.

How are the mitochondrial genomes reorganized in Hexapoda? Differential evolution and the first report of convergences within Hexapoda

The evolution of the Hexapoda mitochondrial genome has been the focus of several genetic and evolutionary studies over the last decades. However, they have concentrated on certain taxonomic orders of economic or health importance. The recent increase of mitochondrial genomes sequencing of diverse taxonomic orders generates an important opportunity to clarify the evolution of this group of organisms. However, there is no comparative study that investigates the evolution of the Hexapoda mitochondrial genome. In order to verify the level of rearrangement and the mitochondrial genome evolution, we performed a comparative genomic analysis of the Hexapoda mitochondrial genome available in the NCBI database. Using a combination of bioinformatics methods to carefully examine the mitochondrial gene rearrangements in 1198 Hexapoda species belonging to 32 taxonomic orders, we determined that there is a great variation in the rate of rearrangement by gene and by taxonomic order. A higher rate of genetic reassortment is observed in Phthiraptera, Thysanoptera, Protura, and Hymenoptera; compared to other taxonomic orders. Twenty-four events of convergence in the genetic order between different taxonomic orders were determined, most of them not previously reported; which proves the great evolutionary dynamics within Hexapoda.

Predicting the Evolution of Syntenies—An Algorithmic Review

Syntenies are genomic segments of consecutive genes identified by a certain conservation in gene content and order. The notion of conservation may vary from one definition to another, the more constrained requiring identical gene contents and gene orders, while more relaxed definitions just require a certain similarity in gene content, and not necessarily in the same order. Regardless of the way they are identified, the goal is to characterize homologous genomic regions, i.e., regions deriving from a common ancestral region, reflecting a certain gene co-evolution that can enlighten important functional properties. In addition of being able to identify them, it is also necessary to infer the evolutionary history that has led from the ancestral segment to the extant ones. In this field, most algorithmic studies address the problem of inferring rearrangement scenarios explaining the disruption in gene order between segments with the same gene content, some of them extending the evolutionary model to gene insertion and deletion. However, syntenies also evolve through other events modifying their content in genes, such as duplications, losses or horizontal gene transfers, i.e., the movement of genes from one species to another. Although the reconciliation approach between a gene tree and a species tree addresses the problem of inferring such events for single-gene families, little effort has been dedicated to the generalization to segmental events and to syntenies. This paper reviews some of the main algorithmic methods for inferring ancestral syntenies and focus on those integrating both gene orders and gene trees.

Sexually antagonistic evolution of mitochondrial and nuclear linkage

Across eukaryotes, genes encoding bioenergetic machinery are located in both mitochondrial and nuclear DNA, and incompatibilities between the two genomes can be devastating. Mitochondria are often inherited maternally, and theory predicts sex-specific fitness effects of mitochondrial mutational diversity. Yet how evolution acts on linkage patterns between mitochondrial and nuclear genomes is poorly understood. Using novel mito-nuclear population-genetic models, we show that the interplay between nuclear and mitochondrial genes maintains mitochondrial haplotype diversity within populations, and selects both for sex-independent segregation of mitochondrion-interacting genes and for paternal leakage. These effects of genetic linkage evolution can eliminate male-harming fitness effects of mtDNA mutational diversity. With maternal mitochondrial inheritance, females maintain a tight mitochondrial-nuclear match, but males accumulate mismatch mutations because of the weak statistical associations between the two genomic components. Sex-independent segregation of mitochondria-interacting loci improves the mito-nuclear match. In a sexually antagonistic evolutionary process, male nuclear alleles evolve to increase the rate of recombination, whereas females evolve to suppress it. Paternal leakage of mitochondria can evolve as an alternative mechanism to improve the mito-nuclear linkage. Our modelling framework provides an evolutionary explanation for the observed paucity of mitochondrion-interacting genes on mammalian sex chromosomes and for paternal leakage in protists, plants, fungi and some animals.

Chloroplast genome assembly of Handroanthus impetiginosus: comparative analysis and molecular evolution in Bignoniaceae

Bignoniaceae species have conserved chloroplast structure, with hotspots of nucleotide diversity. Several genes are under positive selection, and can be targets for evolutionary studies. Bignoniaceae is one of the most species-rich family of woody plants in Neotropical seasonally dry forests. Here we report the assembly of Handroanthus impetiginosus chloroplast genome and evolutionary comparative analyses of ten Bignoniaceae species representing the genera for which whole-genome chloroplast sequences were available. The chloroplast genome of H. impetiginosus is 159,462 bp in size and has a similar structure compared to the other nine species. The total number of genes was slightly variable amongst the Bignoniaceae, ranging from 124 in H. impetiginosus to 144 in Anemopaegma acutifolium. The inverted repeat (IR) size was variable, ranging from 24,657 bp (Tecomaria capensis) to 40,481 bp (A. acutifolium), due to the contraction and retraction at its boundaries. However, gene boundaries were very similar among the ten species. We found 98 forward and palindromic dispersed repeats, and 85 simple sequence repeats (SSRs). In general, chloroplast sequences were highly conserved, with few nucleotide diversity hotspots in the genes accD, clpP, rpoA, ycf1, ycf2. The phylogenetic analysis based on 77 coding genes was highly consistent with Angiosperm Phylogeny Group (APG) IV. Our results also indicate that most genes are under negative selection or neutral evolution. We found no evidence of branch-site selection, implying that H. impetiginosus is not evolving faster than the other species analyzed, notwithstanding we found site positive selection signal in several genes. These genes can provide targets for evolutionary studies in Bignoniaceae and Lamiales species.

The Mitochondrial Genome Assembly of Fennel (Foeniculum vulgare) Reveals Two Different atp6 Gene Sequences in Cytoplasmic Male Sterile Accessions

Cytoplasmic male sterility (CMS) has always aroused interest among researchers and breeders, being a valuable resource widely exploited not only to breed F1 hybrid varieties but also to investigate genes that control stamen and pollen development. With the aim of identifying candidate genes for CMS in fennel, we adopted an effective strategy relying on the comparison between mitochondrial genomes (mtDNA) of both fertile and sterile genotypes. mtDNA raw reads derived from a CMS genotype were assembled in a single molecule (296,483 bp), while a draft mtDNA assembly (166,124 nucleotides, 94 contigs) was performed using male fertile sample (MF) sequences. From their annotation and alignment, two atp6-like sequences were identified. atp6-, the putative mutant copy with a 300 bp truncation at the 5'-end, was found only in the mtDNA of CMS samples, while the wild type copy (atp6+) was detected only in the MF mtDNA. Further analyses (i.e., reads mapping and Sanger sequencing), revealed an atp6+ copy also in CMS samples, probably in the nuclear DNA. However, qPCRs performed on different tissues proved that, despite its availability, atp6+ is expressed only in MF samples, while apt6- mRNA was always detected in CMS individuals. In the light of these findings, the energy deficiency model could explain the pollen deficiency observed in male sterile flower. atp6- could represent a gene whose mRNA is translated into a not-fully functional protein leading to suboptimal ATP production that guarantees essential cellular processes but not a high energy demand process such as pollen development. Our study provides novel insights into the fennel mtDNA genome and its atp6 genes, and paves the way for further studies aimed at understanding their functional roles in the determination of male sterility.

Cytonuclear Genetic Incompatibilities in Plant Speciation

Due to the endosymbiotic origin of organelles, a pattern of coevolution and coadaptation between organellar and nuclear genomes is required for proper cell function. In this review, we focus on the impact of cytonuclear interaction on the reproductive isolation of plant species. We give examples of cases where species exhibit barriers to reproduction which involve plastid-nuclear or mito-nuclear genetic incompatibilities, and describe the evolutionary processes at play. We also discuss potential mechanisms of hybrid fitness recovery such as paternal leakage. Finally, we point out the possible interplay between plant mating systems and cytonuclear coevolution, and its consequence on plant speciation.

RPGeNet v2.0: expanding the universe of retinal disease gene interactions network

RPGeNet offers researchers a user-friendly queriable tool to visualize the interactome network of visual disorder genes, thus enabling the identification of new potential causative genes and the assignment of novel candidates to specific retinal or cellular pathways. This can be highly relevant for clinical applications as retinal dystrophies affect 1:3000 people worldwide, and the causative genes are still unknown for 30% of the patients. RPGeNet is a refined interaction network interface that limits its skeleton network to the shortest paths between each and every known causative gene of inherited syndromic and non-syndromic retinal dystrophies. RPGeNet integrates interaction information from STRING, BioGRID and PPaxe, along with retina-specific expression data and associated genetic variants, over a Cytoscape.js web interface. For the new version, RPGeNet v2.0, the database engine was migrated to Neo4j graph database manager, which speeds up the initial queries and can handle whole interactome data for new ways to query the network. Further, user facilities have been introduced as the capability of saving and restoring a researcher customized network layout or as novel features to facilitate navigation and data projection on the network explorer interface. Responsiveness has been further improved by transferring some functionality to the client side.

Extensive survey of the ycf4 plastid gene throughout the IRLC legumes: Robust evidence of its locus and lineage specific accelerated rate of evolution, pseudogenization and gene loss in the tribe Fabeae

The genome organization and gene content of plastome (plastid genome) are highly conserved among most flowering plant species. Plastome variation (in size and gene order) is rare in photosynthetic species but size variation, rearrangements and gene/intron losses is attributed to groups of seed plants. Fabaceae (legume family), in particular the subfamily Papilionoideae and the inverted repeat lacking clade (IRLC), a largest legume lineage, display the most dramatic and structural change which providing an excellent model for understanding of mechanisms of genomic evolution. The IRLC comprises 52 genera and ca 4000 species divided into seven tribes. In present study, we have sampled several representatives from each tribe across the IRLC from various herbaria and field. The ycf4 gene, which plays a role in regulating and assembly of photosystem I, is more variable in the tribe Fabeae than in other tribes. In certain species of Lathyrus, Pisum and Vavilovia, all belonging to Fabeae, the gene is either absent or a pseudogene. Our study suggests that ycf4 gene has undergone positive selection. Furthermore, the rapid evolution of the gene is locus and lineage specific and is not a shared character of the IRLC in legumes.

Plant organellar RNA editing: what 30 years of research has revealed

The central dogma in biology defines the flow of genetic information from DNA to RNA to protein. Accordingly, RNA molecules generally accurately follow the sequences of the genes from which they are transcribed. This rule is transgressed by RNA editing, which creates RNA products that differ from their DNA templates. Analyses of the RNA landscapes of terrestrial plants have indicated that RNA editing (in the form of C-U base transitions) is highly prevalent within organelles (that is, mitochondria and chloroplasts). Numerous C→U conversions (and in some plants also U→C) alter the coding sequences of many of the organellar transcripts and can also produce translatable mRNAs by creating AUG start sites or eliminating premature stop codons, or affect the RNA structure, influence splicing and alter the stability of RNAs. RNA-binding proteins are at the heart of post-transcriptional RNA expression. The C-to-U RNA editing process in plant mitochondria involves numerous nuclear-encoded factors, many of which have been identified as pentatricopeptide repeat (PPR) proteins that target editing sites in a sequence-specific manner. In this review we report on major discoveries on RNA editing in plant organelles, since it was first documented 30 years ago.

Heteroplasmy and Patterns of Cytonuclear Linkage Disequilibrium in Wild Carrot

Organellar genomes are considered to be strictly uniparentally-inherited. Uniparental inheritance allows for cytonuclear coevolution and the development of highly coordinated cytonuclear interactions. Yet, instances of biparental inheritance have been documented across eukaryotes. Biparental inheritance in otherwise uniparentally-inherited organelles is termed leakage (maternal or paternal) and allows for the presence of multiple variants of the same organellar genome within an individual, called heteroplasmy. It is unclear what, if any, evolutionary consequences are placed on nuclear and/or organellar genomes due to heteroplasmy. One way of accessing cytonuclear interactions and potential coevolution is through calculating cytonuclear linkage disequilibrium (cnLD), or the non-random association of alleles between nuclear and organellar genomes. Patterns of cnLD can indicate positive or negative cytonuclear selection, coevolution between the nuclear and organellar genomes, non-traditional organellar inheritance, or instances of ancestral heteroplasmy. In plants, cytonuclear interactions have been shown to play a role in cytoplasmic male sterility which occurs in gynodioecious species and is associated with leakage. We used the gynodioecious species, Daucus carota L. spp. carota, or wild carrot, to investigate cnLD. We genotyped a total of 265 individuals from two regions of the USA at 15 nuclear microsatellites, the mitochondrial genes cox1 and atp9, and an intergenic region between trnS and trnG (StoG) in the plastid genome to calculate nuclear–nuclear LD (nucLD), cnLD, and organellar LD (i.e., within the mtDNA and between mtDNA and ptDNA) within the two regions. We were further able to identify cox1 and StoG heteroplasmy and calculate some of the same LD measures within heteroplasmic and homoplasmic (non-heteroplasmic) datasets. We used a Z-transformation test to demonstrate that heteroplasmic individuals display significantly higher levels of cnLD within both regions. In spite of this, within and between organellar LD is low to moderate. Given these patterns of LD in two regions of the USA in which gene flow has been shown to occur between crop and wild carrot, we suggest that heteroplasmy is an evolutionary mechanism which permits the maintenance of cnLD while also acting to disrupt organellar LD.

Regulated chloroplast transcription termination

Transcription termination by the RNA polymerase (RNAP) is a fundamental step of gene expression that involves the release of the nascent transcript and dissociation of the RNAP from the DNA template. However, the functional importance of termination extends beyond the mere definition of the gene borders. Chloroplasts originate from cyanobacteria and possess their own gene expression system. Plastids have a unique hybrid transcription system consisting of two different types of RNAPs of dissimilar phylogenetic origin together with several additional nuclear encoded components. Although the basic components involved in chloroplast transcription have been identified, little attention has been paid to the chloroplast transcription termination. Recent identification and functional characterization of novel factors in regulating transcription termination in Arabidopsis chloroplasts via genetic and biochemical approaches have provided insights into the mechanisms and significance of transcription termination in chloroplast gene expression. This review provides an overview of the current knowledge of the transcription termination in chloroplasts.

Cytonuclear integration and co-evolution

The partitioning of genetic material between the nucleus and cytoplasmic (mitochondrial and plastid) genomes within eukaryotic cells necessitates coordinated integration between these genomic compartments, with important evolutionary and biomedical implications. Classic questions persist about the pervasive reduction of cytoplasmic genomes via a combination of gene loss, transfer and functional replacement - and yet why they are almost always retained in some minimal form. One striking consequence of cytonuclear integration is the existence of 'chimeric' enzyme complexes composed of subunits encoded in two different genomes. Advances in structural biology and comparative genomics are yielding important insights into the evolution of such complexes, including correlated sequence changes and recruitment of novel subunits. Thus, chimeric cytonuclear complexes provide a powerful window into the mechanisms of molecular co-evolution.

Comparative analysis of the complete plastomes of garlic Allium sativum and bulb onion Allium cepa

Sequencing and comparative characterization of plant plastid genomes, or plastomes, is an important tool for modern phylogenetic and taxonomic studies, as well as for understanding the plastome evolution. The genusAlliumL. (family Amaryllidaceae) incorporates more than 900 species, includes economically signifi­cant vegetable crops such as garlicA. sativum, onionA. cepa, leekA. porrum, etc. In this work, the plastome of garlicA. sativumhas been completely sequenced. TheA. sativumplastome is 153172 bp in size. It consists of a large unique (LSC, 82035 bp) and small unique (SSC, 18015 bp) copies, separated by inverted repeats (IRa and IRb) of 26561 bp each. In the garlic plastome, 134 genes have been annotated: 82 protein-coding genes, 38 tRNA genes, 8 rRNA genes, and 6 pseudogenes. Comparative analysis ofA. sativumandA. cepaplastomes reveals differences in the sizes of structural elements and spacers at the inverted repeat bound­aries. The total numbers of genes inA. sativumandA. cepaare the same, but the gene composition is dif­ferent: therpl22gene is functional inA. sativum, being a pseudogene inA. cepa; conversely, therps16gene is a pseudogene inA. sativumand a protein-coding gene inA. cepa. In theA. sativumandA. cepaplastomes, 32 SSR sequences have been identified. More than half of them are dinucleotides, and the remaining are tetra-, penta-, and hexanucleotides at the same time, trinucleotides were absent. The compared plastomes differ in the numbers of certain SSRs, and some are present in only one of the species.

Some theoretical insights into the hologenome theory of evolution and the role of microbes in speciation

Research on symbiotic communities (microbiomes) of multicellular organisms seems to be changing our understanding of how species of plants and animals have evolved over millions of years. The quintessence of these discoveries is the emergence of the hologenome theory of evolution, founded on the concept that a holobiont (a host along with all of its associated symbiotic microorganisms) acts a single unit of selection in the process of evolution. Although the hologenome theory has become very popular among certain scientific circles, its principles are still being debated. In this paper, we argue, firstly, that only a very small number of symbiotic microorganisms are sufficiently integrated into multicellular organisms to act in concert with them as units of selection, thus rendering claims that holobionts are units of selection invalid. Secondly, even though holobionts are not units of selection, they can still constitute genuine units from an evolutionary perspective, provided we accept certain constraints: mainly, they should be considered units of co-operation. Thirdly, we propose a reconciliation of the role of symbiotic microorganisms with the theory of speciation through the use of a developed framework. Mainly, we will argue that, in order to understand the role of microorganisms in the speciation of multicellular organisms, it is not necessary to consider holobionts units of selection; it is sufficient to consider them units of co-operation.

Contrasting Patterns of Nucleotide Substitution Rates Provide Insight into Dynamic Evolution of Plastid and Mitochondrial Genomes of Geranium

Geraniaceae have emerged as a model system for investigating the causes and consequences of variation in plastid and mitochondrial genomes. Incredible structural variation in plastid genomes (plastomes) and highly accelerated evolutionary rates have been reported in selected lineages and functional groups of genes in both plastomes and mitochondrial genomes (mitogenomes), and these phenomena have been implicated in cytonuclear incompatibility. Previous organelle genome studies have included limited sampling of Geranium, the largest genus in the family with over 400 species. This study reports on rates and patterns of nucleotide substitutions in plastomes and mitogenomes of 17 species of Geranium and representatives of other Geraniaceae. As detected across other angiosperms, substitution rates in the plastome are 3.5 times higher than the mitogenome in most Geranium. However, in the branch leading to Geranium brycei/Geranium incanum mitochondrial genes experienced significantly higher dN and dS than plastid genes, a pattern that has only been detected in one other angiosperm. Furthermore, rate accelerations differ in the two organelle genomes with plastomes having increased dN and mitogenomes with increased dS. In the Geranium phaeum/Geranium reflexum clade, duplicate copies of clpP and rpoA genes that experienced asymmetric rate divergence were detected in the single copy region of the plastome. In the case of rpoA, the branch leading to G. phaeum/G. reflexum experienced positive selection or relaxation of purifying selection. Finally, the evolution of acetyl-CoA carboxylase is unusual in Geraniaceae because it is only the second angiosperm family where both prokaryotic and eukaryotic ACCases functionally coexist in the plastid.

Nuclear-mitochondrial epistasis: a gene's eye view of genomic conflict

We use population genetic models to investigate the cooperative and conflicting synergistic fitness effects between genes from the nucleus and the mitochondrion. By varying fitness parameters, we examine the scope for conflict relative to cooperation among genomes and the utility of the “gene's eye view” analytical approach, which is based on the marginal average fitness of specific alleles. Because sexual conflict can maintain polymorphism of mitochondrial haplotypes, we can explore two types of evolutionary conflict (genomic and sexual) with one epistatic model. We find that the nuclear genetic architecture (autosomal, X‐linked, or Z‐linked) and the mating system change the regions of parameter space corresponding to the evolution by sexual and genomic conflict. For all models, regardless of conflict or cooperation, we find that population mean fitness increases monotonically as evolution proceeds. Moreover, we find that the process of gene frequency change with positive, synergistic fitnesses is self‐accelerating, as the success of an allele in one genome or in one sex increases the frequency of the interacting allele upon which its success depends. This results in runaway evolutionary dynamics caused by the positive intergenomic associations generated by selection. An inbreeding mating system tends to further accelerate these runaway dynamics because it maintains favorable host–symbiont or male–female gene combinations. In contrast, where conflict predominates, the success of an allele in one genome or in one sex diminishes the frequency of the corresponding allele in the other, resulting in considerably slower evolutionary dynamics. The rate of change of mean fitness is also much faster with positive, synergistic fitnesses and much slower where conflict is predominant. Consequently, selection rapidly fixes cooperative gene combinations, while leaving behind a slowing evolving residue of conflicting gene combinations at mutation–selection balance. We discuss how an emphasis on marginal fitness averages may obscure the interdependence of allelic fitness across genomes, making the evolutionary trajectories appear independent of one another when they are not.

Mitonuclear linkage disequilibrium in human populations

There is extensive evidence from model systems that disrupting associations between co-adapted mitochondrial and nuclear genotypes can lead to deleterious and even lethal consequences. While it is tempting to extrapolate from these observations and make inferences about the human-health effects of altering mitonuclear associations, the importance of such associations may vary greatly among species, depending on population genetics, demographic history and other factors. Remarkably, despite the extensive study of human population genetics, the statistical associations between nuclear and mitochondrial alleles remain largely uninvestigated. We analysed published population genomic data to test for signatures of historical selection to maintain mitonuclear associations, particularly those involving nuclear genes that encode mitochondrial-localized proteins (N-mt genes). We found that significant mitonuclear linkage disequilibrium (LD) exists throughout the human genome, but these associations were generally weak, which is consistent with the paucity of population genetic structure in humans. Although mitonuclear LD varied among genomic regions (with especially high levels on the X chromosome), N-mt genes were statistically indistinguishable from background levels, suggesting that selection on mitonuclear epistasis has not preferentially maintained associations involving this set of loci at a species-wide level. We discuss these findings in the context of the ongoing debate over mitochondrial replacement therapy.

Plastid-Nuclear Interaction and Accelerated Coevolution in Plastid Ribosomal Genes in Geraniaceae

Plastids and mitochondria have many protein complexes that include subunits encoded by organelle and nuclear genomes. In animal cells, compensatory evolution between mitochondrial and nuclear-encoded subunits was identified and the high mitochondrial mutation rates were hypothesized to drive compensatory evolution in nuclear genomes. In plant cells, compensatory evolution between plastid and nucleus has rarely been investigated in a phylogenetic framework. To investigate plastid–nuclear coevolution, we focused on plastid ribosomal protein genes that are encoded by plastid and nuclear genomes from 27 Geraniales species. Substitution rates were compared for five sets of genes representing plastid- and nuclear-encoded ribosomal subunit proteins targeted to the cytosol or the plastid as well as nonribosomal protein controls. We found that nonsynonymous substitution rates (dN) and the ratios of nonsynonymous to synonymous substitution rates (ω) were accelerated in both plastid- (CpRP) and nuclear-encoded subunits (NuCpRP) of the plastid ribosome relative to control sequences. Our analyses revealed strong signals of cytonuclear coevolution between plastid- and nuclear-encoded subunits, in which nonsynonymous substitutions in CpRP and NuCpRP tend to occur along the same branches in the Geraniaceae phylogeny. This coevolution pattern cannot be explained by physical interaction between amino acid residues. The forces driving accelerated coevolution varied with cellular compartment of the sequence. Increased ω in CpRP was mainly due to intensified positive selection whereas increased ω in NuCpRP was caused by relaxed purifying selection. In addition, the many indels identified in plastid rRNA genes in Geraniaceae may have contributed to changes in plastid subunits.

Drastic reduction of plastome size in the mycoheterotrophic Thismia tentaculata relative to that of its autotrophic relative Tacca chantrieri

PREMISE OF THE STUDY

Heterotrophic angiosperms tend to have reduced plastome sizes relative to those of their autotrophic relatives because genes that code for proteins involved in photosynthesis are lost. However, some plastid-encoded proteins may have vital nonphotosynthetic functions, and the plastome therefore may be retained after the loss of photosynthesis.

METHODS

We sequenced the plastome of the mycoheterotrophic species Thismia tentaculata and a representative of its sister genus, Tacca chantrieri, using next-generation technology, and we compared sequences and structures of genes and genomes of these species.

KEY RESULTS

The plastome of Tacca chantrieri is similar to those of other autotrophic taxa of Dioscoreaceae, except in a few local rearrangements and one gene loss. The plastome of Thismia tentaculata is ca. 16 kbp long with a quadripartite structure and is among the smallest known plastomes. Synteny is minimal between the plastomes of Tacca chantrieri and Thismia tentaculata. The latter includes only 12 candidate genes, with all except accD involved in protein synthesis. Of the 12 genes, trnE, trnfM, and accD are frequently among the few that remain in depauperate plastomes.

CONCLUSIONS

The plastome of Thismia tentaculata, like those of most other heterotrophic plants, includes a small number of genes previously suggested to be essential to plastome survival.

The potential role of sexual conflict and sexual selection in shaping the genomic distribution of Mito-nuclear genes

Mitochondrial interactions with the nuclear genome represent one of life's most important co-evolved mutualisms. In many organisms, mitochondria are maternally inherited, and in these cases, co-transmission between the mitochondrial and nuclear genes differs across different parts of the nuclear genome, with genes on the X chromosome having two-third probability of co-transmission, compared with one-half for genes on autosomes. These asymmetrical inheritance patterns of mitochondria and different parts of the nuclear genome have the potential to put certain gene combinations in inter-genomic co-adaptation or conflict. Previous work in mammals found strong evidence that the X chromosome has a dearth of genes that interact with the mitochondria (mito-nuclear genes), suggesting that inter-genomic conflict might drive genes off the X onto the autosomes for their male-beneficial effects. Here, we developed this idea to test coadaptation and conflict between mito-nuclear gene combinations across phylogenetically independent sex chromosomes on a far broader scale. We found that, in addition to therian mammals, only Caenorhabditis elegans showed an under-representation of mito-nuclear genes on the sex chromosomes. The remaining species studied showed no overall bias in their distribution of mito-nuclear genes. We discuss possible factors other than inter-genomic conflict that might drive the genomic distribution of mito-nuclear genes.

Using plants to elucidate the mechanisms of cytonuclear co-evolution

The presence of both cytoplasmic and nuclear genomes within eukaryotic cells raises fascinating questions about co-evolution between genomic compartments that experience fundamentally different mutation rates and modes of inheritance. The highly mutagenic environments found in the mitochondria of some eukaryotes have generated interest in the role that mitochondrial mutation accumulation plays in phenomena such as intracellular gene transfer, compensatory evolution in the nucleus and the evolution of reproductive isolation. Although plant systems have played an important historical role in the study of cytonuclear co-evolution, they remain underutilized in many respects. In particular, the enormous natural variation in DNA substitution rates, gene content and genome architecture in plant mitochondria - much of which has even been found within a single genus – provides opportunities to resolve longstanding evolutionary questions about the consequences of mitochondrial mutation accumulation. This review summarizes some of the classic questions about cytonuclear co-evolution that could be addressed by taking advantage of the variation in plants and highlights a recent analysis of the effect of mitochondrial mutation accumulation on rates of molecular evolution in the nucleus.

Gene interactions in the evolution of genomic imprinting

Numerous evolutionary theories have been developed to explain the epigenetic phenomenon of genomic imprinting. Here, we explore a subset of theories wherein non-additive genetic interactions can favour imprinting. In the simplest genic interaction--the case of underdominance--imprinting can be favoured to hide effectively low-fitness heterozygous genotypes; however, as there is no asymmetry between maternally and paternally inherited alleles in this model, other means of enforcing monoallelic expression may be more plausible evolutionary outcomes than genomic imprinting. By contrast, more successful interaction models of imprinting rely on an asymmetry between the maternally and paternally inherited alleles at a locus that favours the silencing of one allele as a means of coordinating the expression of high-fitness allelic combinations. For example, with interactions between autosomal loci, imprinting functionally preserves high-fitness genotypes that were favoured by selection in the previous generation. In this scenario, once a focal locus becomes imprinted, selection at interacting loci favours a matching imprint. Uniparental transmission generates similar asymmetries for sex chromosomes and cytoplasmic factors interacting with autosomal loci, with selection favouring the expression of either maternal or paternally derived autosomal alleles depending on the pattern of transmission of the uniparentally inherited factor. In a final class of models, asymmetries arise when genes expressed in offspring interact with genes expressed in one of its parents. Under such a scenario, a locus evolves to have imprinted expression in offspring to coordinate the interaction with its parent's genome. We illustrate these models and explore key links and differences using a unified framework.

Cytonuclear interactions and relaxed selection accelerate sequence evolution in organelle ribosomes

Many mitochondrial and plastid protein complexes contain subunits that are encoded in different genomes. In animals, nuclear-encoded mitochondrial proteins often exhibit rapid sequence evolution, which has been hypothesized to result from selection for mutations that compensate for changes in interacting subunits encoded in mutation-prone animal mitochondrial DNA. To test this hypothesis, we analyzed nuclear genes encoding cytosolic and organelle ribosomal proteins in flowering plants. The model angiosperm genus Arabidopsis exhibits low organelle mutation rates, typical of most plants. Nevertheless, we found that (nuclear-encoded) subunits of organelle ribosomes in Arabidopsis have higher amino acid sequence polymorphism and divergence than their counterparts in cytosolic ribosomes, suggesting that organelle ribosomes experience relaxed functional constraint. However, the observed difference between organelle and cytosolic ribosomes was smaller than in animals and could be partially attributed to rapid evolution in N-terminal organelle-targeting peptides that are not involved in ribosome function. To test the role of organelle mutation more directly, we used transcriptomic data from an angiosperm genus (Silene) with highly variable rates of organelle genome evolution. We found that Silene species with unusually fast-evolving mitochondrial and plastid DNA exhibited increased amino acid sequence divergence in ribosomal proteins targeted to the organelles but not in those that function in cytosolic ribosomes. Overall, these findings support the hypothesis that rapid organelle genome evolution has selected for compensatory mutations in nuclear-encoded proteins. We conclude that coevolution between interacting subunits encoded in different genomic compartments within the eukaryotic cell is an important determinant of variation in rates of protein sequence evolution.

Patterns of cyto-nuclear linkage disequilibrium in Silene latifolia: genomic heterogeneity and temporal stability

Non-random association of alleles in the nucleus and cytoplasmic organelles, or cyto-nuclear linkage disequilibrium (LD), is both an important component of a number of evolutionary processes and a statistical indicator of others. The evolutionary significance of cyto-nuclear LD will depend on both its magnitude and how stable those associations are through time. Here, we use a longitudinal population genetic data set to explore the magnitude and temporal dynamics of cyto-nuclear disequilibria through time. We genotyped 135 and 170 individuals from 16 and 17 patches of the plant species Silene latifolia in Southwestern VA, sampled in 1993 and 2008, respectively. Individuals were genotyped at 14 highly polymorphic microsatellite markers and a single-nucleotide polymorphism (SNP) in the mitochondrial gene, atp1. Normalized LD (D') between nuclear and cytoplasmic loci varied considerably depending on which nuclear locus was considered (ranging from 0.005-0.632). Four of the 14 cyto-nuclear associations showed a statistically significant shift over approximately seven generations. However, the overall magnitude of this disequilibrium was largely stable over time. The observed origin and stability of cyto-nuclear LD is most likely caused by the slow admixture between anciently diverged lineages within the species' newly invaded range, and the local spatial structure and metapopulation dynamics that are known to structure genetic variation in this system.

Lineage-specific fragmentation and nuclear relocation of the mitochondrial cox2 gene in chlorophycean green algae (Chlorophyta)

In most eukaryotes the subunit 2 of cytochrome c oxidase (COX2) is encoded in intact mitochondrial genes. Some green algae, however, exhibit split cox2 genes (cox2a and cox2b) encoding two polypeptides (COX2A and COX2B) that form a heterodimeric COX2 subunit. Here, we analyzed the distribution of intact and split cox2 gene sequences in 39 phylogenetically diverse green algae in phylum Chlorophyta obtained from databases (28 sequences from 22 taxa) and from new cox2 data generated in this work (23 sequences from 18 taxa). Our results support previous observations based on a smaller number of taxa, indicating that algae in classes Prasinophyceae, Ulvophyceae, and Trebouxiophyceae contain orthodox, intact mitochondrial cox2 genes. In contrast, all of the algae in Chlorophyceae that we examined exhibited split cox2 genes, and could be separated into two groups: one that has a mitochondrion-localized cox2a gene and a nucleus-localized cox2b gene ("Scenedesmus-like"), and another that has both cox2a and cox2b genes in the nucleus ("Chlamydomonas-like"). The location of the split cox2a and cox2b genes was inferred using five different criteria: differences in amino acid sequences, codon usage (mitochondrial vs. nuclear), codon preference (third position frequencies), presence of nucleotide sequences encoding mitochondrial targeting sequences and presence of spliceosomal introns. Distinct green algae could be grouped according to the form of cox2 gene they contain: intact or fragmented, mitochondrion- or nucleus-localized, and intron-containing or intron-less. We present a model describing the events that led to mitochondrial cox2 gene fragmentation and the independent and sequential migration of cox2a and cox2b genes to the nucleus in chlorophycean green algae. We also suggest that the distribution of the different forms of the cox2 gene provides important insights into the phylogenetic relationships among major groups of Chlorophyceae.

Evidence of a paucity of genes that interact with the mitochondrion on the X in mammals

Mitochondria are essential organelles whose replication, development, and physiology are dependent upon coordinated gene interactions with both the mitochondrial and the nuclear genomes. The evolution of coadapted (CA) nuclear-mitochondrial gene combinations would be facilitated if such nuclear genes were located on the X-chromosome instead of on the autosomes because of the increased probability of cotransmission. Here, we test the prediction of the CA hypothesis by investigating the chromosomal distribution of nuclear genes that interact with mitochondria. Using the online genome database BIOMART, we compared the density of genes that have a mitochondrion cellular component annotation across chromosomes in 16 vertebrates. We find a strong and highly significant genomic pattern against the CA hypothesis: nuclear genes interacting with the mitochondrion are significantly underrepresented on the X-chromosome in mammals but not in birds. We interpret our findings in terms of sexual conflict as a mechanism that may generate the observed pattern. Our finding extends single-gene theory for the evolution of sexually antagonistic genes to nuclear-mitochondrial gene combinations.

The transcription machineries of plant mitochondria and chloroplasts: Composition, function, and regulation

Although genomes of mitochondria and plastids are very small compared to those of their bacterial ancestors, the transcription machineries of these organelles are of surprising complexity. With respect to the number of different RNA polymerases per organelle, the extremes are represented on one hand by chloroplasts of eudicots which use one bacterial-type RNA polymerase and two phage-type RNA polymerases to transcribe their genes, and on the other hand by Physcomitrella possessing three mitochondrial RNA polymerases of the phage type. Transcription of genes/operons is often driven by multiple promoters in both organelles. This review describes the principle components of the transcription machineries (RNA polymerases, transcription factors, promoters) and the division of labor between the different RNA polymerases. While regulation of transcription in mitochondria seems to be only of limited importance, the plastid genes of higher plants respond to exogenous and endogenous cues rather individually by altering their transcriptional activities.

Horizontal transmission rapidly erodes disequilibria between organelle and symbiont genomes

We investigate the generation and decay of interspecific disequilibrium (ID) between organelle and symbiont genomes as a function of the rate of horizontal transmission. We show that rare horizontal transmission greatly diminishes the covariance between organelle and symbiont genomes. This result has two important implications. First, a low level of ID does not indicate low levels of vertical transmission. Second, even with low levels of horizontal transmission, the additive effects of host and symbiont loci will determine the response to selection, while epistatic effects will not be selectable.

Localized hypermutation and associated gene losses in legume chloroplast genomes

Point mutations result from errors made during DNA replication or repair, so they are usually expected to be homogeneous across all regions of a genome. However, we have found a region of chloroplast DNA in plants related to sweetpea (Lathyrus) whose local point mutation rate is at least 20 times higher than elsewhere in the same molecule. There are very few precedents for such heterogeneity in any genome, and we suspect that the hypermutable region may be subject to an unusual process such as repeated DNA breakage and repair. The region is 1.5 kb long and coincides with a gene, ycf4, whose rate of evolution has increased dramatically. The product of ycf4, a photosystem I assembly protein, is more divergent within the single genus Lathyrus than between cyanobacteria and other angiosperms. Moreover, ycf4 has been lost from the chloroplast genome in Lathyrus odoratus and separately in three other groups of legumes. Each of the four consecutive genes ycf4-psaI-accD-rps16 has been lost in at least one member of the legume "inverted repeat loss" clade, despite the rarity of chloroplast gene losses in angiosperms. We established that accD has relocated to the nucleus in Trifolium species, but were unable to find nuclear copies of ycf4 or psaI in Lathyrus. Our results suggest that, as well as accelerating sequence evolution, localized hypermutation has contributed to the phenomenon of gene loss or relocation to the nucleus.