Mitochondrially mediated RNA interference, a retrograde signaling system affecting nuclear gene expression.

Several functional classes of short noncoding RNAs are involved in manifold regulatory processes in eukaryotes, including, among the best characterized, miRNAs. One of the most intriguing regulatory networks in the eukaryotic cell is the mito-nuclear crosstalk: recently, miRNA-like elements of mitochondrial origin, called smithRNAs, were detected in a bivalve species, Ruditapes philippinarum. These RNA molecules originate in the organelle but were shown in vivo to regulate nuclear genes. Since miRNA genes evolve easily de novo with respect to protein-coding genes, in the present work we estimate the probability with which a newly arisen smithRNA finds a suitable target in the nuclear transcriptome. Simulations with transcriptomes of 12 bivalve species suggest that this probability is high and not species specific: one in a hundred million (1 × 10-8) if five mismatches between the smithRNA and the 3' mRNA are allowed, yet many more are allowed in animals. We propose that novel smithRNAs may easily evolve as exaptation of the pre-existing mitochondrial RNAs. In turn, the ability of evolving novel smithRNAs may have played a pivotal role in mito-nuclear interactions during animal evolution, including the intriguing possibility of acting as speciation trigger.

Mitochondrion-to-nucleus communication mediated by RNA export: a survey of potential mechanisms and players across eukaryotes.

The nucleus interacts with the other organelles to perform essential functions of the eukaryotic cell. Mitochondria have their own genome and communicate back to the nucleus in what is known as mitochondrial retrograde response. Information is transferred to the nucleus in many ways, leading to wide-ranging changes in nuclear gene expression and culminating with changes in metabolic, regulatory or stress-related pathways. RNAs are emerging molecules involved in this signalling. RNAs encode precise information and are involved in highly target-specific signalling, through a wide range of processes known as RNA interference. RNA-mediated mitochondrial retrograde response requires these molecules to exit the mitochondrion, a process that is still mostly unknown. We suggest that the proteins/complexes translocases of the inner membrane, polynucleotide phosphorylase, mitochondrial permeability transition pore, and the subunits of oxidative phosphorylation complexes may be responsible for RNA export.

Mitonuclear Coevolution, but not Nuclear Compensation, Drives Evolution of OXPHOS Complexes in Bivalves.

In Metazoa, four out of five complexes involved in oxidative phosphorylation (OXPHOS) are formed by subunits encoded by both the mitochondrial (mtDNA) and nuclear (nuDNA) genomes, leading to the expectation of mitonuclear coevolution. Previous studies have supported coadaptation of mitochondria-encoded (mtOXPHOS) and nuclear-encoded OXPHOS (nuOXPHOS) subunits, often specifically interpreted with regard to the "nuclear compensation hypothesis," a specific form of mitonuclear coevolution where nuclear genes compensate for deleterious mitochondrial mutations due to less efficient mitochondrial selection. In this study, we analyzed patterns of sequence evolution of 79 OXPHOS subunits in 31 bivalve species, a taxon showing extraordinary mtDNA variability and including species with "doubly uniparental" mtDNA inheritance. Our data showed strong and clear signals of mitonuclear coevolution. NuOXPHOS subunits had concordant topologies with mtOXPHOS subunits, contrary to previous phylogenies based on nuclear genes lacking mt interactions. Evolutionary rates between mt and nuOXPHOS subunits were also highly correlated compared with non-OXPHO-interacting nuclear genes. Nuclear subunits of chimeric OXPHOS complexes (I, III, IV, and V) also had higher dN/dS ratios than Complex II, which is formed exclusively by nuDNA-encoded subunits. However, we did not find evidence of nuclear compensation: mitochondria-encoded subunits showed similar dN/dS ratios compared with nuclear-encoded subunits, contrary to most previously studied bilaterian animals. Moreover, no site-specific signals of compensatory positive selection were detected in nuOXPHOS genes. Our analyses extend the evidence for mitonuclear coevolution to a new taxonomic group, but we propose a reconsideration of the nuclear compensation hypothesis.

Phylomitogenomics provides new perspectives on the Euphasmatodea radiation (Insecta: Phasmatodea).

