Dynamics of reductive genome evolution in mitochondria and obligate intracellular microbes

A gene-rich mitochondrion with a unique ancestral protein transport system

Mitochondria originated from an ancient endosymbiosis involving an alphaproteobacterium.1,2,3 Over time, these organelles reduced their gene content massively, with most genes being transferred to the host nucleus before the last eukaryotic common ancestor (LECA).4 This process has yielded varying gene compositions in modern mitogenomes, including the complete loss of this organellar genome in some extreme cases.5,6,7,8,9,10,11,12,13,14 At the other end of the spectrum, jakobids harbor the most gene-rich mitogenomes, encoding 60-66 proteins.8 Here, we introduce the mitogenome of Mantamonas sphyraenae, a protist from the deep-branching CRuMs supergroup.15,16 Remarkably, it boasts the most gene-rich mitogenome outside of jakobids, by housing 91 genes, including 62 protein-coding ones. These include rare homologs of the four subunits of the bacterial-type cytochrome c maturation system I (CcmA, CcmB, CcmC, and CcmF) alongside a unique ribosomal protein S6. During the early evolution of mitochondria, gene transfer from the proto-mitochondrial endosymbiont to the nucleus became possible thanks to systems facilitating the transport of proteins synthesized in the host cytoplasm back to the mitochondrion. In addition to the universally found eukaryotic protein import systems, jakobid mitogenomes were reported to uniquely encode the SecY transmembrane protein of the Sec general secretory pathway, whose evolutionary origin was however unclear. The Mantamonas mitogenome not only encodes SecY but also SecA, SecE, and SecG, making it the sole eukaryote known to house a complete mitochondrial Sec translocation system. Furthermore, our phylogenetic and comparative genomic analyses provide compelling evidence for the alphaproteobacterial origin of this system, establishing its presence in LECA.

Multifaceted Dinoflagellates and the Marine Model Prorocentrum cordatum

BACKGROUND

Dinoflagellates are a monophyletic group within the taxon Alveolata, which comprises unicellular eukaryotes. Dinoflagellates have long been studied for their organismic and morphologic diversity as well as striking cellular features. They have a main size range of 10-100 µm, a complex "cell covering", exceptionally large genomes (∼1-250 Gbp with a mean of 50,000 protein-encoding genes) spread over a variable number of highly condensed chromosomes, and perform a closed mitosis with extranuclear spindles (dinomitosis). Photosynthetic, marine, and free-living Prorocentrum cordatum is a ubiquitously occurring, bloom-forming dinoflagellate, and an emerging model system, particularly with respect to systems biology.

SUMMARY

Focused ion beam/scanning electron microscopy (FIB/SEM) analysis of P. cordatum recently revealed (i) a flattened nucleus with unusual structural features and a total of 62 tightly packed chromosomes, (ii) a single, barrel-shaped chloroplast devoid of grana and harboring multiple starch granules, (iii) a single, highly reticular mitochondrion, and (iv) multiple phosphate and lipid storage bodies. Comprehensive proteomics of subcellular fractions suggested (i) major basic nuclear proteins to participate in chromosome condensation, (ii) composition of nuclear pores to differ from standard knowledge, (iii) photosystems I and II, chloroplast complex I, and chlorophyll a-b binding light-harvesting complex to form a large megacomplex (>1.5 MDa), and (iv) an extraordinary richness in pigment-binding proteins. Systems biology-level investigation of heat stress response demonstrated a concerted down-regulation of CO2-concentrating mechanisms, CO2-fixation, central metabolism, and monomer biosynthesis, which agrees with reduced growth yields.

KEY MESSAGES

FIB/SEM analysis revealed new insights into the remarkable subcellular architecture of P. cordatum, complemented by proteogenomic unraveling of novel nuclear structures and a photosynthetic megacomplex. These recent findings are put in the wider context of current understanding of dinoflagellates.

Endosymbioses Have Shaped the Evolution of Biological Diversity and Complexity Time and Time Again

Life on Earth comprises prokaryotes and a broad assemblage of endosymbioses. The pages of Molecular Biology and Evolution and Genome Biology and Evolution have provided an essential window into how these endosymbiotic interactions have evolved and shaped biological diversity. Here, we provide a current perspective on this knowledge by drawing on decades of revelatory research published in Molecular Biology and Evolution and Genome Biology and Evolution, and insights from the field at large. The accumulated work illustrates how endosymbioses provide hosts with novel phenotypes that allow them to transition between adaptive landscapes to access environmental resources. Such endosymbiotic relationships have shaped and reshaped life on Earth. The early serial establishment of mitochondria and chloroplasts through endosymbioses permitted massive upscaling of cellular energetics, multicellularity, and terrestrial planetary greening. These endosymbioses are also the foundation upon which all later ones are built, including everything from land-plant endosymbioses with fungi and bacteria to nutritional endosymbioses found in invertebrate animals. Common evolutionary mechanisms have shaped this broad range of interactions. Endosymbionts generally experience adaptive and stochastic genome streamlining, the extent of which depends on several key factors (e.g. mode of transmission). Hosts, in contrast, adapt complex mechanisms of resource exchange, cellular integration and regulation, and genetic support mechanisms to prop up degraded symbionts. However, there are significant differences between endosymbiotic interactions not only in how partners have evolved with each other but also in the scope of their influence on biological diversity. These differences are important considerations for predicting how endosymbioses will persist and adapt to a changing planet.

Comparative Mitogenomics Reveals Cryptic Species in Sillago ingenuua McKay, 1985 (Perciformes: Sillaginidae)

It is unreliable to identify marine fishes only by external morphological features. Species misidentification brings great challenges to fishery research, resource monitoring and ecomanagement. Sillago ingenuua is an important part of commercial marine fishes, and in which, the morphological differences between different groups are not obvious. Here, we compared different geographical groups of S. ingenuua which were collected from Xiamen, Dongshan, Keelung, Songkhla and Java. The results showed that all samples of S. ingenuua were similar in external morphological characteristics and the shape of the swim bladder, but there were two distinctive lineages which were flagged as cryptic species based on DNA barcoding. The comparative mitogenomic results showed that S. ingenuua A and S. ingenuua B were identical in structural organization and gene arrangement. Their nucleotide composition and codon usage were also similar. A phylogenetic analysis was performed based on 13 concatenated PCGs from eight Sillago species. The results showed that the genetic distance between S. ingenuua A and S. ingenuua B was large (D = 0.069), and this genetic distance was large enough to reveal that S. ingenuua A and S. ingenuua B might be different species.

Microbiota in disease-transmitting vectors

Haematophagous arthropods, including mosquitoes, ticks, flies, triatomine bugs and lice (here referred to as vectors), are involved in the transmission of various pathogens to mammals on whom they blood feed. The diseases caused by these pathogens, collectively known as vector-borne diseases (VBDs), threaten the health of humans and animals. Although the vector arthropods differ in life histories, feeding behaviour as well as reproductive strategies, they all harbour symbiotic microorganisms, known as microbiota, on which they depend for completing essential aspects of their biology, such as development and reproduction. In this Review, we summarize the shared and unique key features of the symbiotic associations that have been characterized in the major vector taxa. We discuss the crosstalks between microbiota and their arthropod hosts that influence vector metabolism and immune responses relevant for pathogen transmission success, known as vector competence. Finally, we highlight how current knowledge on symbiotic associations is being explored to develop non-chemical-based alternative control methods that aim to reduce vector populations, or reduce vector competence. We conclude by highlighting the remaining knowledge gaps that stand to advance basic and translational aspects of vector-microbiota interactions.

Single-Cell Genomics Reveals the Divergent Mitochondrial Genomes of Retaria (Foraminifera and Radiolaria)

Mitochondria originated from an ancient bacterial endosymbiont that underwent reductive evolution by gene loss and endosymbiont gene transfer to the nuclear genome. The diversity of mitochondrial genomes published to date has revealed that gene loss and transfer processes are ongoing in many lineages. Most well-studied eukaryotic lineages are represented in mitochondrial genome databases, except for the superphylum Retaria-the lineage comprising Foraminifera and Radiolaria. Using single-cell approaches, we determined two complete mitochondrial genomes of Foraminifera and two nearly complete mitochondrial genomes of radiolarians. We report the complete coding content of an additional 14 foram species. We show that foraminiferan and radiolarian mitochondrial genomes contain a nearly fully overlapping but reduced mitochondrial gene complement compared to other sequenced rhizarians. In contrast to animals and fungi, many protists encode a diverse set of proteins on their mitochondrial genomes, including several ribosomal genes; however, some aerobic eukaryotic lineages (euglenids, myzozoans, and chlamydomonas-like algae) have reduced mitochondrial gene content and lack all ribosomal genes. Similar to these reduced outliers, we show that retarian mitochondrial genomes lack ribosomal protein and tRNA genes, contain truncated and divergent small and large rRNA genes, and contain only 14 or 15 protein-coding genes, including nad1, -3, -4, -4L, -5, and -7, cob, cox1, -2, and -3, and atp1, -6, and -9, with forams and radiolarians additionally carrying nad2 and nad6, respectively. In radiolarian mitogenomes, a noncanonical genetic code was identified in which all three stop codons encode amino acids. Collectively, these results add to our understanding of mitochondrial genome evolution and fill in one of the last major gaps in mitochondrial sequence databases. IMPORTANCE We present the reduced mitochondrial genomes of Retaria, the rhizarian lineage comprising the phyla Foraminifera and Radiolaria. By applying single-cell genomic approaches, we found that foraminiferan and radiolarian mitochondrial genomes contain an overlapping but reduced mitochondrial gene complement compared to other sequenced rhizarians. An alternative genetic code was identified in radiolarian mitogenomes in which all three stop codons encode amino acids. Collectively, these results shed light on the divergent nature of the mitochondrial genomes from an ecologically important group, warranting further questions into the biological underpinnings of gene content variability and genetic code variation between mitochondrial genomes.

Mitochondrial genomes in the iconic reindeer lichens: Architecture, variation, and synteny across multiple evolutionary scales