Phasmatodea species diversity lies almost entirely within its suborder Euphasmatodea, which exhibits a pantropical distribution and is considered to derive from a recent and rapid evolutionary radiation. To shed light on Euphasmatodea origins and diversification, we assembled the mitogenomes of 17 species from transcriptomic sequencing data and analysed them along with 22 already available Phasmatodea mitogenomes and 33 mitogenomes representing most of the Polyneoptera lineages. Maximum Likelihood and Bayesian Inference approaches retrieved consistent topologies, both showing the widespread conflict between phylogenetic approaches and traditional systematics. We performed a divergence time analysis leveraging ten fossil specimens representative of most polyneopteran lineages: the time tree obtained supports an older radiation of the clade with respect to previous hypotheses. Euphasmatodea diversification is inferred to have started ~ 187 million years ago, suggesting that the Triassic-Jurassic mass extinction and the breakup of Pangea could have contributed to the process. We also investigated Euphasmatodea mitogenomes patterns of dN, dS and dN/dS ratio throughout our time-tree, trying to characterize the selective regime which may have shaped the clade evolution.

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.

The Evolution and Characterization of the RNA Interference Pathways in Lophotrochozoa.

In animals, three main RNA interference mechanisms have been described so far, which respectively maturate three types of small noncoding RNAs (sncRNAs): miRNAs, piRNAs, and endo-siRNAs. The diversification of these mechanisms is deeply linked with the evolution of the Argonaute gene superfamily since each type of sncRNA is typically loaded by a specific Argonaute homolog. Moreover, other protein families play pivotal roles in the maturation of sncRNAs, like the DICER ribonuclease family, whose DICER1 and DICER2 paralogs maturate respectively miRNAs and endo-siRNAs. Within Metazoa, the distribution of these families has been only studied in major groups, and there are very few data for clades like Lophotrochozoa. Thus, we here inferred the evolutionary history of the animal Argonaute and DICER families including 43 lophotrochozoan species. Phylogenetic analyses along with newly sequenced sncRNA libraries suggested that in all Trochozoa, the proteins related to the endo-siRNA pathway have been lost, a part of them in some phyla (i.e. Nemertea, Bryozoa, Entoprocta), while all of them in all the others. On the contrary, early diverging phyla, Platyhelminthes and Syndermata, showed a complete endo-siRNA pathway. On the other hand, miRNAs were revealed the most conserved and ubiquitous mechanism of the metazoan RNA interference machinery, confirming their pivotal role in animal cell regulation.

The complete mitochondrial genome of Solemya velum (Mollusca: Bivalvia) and its relationships with conchifera.

BACKGROUND: Bivalve mitochondrial genomes exhibit a wide array of uncommon features, like extensive gene rearrangements, large sizes, and unusual ways of inheritance. Species pertaining to the order Solemyida (subclass Opponobranchia) show many peculiar evolutionary adaptations, f.i. extensive symbiosis with chemoautotrophic bacteria. Despite Opponobranchia are central in bivalve phylogeny, being considered the sister group of all Autobranchia, a complete mitochondrial genome has not been sequenced yet. RESULTS: In this paper, we characterized the complete mitochondrial genome of the Atlantic awning clam Solemya velum: A-T content, gene arrangement and other features are more similar to putative ancestral mollusks than to other bivalves. Two supranumerary open reading frames are present in a large, otherwise unassigned, region, while the origin of replication could be located in a region upstream to the cox3 gene. CONCLUSIONS: We show that S. velum mitogenome retains most of the ancestral conchiferan features, which is unusual among bivalve mollusks, and we discuss main peculiarities of this first example of an organellar genome coming from the subclass Opponobranchia. Mitochondrial genomes of Solemya (for bivalves) and Haliotis (for gastropods) seem to retain the original condition of mollusks, as most probably exemplified by Katharina.

Biodiversity of extant snails (Gastropoda, Mollusca) in the Pliocene Mountain Spur Natural Reserve (Northern Apennine, Italy).

Background: To date, there is a substantial lack of information about gastropods from the Tuscan-Emilian Apennine in the north of Italy, notwithstanding the availability of detailed and comprehensive literature on this molluscan class. We present a gastropod fauna from the Natural Reserve of the Pliocene Mountain Spur: to our knowledge, this is the first investigation of the extant gastropod fauna in the Tuscan-Emilian Apennine and one of the few in the Apennine's mountain chain as a whole. New information: We describe a gastropod fauna comprised by 25 species, belonging to 18 genera and 10 families: the general figure which is emerging is an assemblage of European and Mediterranean-European species, with a single Asian contribution. Several shells were collected as embedded in sandy-clayey soils and showed fossilisation traces: therefore, we regard these shells as subfossil samples. Namely, subfossil shells are from the species Pomatiaselegans, Granariafrumentum, Retinellaolivetorum, Xerolentaobviaobvia, Oxychiluscf.draparnaudi, Monachacartusiana and Monachacantiana. The present checklist is the first report in the Tuscan-Emilian Apennine and Emilia-Romagna areas for six taxa: Morlinaglabraglabra, Oxychiliusalliarius, Xerosectacespitum, Fruticicolafruticum, Xerogyraspadae and Xerolentaobviaobvia.