Variation in mitochondrial genome composition across intraspecific, interspecific, and higher taxonomic scales has been little studied in lichen obligate symbioses. Cladonia is one of the most diverse and ecologically important lichen genera, with over 500 species representing an array of unique morphologies and chemical profiles. Here, we assess mitochondrial genome diversity and variation in this flagship genus, with focused sampling of two clades of the "true" reindeer lichens, Cladonia subgenus Cladina, and additional genomes from nine outgroup taxa. We describe composition and architecture at the gene and the genome scale, examining patterns in organellar genome size in larger taxonomic groups in Ascomycota. Mitochondrial genomes of Cladonia, Pilophorus, and Stereocaulon were consistently larger than those of Lepraria and contained more introns, suggesting a selective pressure in asexual morphology in Lepraria driving it toward genomic simplification. Collectively, lichen mitochondrial genomes were larger than most other fungal life strategies, reaffirming the notion that coevolutionary streamlining does not correlate to genome size reductions. Genomes from Cladonia ravenelii and Stereocaulon pileatum exhibited ATP9 duplication, bearing paralogs that may still be functional. Homing endonuclease genes (HEGs), though scarce in Lepraria, were diverse and abundant in Cladonia, exhibiting variable evolutionary histories that were sometimes independent of the mitochondrial evolutionary history. Intraspecific HEG diversity was also high, with C. rangiferina especially bearing a range of HEGs with one unique to the species. This study reveals a rich history of events that have transformed mitochondrial genomes of Cladonia and related genera, allowing future study alongside a wealth of assembled genomes.

The colonization of land was a likely driving force for the evolution of mitochondrial retrograde signalling in plants

Most retrograde signalling research in plants was performed using Arabidopsis, so an evolutionary perspective on mitochondrial retrograde regulation (MRR) is largely missing. Here, we used phylogenetics to track the evolutionary origins of factors involved in plant MRR. In all cases, the gene families can be traced to ancestral green algae or earlier. However, the specific subfamilies containing factors involved in plant MRR in many cases arose during the transition to land. NAC transcription factors with C-terminal transmembrane domains, as observed in the key regulator ANAC017, can first be observed in non-vascular mosses, and close homologs to ANAC017 can be found in seed plants. Cyclin-dependent kinases (CDKs) are common to eukaryotes, but E-type CDKs that control MRR also diverged in conjunction with plant colonization of land. AtWRKY15 can be traced to the earliest land plants, while AtWRKY40 only arose in angiosperms and AtWRKY63 even more recently in Brassicaceae. Apetala 2 (AP2) transcription factors are traceable to algae, but the ABI4 type again only appeared in seed plants. This strongly suggests that the transition to land was a major driver for developing plant MRR pathways, while additional fine-tuning events have appeared in seed plants or later. Finally, we discuss how MRR may have contributed to meeting the specific challenges that early land plants faced during terrestrialization.

RibU is an essential determinant of Listeria pathogenesis that mediates acquisition of FMN and FAD during intracellular growth

Flavin mononucleotide (FMN) and flavin adenine dinucleotide (FAD) are essential riboflavin-derived cofactors involved in a myriad of redox reactions across all forms of life. Nevertheless, the basis of flavin acquisition strategies by riboflavin auxotrophic pathogens remains poorly defined. In this study, we examined how the facultative intracellular pathogen Listeria monocytogenes, a riboflavin auxotroph, acquires flavins during infection. A L. monocytogenes mutant lacking the putative riboflavin transporter (RibU) was completely avirulent in mice but had no detectable growth defect in nutrient-rich media. However, unlike wild type, the RibU mutant was unable to grow in defined media supplemented with FMN or FAD or to replicate in macrophages starved for riboflavin. Consistent with RibU functioning to scavenge FMN and FAD inside host cells, a mutant unable to convert riboflavin to FMN or FAD retained virulence and grew in cultured macrophages and in spleens and livers of infected mice. However, this FMN- and FAD-requiring strain was unable to grow in the gallbladder or intestines, where L. monocytogenes normally grows extracellularly, suggesting that these sites do not contain sufficient flavin cofactors to promote replication. Thus, by deleting genes required to synthesize FMN and FAD, we converted L. monocytogenes from a facultative to an obligate intracellular pathogen. Collectively, these data indicate that L. monocytogenes requires riboflavin to grow extracellularly in vivo but scavenges FMN and FAD to grow in host cells.

Genomic evolution and adaptation of arthropod-associated Rickettsia

Rickettsia species are endosymbionts hosted by arthropods and are known to cause mild to fatal diseases in humans. Here, we analyse the evolution and diversity of 34 Rickettsia species using a pangenomic meta-analysis (80 genomes/41 plasmids). Phylogenomic trees showed that Rickettsia spp. diverged into two Spotted Fever groups, a Typhus group, a Canadensis group and a Bellii group, and may have inherited their plasmids from an ancestral plasmid that persisted in some strains or may have been lost by others. The results suggested that the ancestors of Rickettsia spp. might have infected Acari and/or Insecta and probably diverged by persisting inside and/or switching hosts. Pangenomic analysis revealed that the Rickettsia genus evolved through a strong interplay between genome degradation/reduction and/or expansion leading to possible distinct adaptive trajectories. The genus mainly shared evolutionary relationships with α-proteobacteria, and also with γ/β/δ-proteobacteria, cytophagia, actinobacteria, cyanobacteria, chlamydiia and viruses, suggesting lateral exchanges of several critical genes. These evolutionary processes have probably been orchestrated by an abundance of mobile genetic elements, especially in the Spotted Fever and Bellii groups. In this study, we provided a global evolutionary genomic view of the intracellular Rickettsia that may help our understanding of their diversity, adaptation and fitness.