Metal and metalloid concentrations in wild mammals from SW Europe: European hedgehog (Erinaceus europaeus) and badger (Meles meles).

In recent years, there have been increasing ecological and global concerns associated to Potentially Toxic Elements (PTEs). Thus, the relevance of wild mammals as biomonitors has been globally recognised. In the present study, Cd, Pb, Hg, Zn and As concentrations were quantified in European hedgehog and badger inhabiting SW Europe, and cumulative trends in relation to age and sex were evaluated. Liver and kidney samples were collected, mineralised and PTE content was determined by ICP-MS. Zn was the most abundant element quantified in both organs (239 and 89.8 mg kg-1 for hedgehogs and 179 and 164 mg kg-1 dw for badgers). In hedgehogs, very high Hg concentration were quantified (4.35 and 15.5 mg kg-1 dw in liver and kidney), and Cd was the most abundant for badgers (4.70 and 7.61 mg kg-1 dw in liver and kidney). Positive correlations were observed for the concentrations of PTE in the organs of both species. Age-dependence increased only Cd concentration, with levels in adult kidneys being significantly higher. In this study, European hedgehog and badger were used as biomonitors for the determination of PTEs to provide current reference values in relatively non-polluted areas of SW Europe, and to enhance the use of these species for future ecotoxicological studies.

Mito-nuclear coevolution and phylogenetic artifacts: the case of bivalve mollusks.

Mito-nuclear phylogenetic discordance in Bivalvia is well known. In particular, the monophyly of Amarsipobranchia (Heterodonta + Pteriomorphia), retrieved from mitochondrial markers, contrasts with the monophyly of Heteroconchia (Heterodonta + Palaeoheterodonta), retrieved from nuclear markers. However, since oxidative phosphorylation nuclear markers support the Amarsipobranchia hypothesis instead of the Heteroconchia one, interacting subunits of the mitochondrial complexes ought to share the same phylogenetic signal notwithstanding the genomic source, which is different from the signal obtained from other nuclear markers. This may be a clue of coevolution between nuclear and mitochondrial genes. In this work we inferred the phylogenetic signal from mitochondrial and nuclear oxidative phosphorylation markers exploiting different phylogenetic approaches and added two more datasets for comparison: genes of the glycolytic pathway and genes related to the biogenesis of regulative small noncoding RNAs. All trees inferred from mitochondrial and nuclear subunits of the mitochondrial complexes support the monophyly of Amarsipobranchia, regardless of the phylogenetic pipeline. However, not every single marker agrees with this topology: this is clearly visible in nuclear subunits that do not directly interact with the mitochondrial counterparts. Overall, our data support the hypothesis of a coevolution between nuclear and mitochondrial genes for the oxidative phosphorylation. Moreover, we suggest a relationship between mitochondrial topology and different nucleotide composition between clades, which could be associated to the highly variable gene arrangement in Bivalvia.

The mitochondrial genome of Bacillus stick insects (Phasmatodea) and the phylogeny of orthopteroid insects.

The Order Phasmatodea (stick and leaf insects) includes many well-known species of cryptic phytophagous insects. In this work, we sequenced the almost complete mitochondrial genomes of two stick insect species of the genus Bacillus. Phasmatodea pertain to the Polyneoptera, and represent one of the major clades of heterometabolous insects. Orthopteroid insect lineages arose through rapid evolutionary radiation events, which likely blurred the phylogenetic reconstructions obtained so far; we therefore performed a phylogenetic analysis to resolve and date all major splits of orthopteroid phylogeny, including the relationships between Phasmatodea and other polyneopterans. We explored several molecular models, with special reference to data partitioning, to correctly detect any phylogenetic signal lying in rough data. Phylogenetic Informativeness analysis showed that the maximum resolving power on the orthopteroid mtDNA dataset is expected for the Upper Cretaceous, about 80millionyears ago (Mya), but at least 70% of the maximum informativeness is also expected for the 150-200 Mya timespan, which makes mtDNA a suitable marker to study orthopteroid splits. A complete chronological calibration has also been computed following a Penalized Likelihood method. In summary, our analysis confirmed the monophyly of Phasmatodea, Dictyoptera and Orthoptera, and retrieved Mantophasmatodea as sister group of Phasmatodea. The origin of orthopteroid insects was also estimated to be in the Middle Triassic, while the order Phasmatodea seems to appear in the Upper Jurassic. The obtained results evidenced that mtDNA is a suitable marker to unravel the ancient splits leading to the orthopteroid orders, given a proper methodological approach.