The glycan alphabet is not universal: a hypothesis

Several monosaccharides constitute naturally occurring glycans, but it is uncertain whether they constitute a universal set like the alphabets of proteins and DNA. Based on the available experimental observations, it is hypothesized herein that the glycan alphabet is not universal. Data on the presence/absence of pathways for the biosynthesis of 55 monosaccharides in 12 939 completely sequenced archaeal and bacterial genomes are presented in support of this hypothesis. Pathways were identified by searching for homologues of biosynthesis pathway enzymes. Substantial variations were observed in the set of monosaccharides used by organisms belonging to the same phylum, genera and even species. Monosaccharides were grouped as common, less common and rare based on their prevalence in Archaea and Bacteria. It was observed that fewer enzymes are sufficient to biosynthesize monosaccharides in the common group. It appears that the common group originated before the formation of the three domains of life. In contrast, the rare group is confined to a few species in a few phyla, suggesting that these monosaccharides evolved much later. Fold conservation, as observed in aminotransferases and SDR (short-chain dehydrogenase reductase) superfamily members involved in monosaccharide biosynthesis, suggests neo- and sub-functionalization of genes led to the formation of the rare group monosaccharides. The non-universality of the glycan alphabet begets questions about the role of different monosaccharides in determining an organism’s fitness.

A Biochemical and Structural Understanding of TOM Complex Interactions and Implications for Human Health and Disease

The central role mitochondria play in cellular homeostasis has made its study critical to our understanding of various aspects of human health and disease. Mitochondria rely on the translocase of the outer membrane (TOM) complex for the bulk of mitochondrial protein import. In addition to its role as the major entry point for mitochondrial proteins, the TOM complex serves as an entry pathway for viral proteins. TOM complex subunits also participate in a host of interactions that have been studied extensively for their function in neurodegenerative diseases, cardiovascular diseases, innate immunity, cancer, metabolism, mitophagy and autophagy. Recent advances in our structural understanding of the TOM complex and the protein import machinery of the outer mitochondrial membrane have made structure-based therapeutics targeting outer mitochondrial membrane proteins during mitochondrial dysfunction an exciting prospect. Here, we describe advances in understanding the TOM complex, the interactome of the TOM complex subunits, the implications for the development of therapeutics, and our understanding of the structure/function relationship between components of the TOM complex and mitochondrial homeostasis.

HERMES: An improved method to test mitochondrial genome molecular synapomorphies among clades

Mitochondrial chromosomes have diversified among eukaryotes and many different architectures and features are now acknowledged for this genome. Here we present the improved HERMES index, which can measure and quantify the amount of molecular change experienced by mitochondrial genomes. We test the improved approach with ten molecular phylogenetic studies based on complete mitochondrial genomes, representing six bilaterian Phyla. In most cases, HERMES analysis spotted out clades or single species with peculiar molecular synapomorphies, allowing to identify phylogenetic and ecological patterns. The software presented herein handles linear, circular, and multi-chromosome genomes, thus widening the HERMES scope to the complete eukaryotic domain.

Rescuing stalled mammalian mitoribosomes - what can we learn from bacteria?

In the canonical process of translation, newly completed proteins escape from the ribosome following cleavage of the ester bond that anchors the polypeptide to the P-site tRNA, after which the ribosome can be recycled to initiate a new round of translation. Not all protein synthesis runs to completion as various factors can impede the progression of ribosomes. Rescuing of stalled ribosomes in mammalian mitochondria, however, does not share the same mechanisms that many bacteria use. The classic method for rescuing bacterial ribosomes is trans-translation. The key components of this system are absent from mammalian mitochondria; however, four members of a translation termination factor family are present, with some evidence of homology to members of a bacterial back-up rescue system. To date, there is no definitive demonstration of any other member of this family functioning in mitoribosome rescue. Here, we provide an overview of the processes and key players of canonical translation termination in both bacteria and mammalian mitochondria, followed by a perspective of the bacterial systems used to rescue stalled ribosomes. We highlight any similarities or differences with the mitochondrial translation release factors, and suggest potential roles for these proteins in ribosome rescue in mammalian mitochondria.

The functional ClpXP protease of Chlamydia trachomatis requires distinct clpP genes from separate genetic loci

Clp proteases play a central role in bacterial physiology and, for some bacterial species, are even essential for survival. Also due to their conservation among bacteria including important human pathogens, Clp proteases have recently attracted considerable attention as antibiotic targets. Here, we functionally reconstituted and characterized the ClpXP protease of Chlamydia trachomatis (ctClpXP), an obligate intracellular pathogen and the causative agent of widespread sexually transmitted diseases in humans. Our in vitro data show that ctClpXP is formed by a hetero-tetradecameric proteolytic core, composed of two distinct homologs of ClpP (ctClpP1 and ctClpP2), that associates with the unfoldase ctClpX via ctClpP2 for regulated protein degradation. Antibiotics of the ADEP class interfere with protease functions by both preventing the interaction of ctClpX with ctClpP1P2 and activating the otherwise dormant proteolytic core for unregulated proteolysis. Thus, our results reveal molecular insight into ctClpXP function, validating this protease as an antibacterial target.