Molecular Phylogenetics and Mitochondrial Evolution.

The myth of a "typical" mitochondrial genome (mtDNA) is a rock-hard belief in the field of genetics, at least for the animal kingdom [...].

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.

Mitochondrial Genomic Landscape: A Portrait of the Mitochondrial Genome 40 Years after the First Complete Sequence.

Notwithstanding the initial claims of general conservation, mitochondrial genomes are a largely heterogeneous set of organellar chromosomes which displays a bewildering diversity in terms of structure, architecture, gene content, and functionality. The mitochondrial genome is typically described as a single chromosome, yet many examples of multipartite genomes have been found (for example, among sponges and diplonemeans); the mitochondrial genome is typically depicted as circular, yet many linear genomes are known (for example, among jellyfish, alveolates, and apicomplexans); the chromosome is normally said to be "small", yet there is a huge variation between the smallest and the largest known genomes (found, for example, in ctenophores and vascular plants, respectively); even the gene content is highly unconserved, ranging from the 13 oxidative phosphorylation-related enzymatic subunits encoded by animal mitochondria to the wider set of mitochondrial genes found in jakobids. In the present paper, we compile and describe a large database of 27,873 mitochondrial genomes currently available in GenBank, encompassing the whole eukaryotic domain. We discuss the major features of mitochondrial molecular diversity, with special reference to nucleotide composition and compositional biases; moreover, the database is made publicly available for future analyses on the MoZoo Lab GitHub page.

One in a Million: Genetic Diversity and Conservation of the Reference Crassostrea angulata Population in Europe from the Sado Estuary (Portugal).

The production of cupped oysters is an important component of European aquaculture. Most of the production relies on the cultivation of the Pacific oyster Crassostrea gigas, although the Portuguese oyster Crassostrea angulata represents a valuable product with both cultural and economic relevance, especially in Portugal. The authors of the present study investigated the genetic diversity of Portuguese oyster populations of the Sado estuary, both from natural oyster beds and aquaculture facilities, through cox1 gene fragment sequencing. Then, a comparison with a wide dataset of cupped oyster sequences obtained from GenBank (up to now the widest available dataset in literature for the Portuguese oyster) was performed. Genetic data obtained from this work confirmed that the Pacific oyster does not occur in the natural oyster beds of the Sado estuary but showed that the species occasionally occurs in the oyster hatcheries. Moreover, the results showed that despite the founder effect and the bottleneck events that the Sado populations have experienced, they still exhibit high haplotype diversity. Risks are arising for the conservation of the Portuguese oyster reference populations of the Sado estuary due to the occurrence of the Pacific oyster in the local hatcheries. Therefore, researchers, local authorities, and oyster producers should work together to avoid the loss of this valuable resource.

Mitochondrial phylogeny and taxonomic revision of Italian and Slovenian fluvio-lacustrine barbels, Barbus sp. (Cypriniformes, Cyprinidae).