Genome streamlining via complete loss of introns has occurred multiple times in lichenized fungal mitochondria

Reductions in genome size and complexity are a hallmark of obligate symbioses. The mitochondrial genome displays clear examples of these reductions, with the ancestral alpha-proteobacterial genome size and gene number having been reduced by orders of magnitude in most descendent modern mitochondrial genomes. Here, we examine patterns of mitochondrial evolution specifically looking at intron size, number, and position across 58 species from 21 genera of lichenized Ascomycete fungi, representing a broad range of fungal diversity and niches. Our results show that the cox1gene always contained the highest number of introns out of all the mitochondrial protein-coding genes, that high intron sequence similarity (>90%) can be maintained between different genera, and that lichens have undergone at least two instances of complete, genome-wide intron loss consistent with evidence for genome streamlining via loss of parasitic, noncoding DNA, in Phlyctis boliviensisand Graphis lineola. Notably, however, lichenized fungi have not only undergone intron loss but in some instances have expanded considerably in size due to intron proliferation (e.g., Alectoria fallacina and Parmotrema neotropicum), even between closely related sister species (e.g., Cladonia). These results shed light on the highly dynamic mitochondrial evolution that is occurring in lichens and suggest that these obligate symbiotic organisms are in some cases undergoing recent, broad-scale genome streamlining via loss of protein-coding genes as well as noncoding, parasitic DNA elements.

Bacterial endosymbiont Cardinium cSfur genome sequence provides insights for understanding the symbiotic relationship in Sogatella furcifera host

BACKGROUND

Sogatella furcifera is a migratory pest that damages rice plants and causes severe economic losses. Due to its ability to annually migrate long distances, S. furcifera has emerged as a major pest of rice in several Asian countries. Symbiotic relationships of inherited bacteria with terrestrial arthropods have significant implications. The genus Cardinium is present in many types of arthropods, where it influences some host characteristics. We present a report of a newly identified strain of the bacterial endosymbiont Cardinium cSfur in S. furcifera.

RESULT

From the whole genome of S. furcifera previously sequenced by our laboratory, we assembled the whole genome sequence of Cardinium cSfur. The sequence comprised 1,103,593 bp with a GC content of 39.2%. The phylogenetic tree of the Bacteroides phylum to which Cardinium cSfur belongs suggests that Cardinium cSfur is closely related to the other strains (Cardinium cBtQ1 and cEper1) that are members of the Amoebophilaceae family. Genome comparison between the host-dependent endosymbiont including Cardinium cSfur and free-living bacteria revealed that the endosymbiont has a smaller genome size and lower GC content, and has lost some genes related to metabolism because of its special environment, which is similar to the genome pattern observed in other insect symbionts. Cardinium cSfur has limited metabolic capability, which makes it less contributive to metabolic and biosynthetic processes in its host. From our findings, we inferred that, to compensate for its limited metabolic capability, Cardinium cSfur harbors a relatively high proportion of transport proteins, which might act as the hub between it and its host. With its acquisition of the whole operon related to biotin synthesis and glycolysis related genes through HGT event, Cardinium cSfur seems to be undergoing changes while establishing a symbiotic relationship with its host.

CONCLUSION

A novel bacterial endosymbiont strain (Cardinium cSfur) has been discovered. A genomic analysis of the endosymbiont in S. furcifera suggests that its genome has undergone certain changes to facilitate its settlement in the host. The envisaged potential reproduction manipulative ability of the new endosymbiont strain in its S. furcifera host has vital implications in designing eco-friendly approaches to combat the insect pest.

Unravelling the relationship between the tsetse fly and its obligate symbiont Wigglesworthia: transcriptomic and metabolomic landscapes reveal highly integrated physiological networks

Insects with restricted diets rely on obligate microbes to fulfil nutritional requirements essential for biological function. Tsetse flies, vectors of African trypanosome parasites, feed exclusively on vertebrate blood and harbour the obligate endosymbiont Wigglesworthia glossinidia. Without Wigglesworthia, tsetse are unable to reproduce. These symbionts are sheltered within specialized cells (bacteriocytes) that form the midgut-associated bacteriome organ. To decipher the core functions of this symbiosis essential for tsetse's survival, we performed dual-RNA-seq analysis of the bacteriome, coupled with metabolomic analysis of bacteriome and haemolymph collected from normal and symbiont-cured (sterile) females. Bacteriocytes produce immune regulatory peptidoglycan recognition protein (pgrp-lb) that protects Wigglesworthia, and a multivitamin transporter (smvt) that can aid in nutrient dissemination. Wigglesworthia overexpress a molecular chaperone (GroEL) to augment their translational/transport machinery and biosynthesize an abundance of B vitamins (specifically B₁-, B₂-, B₃- and B₆-associated metabolites) to supplement the host's nutritionally deficient diet. The absence of Wigglesworthia's contributions disrupts multiple metabolic pathways impacting carbohydrate and amino acid metabolism. These disruptions affect the dependent downstream processes of nucleotide biosynthesis and metabolism and biosynthesis of S-adenosyl methionine (SAM), an essential cofactor. This holistic fundamental knowledge of the symbiotic dialogue highlights new biological targets for the development of innovative vector control methods.

SmithRNAs: Could Mitochondria "Bend" Nuclear Regulation?

Typically, animal mitochondria have very compact genomes, with few short intergenic regions, and no introns. Hence, it may seem that there is little space for unknown functions in mitochondrial DNA (mtDNA). However, mtDNA can also operate through RNA interference, as small non coding RNAs (sncRNAs) produced by mtDNA have already been proposed for humans. We sequenced sncRNA libraries from isolated mitochondria of Ruditapes philippinarum (Mollusca Bivalvia) gonads, a species with doubly uniparental inheritance of mitochondria, and identified several putative sncRNAs of mitochondrial origin. Some sncRNAs are transcribed by intergenic regions that form stable stem-hairpin structures, which makes them good miRNA-like candidates. We decided to name them small mitochondrial highly-transcribed RNAs (smithRNAs). Many concurrent data support that we have recovered sncRNAs of mitochondrial origin that might be involved in gonad formation and able to affect nuclear gene expression. This possibility has been never suggested before. If mtDNA can affect nuclear gene expression through RNA interference, this opens a plethora of new possibilities for it to interact with the nucleus, and makes metazoan mtDNA a much more complex genome than previously thought.

Burrowers from the Past: Mitochondrial Signatures of Ordovician Bivalve Infaunalization

Bivalves and gastropods are the two largest classes of extant molluscs. Despite sharing a huge number of features, they do not share a key ecological one: gastropods are essentially epibenthic, although most bivalves are infaunal. However, this is not the ancestral bivalve condition; Cambrian forms were surface crawlers and only during the Ordovician a fundamental infaunalization process took place, leading to bivalves as we currently know them. This major ecological shift is linked to the exposure to a different redox environoments (hypoxic or anoxic) and with the Lower Devonian oxygenation event. We investigated selective signatures on bivalve and gastropod mitochondrial genomes with respect to a time calibrated mitochondrial phylogeny by means of dN/dS ratios. We were able to detect 1) a major signal of directional selection between the Ordovician and the Lower Devonian for bivalve mitochondrial Complex I, and 2) an overall higher directional selective pressure on bivalve Complex V with respect to gastropods. These and other minor dN/dS patterns and timings are discussed, showing that the Ordovician infaunalization event left heavy traces in bivalve mitochondrial genomes.

Comparative Large-Scale Mitogenomics Evidences Clade-Specific Evolutionary Trends in Mitochondrial DNAs of Bivalvia

Despite the figure of complete bivalve mitochondrial genomes keeps growing, an assessment of the general features of these genomes in a phylogenetic framework is still lacking, despite the fact that bivalve mitochondrial genomes are unusual under different aspects. In this work, we constructed a dataset of one hundred mitochondrial genomes of bivalves to perform the first systematic comparative mitogenomic analysis, developing a phylogenetic background to scaffold the evolutionary history of the class' mitochondrial genomes. Highly conserved domains were identified in all protein coding genes; however, four genes (namely, atp6, nad2, nad4L, and nad6) were found to be very divergent for many respects, notwithstanding the overall purifying selection working on those genomes. Moreover, the atp8 gene was newly annotated in 20 mitochondrial genomes, where it was previously declared as lacking or only signaled. Supernumerary mitochondrial proteins were compared, but it was possible to find homologies only among strictly related species. The rearrangement rate on the molecule is too high to be used as a phylogenetic marker, but here we demonstrate for the first time in mollusks that there is correlation between rearrangement rates and evolutionary rates. We also developed a new index (HERMES) to estimate the amount of mitochondrial evolution. Many genomic features are phylogenetically congruent and this allowed us to highlight three main phases in bivalve history: the origin, the branching of palaeoheterodonts, and the second radiation leading to the present-day biodiversity.

Predicting the minimal translation apparatus: lessons from the reductive evolution of mollicutes