BACKGROUND: Barbels are ray finned cyprinid fishes of the Old-World with partially unresolved, intricate taxonomy. Within the Barbus sensu lato paraphyletic assemblage, Barbus sensu stricto is a monophyletic tetraploid lineage of Europe, northern Africa and Middle East, including two monophyletic sibling genera: Barbus and Luciobarbus. Italy, Slovenia and northern Croatia are natively inhabited by several entities of the genus Barbus, whose relationships and taxonomic ranks are still unclear. Aim of the present work is to focus on phylogeography of Italian and Slovenian barbels, with an appraisal of their current taxonomy. RESULTS: One hundred fifty specimens were collected in 78 sampling sites from 33 main watersheds, widely distributed along Italian and Slovenian ichthyogeographic districts. We amplified two mitochondrial markers, cytochrome b (cytb) and control region (D-loop), to infer a robust phylogeny for our sample and investigate on species delimitation. Our results strongly indicate all Italian and Adriatic Slovenian fluvio-lacustrine barbels to be comprised into at least three distinct species. We provide a proposal of taxonomic revision and a list of synonymies for two of them and a new description under the International Code of Zoological Nomenclature rules for the third one. CONCLUSIONS: If nuclear data will confirm our findings, at least three specific entities should be acknowledged across our sampling area. Namely, the three species are (i) Barbus plebejus, in the Padano-Venetian district; (ii) Barbus tyberinus, in the Tuscany-Latium district; (iii) Barbus oscensis Rossi & Plazzi sp. nov., in the Tyrrhenian and southernmost-Adriatic parts of Apulia-Campania district. Finally, we briefly discuss the implications of such a taxonomic scenario on conservation policies.

Phylogenetic representativeness: a new method for evaluating taxon sampling in evolutionary studies.

BACKGROUND: Taxon sampling is a major concern in phylogenetic studies. Incomplete, biased, or improper taxon sampling can lead to misleading results in reconstructing evolutionary relationships. Several theoretical methods are available to optimize taxon choice in phylogenetic analyses. However, most involve some knowledge about the genetic relationships of the group of interest (i.e., the ingroup), or even a well-established phylogeny itself; these data are not always available in general phylogenetic applications. RESULTS: We propose a new method to assess taxon sampling developing Clarke and Warwick statistics. This method aims to measure the "phylogenetic representativeness" of a given sample or set of samples and it is based entirely on the pre-existing available taxonomy of the ingroup, which is commonly known to investigators. Moreover, our method also accounts for instability and discordance in taxonomies. A Python-based script suite, called PhyRe, has been developed to implement all analyses we describe in this paper. CONCLUSIONS: We show that this method is sensitive and allows direct discrimination between representative and unrepresentative samples. It is also informative about the addition of taxa to improve taxonomic coverage of the ingroup. Provided that the investigators' expertise is mandatory in this field, phylogenetic representativeness makes up an objective touchstone in planning phylogenetic studies.

Towards a molecular phylogeny of Mollusks: bivalves' early evolution as revealed by mitochondrial genes.

Despite huge fossil, morphological and molecular data, bivalves' early evolutionary history is still a matter of debate: recently, established phylogeny has been mostly challenged by DNA studies, and little agreement has been reached in literature, because of a substantial lack of widely-accepted methodological approaches to retrieve and analyze bivalves' molecular data. Here we present a molecular phylogeny of the class based on four mitochondrial genes (12s, 16s, cox1, cytb) and a methodological pipeline that proved to be useful to obtain robust results. Actually, best-performing taxon sampling and alignment strategies were tested, and several data partitioning and molecular evolution models were analyzed, thus demonstrating the utility of Bayesian inference and the importance of molding and implementing non-trivial evolutionary models. Therefore, our analysis allowed to target many taxonomic questions of Bivalvia, and to obtain a complete time calibration of the tree depicting bivalves' earlier natural history main events, which mostly dated in the late Cambrian.

Clues of in vivo nuclear gene regulation by mitochondrial short non-coding RNAs.

Gene expression involves multiple processes, from transcription to translation to the mature, functional peptide, and it is regulated at multiple levels. Small RNA molecules are known to bind RNA messengers affecting their fate in the cytoplasm (a process generically termed 'RNA interference'). Such small regulatory RNAs are well-known to be originated from the nuclear genome, while the role of mitochondrial genome in RNA interference was largely overlooked. However, evidence is growing that mitochondrial DNA does provide the cell a source of interfering RNAs. Small mitochondrial highly transcribed RNAs (smithRNAs) have been proposed to be transcribed from the mitochondrion and predicted to regulate nuclear genes. Here, for the first time, we show in vivo clues of the activity of two smithRNAs in the Manila clam, Ruditapes philippinarum. Moreover, we show that smithRNAs are present and can be annotated in representatives of the three main bilaterian lineages; in some cases, they were already described and assigned to a small RNA category (e.g., piRNAs) given their biogenesis, while in other cases their biogenesis remains unclear. If mitochondria may affect nuclear gene expression through RNA interference, this opens a plethora of new possibilities for them to interact with the nucleus and makes metazoan mitochondrial DNA a much more complex genome than previously thought.