Mollicutes is a class of parasitic bacteria that have evolved from a common Firmicutes ancestor mostly by massive genome reduction. With genomes under 1 Mbp in size, most Mollicutes species retain the capacity to replicate and grow autonomously. The major goal of this work was to identify the minimal set of proteins that can sustain ribosome biogenesis and translation of the genetic code in these bacteria. Using the experimentally validated genes from the model bacteria Escherichia coli and Bacillus subtilis as input, genes encoding proteins of the core translation machinery were predicted in 39 distinct Mollicutes species, 33 of which are culturable. The set of 260 input genes encodes proteins involved in ribosome biogenesis, tRNA maturation and aminoacylation, as well as proteins cofactors required for mRNA translation and RNA decay. A core set of 104 of these proteins is found in all species analyzed. Genes encoding proteins involved in post-translational modifications of ribosomal proteins and translation cofactors, post-transcriptional modifications of t+rRNA, in ribosome assembly and RNA degradation are the most frequently lost. As expected, genes coding for aminoacyl-tRNA synthetases, ribosomal proteins and initiation, elongation and termination factors are the most persistent (i.e. conserved in a majority of genomes). Enzymes introducing nucleotides modifications in the anticodon loop of tRNA, in helix 44 of 16S rRNA and in helices 69 and 80 of 23S rRNA, all essential for decoding and facilitating peptidyl transfer, are maintained in all species. Reconstruction of genome evolution in Mollicutes revealed that, beside many gene losses, occasional gains by horizontal gene transfer also occurred. This analysis not only showed that slightly different solutions for preserving a functional, albeit minimal, protein synthetizing machinery have emerged in these successive rounds of reductive evolution but also has broad implications in guiding the reconstruction of a minimal cell by synthetic biology approaches.

Investigating the role of site specific synonymous variation in disease association studies

Synonymous codon changes may not always be neutral indicating their significance in disease association studies, which is almost always overlooked. Synonymous substitutions may affect protein-folding rates leading to protein misfolding and aggregation. Genome wide analysis of 2301 mitochondrial genomes is performed to evaluate the significance of synonymous codons in disease association studies. The analysis revealed usage of rare codons at several sites in mitochondrial genes with rare codon usage higher for hydrophobic amino acids. The analysis suggests that variation data in association studies should be analyzed using site-specific codon usage values to infer the potential phenotypic impact of synonymous changes.

Genomes of the class Erysipelotrichia clarify the firmicute origin of the class Mollicutes

The tree of life is paramount for achieving an integrated understanding of microbial evolution and the relationships between physiology, genealogy and genomics. It provides the framework for interpreting environmental sequence data, whether applied to microbial ecology or to human health. However, there remain many instances where there is ambiguity in our understanding of the phylogeny of major lineages, and/or confounding nomenclature. Here we apply recent genomic sequence data to examine the evolutionary history of members of the classes Mollicutes (phylum Tenericutes) and Erysipelotrichia (phylum Firmicutes). Consistent with previous analyses, we find evidence of a specific relationship between them in molecular phylogenies and signatures of the 16S rRNA, 23S rRNA, ribosomal proteins and aminoacyl-tRNA synthetase proteins. Furthermore, by mapping functions over the phylogenetic tree we find that the erysipelotrichia lineages are involved in various stages of genomic reduction, having lost (often repeatedly) a variety of metabolic functions and the ability to form endospores. Although molecular phylogeny has driven numerous taxonomic revisions, we find it puzzling that the most recent taxonomic revision of the phyla Firmicutes and Tenericutes has further separated them into distinct phyla, rather than reflecting their common roots.

Structure, transcription, and variability of metazoan mitochondrial genome: perspectives from an unusual mitochondrial inheritance system

Despite its functional conservation, the mitochondrial genome (mtDNA) presents strikingly different features among eukaryotes, such as size, rearrangements, and amount of intergenic regions. Nonadaptive processes such as random genetic drift and mutation rate play a fundamental role in shaping mtDNA: the mitochondrial bottleneck and the number of germ line replications are critical factors, and different patterns of germ line differentiation could be responsible for the mtDNA diversity observed in eukaryotes. Among metazoan, bivalve mollusc mtDNAs show unusual features, like hypervariable gene arrangements, high mutation rates, large amount of intergenic regions, and, in some species, an unique inheritance system, the doubly uniparental inheritance (DUI). The DUI system offers the possibility to study the evolutionary dynamics of mtDNAs that, despite being in the same organism, experience different genetic drift and selective pressures. We used the DUI species Ruditapes philippinarum to study intergenic mtDNA functions, mitochondrial transcription, and polymorphism in gonads. We observed: 1) the presence of conserved functional elements and novel open reading frames (ORFs) that could explain the evolutionary persistence of intergenic regions and may be involved in DUI-specific features; 2) that mtDNA transcription is lineage-specific and independent from the nuclear background; and 3) that male-transmitted and female-transmitted mtDNAs have a similar amount of polymorphism but of different kinds, due to different population size and selection efficiency. Our results are consistent with the hypotheses that mtDNA evolution is strongly dependent on the dynamics of germ line formation, and that the establishment of a male-transmitted mtDNA lineage can increase male fitness through selection on sperm function.

Virtual genomes in flux: an interplay of neutrality and adaptability explains genome expansion and streamlining

The picture that emerges from phylogenetic gene content reconstructions is that genomes evolve in a dynamic pattern of rapid expansion and gradual streamlining. Ancestral organisms have been estimated to possess remarkably rich gene complements, although gene loss is a driving force in subsequent lineage adaptation and diversification. Here, we study genome dynamics in a model of virtual cells evolving to maintain homeostasis. We observe a pattern of an initial rapid expansion of the genome and a prolonged phase of mutational load reduction. Generally, load reduction is achieved by the deletion of redundant genes, generating a streamlining pattern. Load reduction can also occur as a result of the generation of highly neutral genomic regions. These regions can expand and contract in a neutral fashion. Our study suggests that genome expansion and streamlining are generic patterns of evolving systems. We propose that the complex genotype to phenotype mapping in virtual cells as well as in their biological counterparts drives genome size dynamics, due to an emerging interplay between adaptation, neutrality, and evolvability.

Patterns and rates of nucleotide substitution, insertion and deletion in the endosymbiont of ants Blochmannia floridanus

Genome reduction is a general process that has been studied in numerous symbiotic bacteria associated with insects. We investigated the last stages of genome degradation in Blochmannia floridanus, a mutualistic bacterial endosymbiont of the ant Camponotus floridanus. We determined the tempo (rates of insertion and deletion) and mode (size and number of insertion-deletion events) of the process in the last 200,000 years by analysing a total of 16 intergenic regions in several strains of this endosymbiont from different ant populations. We provide the first calculation of the reduction rate for noncoding DNA in this endosymbiont (2.2 x 10(-8) lost nucleotides/site/year) and compare it with the rate of loss in other species. Our results confirm, as it has been observed in other organisms like Buchnera aphidicola or Rickettsia spp., that deletions larger than one nucleotide can still appear in advanced stages of genome reduction and that a substitutional deletion bias exists. However, this bias is not due to a higher proportion of deletion over insertion events but to a few deletion events being larger than the rest. Moreover, we detected a substitutional AT bias that is probably responsible for the increase in the number of the small and moderate indel events in the last stages of genome reduction. Accordingly, we found intrapopulational polymorphisms for the detected microsatellites in contrast to the stability associated with these in free-living bacteria such as Escherichia coli.

The consequences of genetic drift for bacterial genome complexity

Genetic drift, which is particularly effective within small populations, can shape the size and complexity of genomes by affecting the fixation of deleterious mutations. In Bacteria, assessing the contribution of genetic drift to genome evolution is problematic because the usual methods, based on intraspecific polymorphisms, can be thwarted by difficulties in delineating species' boundaries. The increased availability of sequenced bacterial genomes allows application of an alternative estimator of drift, the genome-wide ratio of replacement to silent substitutions in protein-coding sequences. This ratio, which reflects the action of purifying selection across the entire genome, shows a strong inverse relationship with genome size, indicating that drift promotes genome reduction in bacteria.

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

The transfer of mitochondrial genes to the nucleus is a recurrent and consistent feature of eukaryotic genome evolution. Although many theories have been proposed to explain such transfers, little relevant data exist. The observation that clonal and self-fertilizing plants transfer more mitochondrial genes to their nuclei than do outcrossing plants contradicts predictions of major theories based on nuclear recombination and leaves a gap in our conceptual understanding how the observed pattern of gene transfer could arise. Here, with a series of deterministic and stochastic simulations, we show how epistatic selection and relative mutation rates of mitochondrial and nuclear genes influence mitochondrial-to-nuclear gene transfer. Specifically, we show that when there is a benefit to having a mitochondrial gene present in the nucleus, but absent in the mitochondria, self-fertilization dramatically increases both the rate and the probability of gene transfer. However, absent such a benefit, when mitochondrial mutation rates exceed those of the nucleus, self-fertilization decreases the rate and probability of transfer. This latter effect, however, is much weaker than the former. Our results are relevant to understanding the probabilities of fixation when loci in different genomes interact.

The cyanobacterial endosymbiont of the unicellular algae Rhopalodia gibba shows reductive genome evolution

BACKGROUND

Bacteria occur in facultative association and intracellular symbiosis with a diversity of eukaryotic hosts. Recently, we have helped to characterise an intracellular nitrogen fixing bacterium, the so-called spheroid body, located within the diatom Rhopalodia gibba. Spheroid bodies are of cyanobacterial origin and exhibit features that suggest physiological adaptation to their intracellular life style. To investigate the genome modifications that have accompanied the process of endosymbiosis, here we compare gene structure, content and organisation in spheroid body and cyanobacterial genomes.

RESULTS

Comparison of the spheroid body's genome sequence with corresponding regions of near free-living relatives indicates that multiple modifications have occurred in the endosymbiont's genome. These include localised changes that have led to elimination of some genes. This gene loss has been accompanied either by deletion of the respective DNA region or replacement with non-coding DNA that is AT rich in composition. In addition, genome modifications have led to the fusion and truncation of genes. We also report that in the spheroid body's genome there is an accumulation of deleterious mutations in genes for cell wall biosynthesis and processes controlled by transposases. Interestingly, the formation of pseudogenes in the spheroid body has occurred in the presence of intact, and presumably functional, recA and recF genes. This is in contrast to the situation in most investigated obligate intracellular bacterium-eukaryote symbioses, where at least either recA or recF has been eliminated.

CONCLUSION

Our analyses suggest highly specific targeting/loss of individual genes during the process of genome reduction and establishment of a cyanobacterial endosymbiont inside a eukaryotic cell. Our findings confirm, at the genome level, earlier speculation on the obligate intracellular status of the spheroid body in Rhopalodia gibba. This association is the first example of an obligate cyanobacterial symbiosis involving nitrogen fixation for which genomic data are available. It represents a new model system to study molecular adaptations of genome evolution that accompany a switch from free-living to intracellular existence.