Gene annotation errors are common in the mammalian mitochondrial genomes database

Museomics help resolving the phylogeny of snowfinches (Aves, Passeridae, Montifringilla and allies)

Historical specimens from museum collections provide a valuable source of material also from remote areas or regions of conflict that are not easily accessible to scientists today. With this study, we are providing a taxon-complete phylogeny of snowfinches using historical DNA from whole skins of an endemic species from Afghanistan, the Afghan snowfinch, Pyrgilauda theresae. To resolve the strong conflict between previous phylogenetic hypotheses, we generated novel mitogenome sequences for selected taxa and genome-wide SNP data using ddRAD sequencing for all extant snowfinch species endemic to the Qinghai-Tibet Plateau (QTP) and for an extended intraspecific sampling of the sole Central and Western Palearctic snowfinch species (Montifringilla nivalis). Our phylogenetic reconstructions unanimously refuted the previously suggested paraphyly of genus Pyrgilauda. Misplacement of one species-level taxon (Onychostruthus tazcanowskii) in previous snowfinch phylogenies was undoubtedly inferred from chimeric mitogenomes that included heterospecific sequence information. Furthermore, comparison of novel and previously generated sequence data showed that the presumed sister-group relationship between M. nivalis and the QTP endemic M. henrici was suggested based on flawed taxonomy. Our phylogenetic reconstructions based on genome-wide SNP data and on mitogenomes were largely congruent and supported reciprocal monophyly of genera Montifringilla and Pyrgilauda with monotypic Onychostruthus being sister to the latter. The Afghan endemic P. theresae likely originated from a rather ancient Pliocene out-of-Tibet dispersal probably from a common ancestor with P. ruficollis. Our extended trans-Palearctic sampling for the white-winged snowfinch, M. nivalis, confirmed strong lineage divergence between an Asian and a European clade dated to 1.5 - 2.7 million years ago (mya). Genome-wide SNP data suggested subtle divergence among European samples from the Alps and from the Cantabrian mountains.

Complete mitochondrial genomes of three vulnerable cave bat species and their phylogenetic relationships within the order Chiroptera

The IUCN Red List of Threatened Species contains 175 Brazilian bat species that are threatened by extinction in some degree. From this perspective, it is essential to expand the knowledge about the genetic diversity of vulnerable bats. Genomic sequencing can be useful to generate robust and informative genetic references, increasing resolution when analyzing relationships among populations, species, or higher taxonomic levels. In this study, we sequenced and characterized in detail the first complete mitochondrial genomes of Furipterus horrens, Lonchorhina aurita, and Natalus macrourus, and investigated their phylogenetic position based on amino acid sequences of protein-coding genes (PCGs). The mitogenomes of these species are 16,516, 16,697, and 16,668 bp in length, respectively, and each comprises 13 PCGs, 22 tRNA genes, two rRNA genes, and a putative control region (CR). In the three species, genes were arranged similarly to all other previously described bat mitogenomes, and nucleotide composition was also consistent with the reported range. The length and arrangement of rrnS and rrnL were also consistent with those of other bat species, showing a positive AT-skew and a negative GC-skew. Except for trnS1, for which we did not observe the DHU arm, all other tRNAs showed the cloverleaf secondary structure in the three species. In addition, the mitogenomes showed minor differences in start and stop codons, and in all PCGs, codons ending in adenine were more common compared to those ending in guanine. We found that PCGs of the three species use multiple codons to encode each amino acid, following the previously documented pattern. Furthermore, all PCGs are under purifying selection, with atp8 experiencing the most relaxed purifying selection. Considering the phylogenetic reconstruction, F. horrens was recovered as sister to Noctilio leporinus, L. aurita and Tonatia bidens shared a node within Phyllostomidae, and N. macrourus appeared as sister to Molossidae and Vespertilionidae.

Variability and evolution of gene order rearrangement in mitochondrial genomes of arthropods (except Hexapoda)

In the species-rich Phylum Arthropoda, the mitochondrial genome is relatively well conserved both in terms of number and order of genes. However, specific clades have a 'typical' gene order that differs from the putative arthropod ancestral arrangement. The aim of this work was to compare the rate of mitochondrial gene rearrangements at inter- and intra-taxonomic levels in the Arthropoda and to postulate the most parsimonious ancestral orders representing the four major arthropod lineages. For this purpose, we performed a comparative genomic analysis of arthropod mitochondrial genomes available in the NCBI database. Using a combination of bioinformatics methods that examined mitochondrial gene rearrangements in 464 species of arthropods from three subphyla (Chelicerata, Myriapoda, and Crustacea [except Hexapoda, previously analyzed]), we observed differences in the rate of rearrangement within major lineages. A higher rate of mitochondrial genome rearrangement was observed in Crustacea and Chelicerata compared to Myriapoda. Likewise, early branching clades exhibit less variability in mitochondrial genome order than late branching clades, within each subphylum. We identified 'hot regions' in the mitochondrial genome of each studied subphylum, and postulated the most likely ancestral gene order in each subphylum and taxonomic order. Our work provides new evidence on the evolutionary dynamics of mitochondrial genome gene order in arthropods and new mitochondrial genome architectures in different taxonomic divisions within each major lineage of arthropods.

Complete mitochondrial genome of Phoneutria depilata (Araneae, Ctenidae): New Insights into the Phylogeny and Evolution of Spiders

Spiders (Araneae) are the most abundant terrestrial predators and megadiverse on earth. In recent years, the mitochondrial genome of a great diversity of species has been sequenced, mainly for ecological and commercial purposes. These studies have uncovered the existence of a variety of mitochondrial genome rearrangements. However, there is poor genetic information in several taxonomic families of spiders. We have sequenced the complete genome of Phoneutria depilata (Ctenidae) and, based on this, extract the mitogenomes of other ctenid species from published transcriptomes to perform a comparative study among spider species to determine the relationship between the level of mitochondrial rearrangements and its possible relationship with molecular variability in spiders. Complete mitochondrial genomes of eighteen spiders (including eight Ctenidae species) were obtained by two different methodologies (sequencing and transcriptome extraction). Fifty-eight spider mitochondrial genomes were downloaded from the NCBI database for gene order analysis. After verifying the annotation of each mitochondrial gene, a phylogenetic and a gene order analysis from 76 spider mitochondrial genomes were carried out. Our results show a high rate of annotation error in the published spider mitochondrial genomes, which could lead to errors in phylogenetic inference. Moreover, to provide new mitochondrial genomes in spiders by two different methodologies to obtain them, our analysis identifies six different mitochondrial architectures among all spiders. Translocation or tandem duplication random loss (TDRL) events in tRNA genes were identified to explain the evolution of the spider mitochondrial genome. In addition, our findings provide new insights into spider mitochondrial evolution.

Evolution of Helicobacter spp: variability of virulence factors and their relationship to pathogenicity

Background

Virulence factors (VF) are bacteria-associated molecules that assist to colonize the host at the cellular level. Bacterial virulence is highly dynamic and specific pathogens have a broad array of VFs. The genus Helicobacter is gram-negative, microaerobic, flagellated, and mucus-inhabiting bacteria associated with gastrointestinal inflammation. To investigate about their pathogenicity, several Helicobacter species have been characterized and sequenced. Since the variability and possible origin of VF in the genus are not clear, our goal was to perform a comparative analysis of Helicobacter species in order to investigate VF variability and their evolutionary origin.

Methods

The complete genomes of 22 Helicobacter species available in NCBI were analyzed, using computational tools. We identifyed gain and loss events in VF genes, which were categorized in seven functional groups to determine their most parsimonious evolutionary origin. After verifying the annotation of all VF genes, a phylogeny from conserved VF organized by Helicobacter species according to gastric Helicobacter species (GHS) or enterohepatic (EHS) classification was obtained.

Results

Gain and loss analysis of VF orthologous in Helicobacter ssp revealed the most possible evolutionary origin for each gene set. Microevolutionary events in urease and flagella genes were detected during the evolution of the genus. Our results pointed that acquisition of ureases and adherence genes and deletion of cytotoxins in some lineages, as well as variation in VF genes copy number, would be related to host adaptation during evolution of the Helicobacter genus. Our findings provided new insights about the genetic differences between GHS and EHS and their relationship with pathogenicity.

The complete mitochondrial genome of giant cricket, Tarbinskiellus portentosus (Orthoptera: Gryllidae) and its curation

Tarbinskiellus portentosus, commonly known as giant cricket one of the important edible cricket species. However, the genetic information of these species is still limited. Therefore, we have assembled and annotated the first mitochondrial genome of T. portentosus. The mitogenome is 15710 bp long and has GC content of 27.19%. The nucleotide composition is similar with other insect mitogenomes (A 40.6%; T 32.2%; C 17.3%; G 9.9%). The gene organization in the mitogenome of T. portentosus is identical to the mitogenome of other cricket species. The complete mitogenome of T. portentosus consisted 37 genes including 13 protein coding genes, 22 tRNA genes, and two rRNA genes. The newly assembled mitogenome will help molecular biology research on edible crickets. Since mitogenome genes are traditionally used for DNA barcoding and phylogenetic analysis, comparative analysis of

T. portentosus mitogenome with other related cricket species will also aid researchers in developing universal primers for species identification toward food security. Apart from the main goal of providing full mitogenome of T. portentosus, paper also provides conceptual workflow based on de novo assembly and its correction for final mitogenome construction.

Mitogenome of the extinct Desert 'rat-kangaroo' times the adaptation to aridity in macropodoids

The evolution of Australia's distinctive marsupial fauna has long been linked to the onset of continent-wide aridity. However, how this profound climate change event affected the diversification of extant lineages is still hotly debated. Here, we assemble a DNA sequence dataset of Macropodoidea-the clade comprising kangaroos and their relatives-that incorporates a complete mitogenome for the Desert 'rat-kangaroo', Caloprymnus campestris. This enigmatic species went extinct nearly 90 years ago and is known from a handful of museum specimens. Caloprymnus is significant because it was the only macropodoid restricted to extreme desert environments, and therefore calibrates the group's specialisation for increasingly arid conditions. Our robustly supported phylogenies nest Caloprymnus amongst the bettongs Aepyprymnus and Bettongia. Dated ancestral range estimations further reveal that the Caloprymnus-Bettongia lineage originated in nascent xeric settings during the middle to late Miocene, ~ 12 million years ago (Ma), but subsequently radiated into fragmenting mesic habitats after the Pliocene to mid-Pleistocene. This timeframe parallels the ancestral divergences of kangaroos in woodlands and forests, but predates their adaptive dispersal into proliferating dry shrublands and grasslands from the late Miocene to mid-Pleistocene, after ~ 7 Ma. We thus demonstrate that protracted changes in both climate and vegetation likely staged the emergence of modern arid zone macropodoids.

Comparative genomic analysis of vertebrate mitochondrial reveals a differential of rearrangements rate between taxonomic class

Vertebrate mitochondrial genomes have been extensively studied for genetic and evolutionary purposes, these are normally believed to be extremely conserved, however, different cases of gene rearrangements have been reported. To verify the level of rearrangement and the mitogenome evolution, we performed a comparative genomic analysis of the 2831 vertebrate mitochondrial genomes representing 12 classes available in the NCBI database. Using a combination of bioinformatics methods, we determined there is a high number of errors in the annotation of mitochondrial genes, especially in tRNAs. We determined there is a large variation in the proportion of rearrangements per gene and per taxonomic class, with higher values observed in Actinopteri, Amphibia and Reptilia. We highlight that these are results for currently available vertebrate sequences, so an increase in sequence representativeness in some groups may alter the rearrangement rates, so in a few years it would be interesting to see if these rates are maintained or altered with the new mitogenome sequences. In addition, within each vertebrate class, different patterns in rearrangement proportion with distinct hotspots in the mitochondrial genome were found. We also determined that there are eleven convergence events in gene rearrangement, nine of which are new reports to the scientific community.

Molecular relationships of the red-bellied dasyure (Phascolosorex doriae) – a rare marsupial from western New Guinea

The mitochondrial genome of the rare endemic New Guinean dasyurid Phascolosorex doriae (Thomas 1886) has been used to clarify relationships within ‘phascolosoricinae’. The mitochondrial genome has the typical gene arrangement seen in other marsupials. Molecular analyses using complete mitogenomes of other dasyurids resolve the red-bellied dasyure as sister to the narrow-striped dasyure Phascolosorex dorsalis and show that these two species diverged in the early Pliocene. The invasion of emergent New Guinean rainforest habitats (in the late Miocene) by the common ancestor of Ph. doriae, Ph. dorsalis and Neophascogale lorentzii represents one of three separate such invasions by dasyurid lineages.

Helicobacter pylori virulence factors: relationship between genetic variability and phylogeographic origin

Background

Helicobacter pylori is a pathogenic bacteria that colonize the gastrointestinal tract from human stomachs and causes diseases including gastritis, peptic ulcers, gastric lymphoma (MALT), and gastric cancer, with a higher prevalence in developing countries. Its high genetic diversity among strains is caused by a high mutation rate, observing virulence factors (VFs) variations in different geographic lineages. This study aimed to postulate the genetic variability associated with virulence factors present in the Helicobacter pylori strains, to identify the relationship of these genes with their phylogeographic origin.

Methods

The complete genomes of 135 strains available in NCBI, from different population origins, were analyzed using bioinformatics tools, identifying a high rate; as well as reorganization events in 87 virulence factor genes, divided into seven functional groups, to determine changes in position, number of copies, nucleotide identity and size, contrasting them with their geographical lineage and pathogenic phenotype.

Results

Bioinformatics analyses show a high rate of gene annotation errors in VF. Analysis of genetic variability of VFs shown that there is not a direct relationship between the reorganization and geographic lineage. However, regarding the pathogenic phenotype demonstrated in the analysis of many copies, size, and similarity when dividing the strains that possess and not the cag pathogenicity island (cagPAI), having a higher risk of developing gastritis and peptic ulcer was evidenced. Our data has shown that the analysis of the overall genetic variability of all VFs present in each strain of H. pylori is key information in understanding its pathogenic behavior.

Aln2tbl: building a mitochondrial features table from a assembly alignment in fasta format

The sequencing, annotation and analysis of complete mitochondrial genomes is an important research tool in phylogeny and evolution. Starting with the primary sequence, genes/features are generally annotated automatically to obtain preliminary annotations in the form of a feature table. Further manual curation in a graphic alignment editor is nevertheless necessary to revise annotations. As such, the automatically generated feature table is invalidated and has to be modified manually before submission to data banks. We developed aln2tbl.py, a python script that recreates a feature table from a manually refined alignment of genes mapped on the mitochondrial genome in fasta format. The feature table is populated with notes and annotations specific to mitochondrial genomes. The table can be used to create a sqn file to be submitted directly to data banks. In summary, our scripts fills one gap in the available toolbox and, combined with other software, allows the automation of the entire process, from primary sequence to annotated genome submission, even if a manual curation step is conducted in a visual sequence editor.

airpg: automatically accessing the inverted repeats of archived plastid genomes

BACKGROUND

In most flowering plants, the plastid genome exhibits a quadripartite genome structure, comprising a large and a small single copy as well as two inverted repeat regions. Thousands of plastid genomes have been sequenced and submitted to public sequence repositories in recent years. The quality of sequence annotations in many of these submissions is known to be problematic, especially regarding annotations that specify the length and location of the inverted repeats: such annotations are either missing or portray the length or location of the repeats incorrectly. However, many biological investigations employ publicly available plastid genomes at face value and implicitly assume the correctness of their sequence annotations.

RESULTS

We introduce airpg, a Python package that automatically assesses the frequency of incomplete or incorrect annotations of the inverted repeats among publicly available plastid genomes. Specifically, the tool automatically retrieves plastid genomes from NCBI Nucleotide under variable search parameters, surveys them for length and location specifications of inverted repeats, and confirms any inverted repeat annotations through self-comparisons of the genome sequences. The package also includes functionality for automatic identification and removal of duplicate genome records and accounts for taxa that genuinely lack inverted repeats. A survey of the presence of inverted repeat annotations among all plastid genomes of flowering plants submitted to NCBI Nucleotide until the end of 2020 using airpg, followed by a statistical analysis of potential associations with record metadata, highlights that release year and publication status of the genome records have a significant effect on the frequency of complete and equal-length inverted repeat annotations.

CONCLUSION

The number of plastid genomes on NCBI Nucleotide has increased dramatically in recent years, and many more genomes will likely be submitted over the next decade. airpg enables researchers to automatically access and evaluate the inverted repeats of these plastid genomes as well as their sequence annotations and, thus, contributes to increasing the reliability of publicly available plastid genomes. The software is freely available via the Python package index at http://pypi.python.org/pypi/airpg .

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.

Molluscan mitochondrial genomes break the rules

The first animal mitochondrial genomes to be sequenced were of several vertebrates and model organisms, and the consistency of genomic features found has led to a 'textbook description'. However, a more broad phylogenetic sampling of complete animal mitochondrial genomes has found many cases where these features do not exist, and the phylum Mollusca is especially replete with these exceptions. The characterization of full mollusc mitogenomes required considerable effort involving challenging molecular biology, but has created an enormous catalogue of surprising deviations from that textbook description, including wide variation in size, radical genome rearrangements, gene duplications and losses, the introduction of novel genes, and a complex system of inheritance dubbed 'doubly uniparental inheritance'. Here, we review the extraordinary variation in architecture, molecular functioning and intergenerational transmission of molluscan mitochondrial genomes. Such features represent a great potential for the discovery of biological history, processes and functions that are novel for animal mitochondrial genomes. This provides a model system for studying the evolution and the manifold roles that mitochondria play in organismal physiology, and many ways that the study of mitochondrial genomes are useful for phylogeny and population biology. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Nanopore Sequencing Resolves Elusive Long Tandem-Repeat Regions in Mitochondrial Genomes

Long non-coding, tandem-repetitive regions in mitochondrial (mt) genomes of many metazoans have been notoriously difficult to characterise accurately using conventional sequencing methods. Here, we show how the use of a third-generation (long-read) sequencing and informatic approach can overcome this problem. We employed Oxford Nanopore technology to sequence genomic DNAs from a pool of adult worms of the carcinogenic parasite, Schistosoma haematobium, and used an informatic workflow to define the complete mt non-coding region(s). Using long-read data of high coverage, we defined six dominant mt genomes of 33.4 kb to 22.6 kb. Although no variation was detected in the order or lengths of the protein-coding genes, there was marked length (18.5 kb to 7.6 kb) and structural variation in the non-coding region, raising questions about the evolution and function of what might be a control region that regulates mt transcription and/or replication. The discovery here of the largest tandem-repetitive, non-coding region (18.5 kb) in a metazoan organism also raises a question about the completeness of some of the mt genomes of animals reported to date, and stimulates further explorations using a Nanopore-informatic workflow.

The evolutionary history of manatees told by their mitogenomes

The manatee family encompasses three extant congeneric species: Trichechus senegalensis (African manatee), T. inunguis (Amazonian manatee), and T. manatus (West Indian manatee). The fossil record for manatees is scant, and few phylogenetic studies have focused on their evolutionary history. We use full mitogenomes of all extant manatee species to infer the divergence dates and biogeographical histories of these species and the effect of natural selection on their mitogenomes. The complete mitochondrial genomes of T. inunguis (16,851 bp), T. senegalensis (16,882 bp), and T. manatus (16,882 bp), comprise 13 protein-coding genes, 2 ribosomal RNA genes (rRNA - 12S and 16S), and 22 transfer RNA genes (tRNA), and (D-loop/CR). Our analyses show that the first split within Trichechus occurred during the Late Miocene (posterior mean 6.56 Ma and 95% HPD 3.81–10.66 Ma), followed by a diversification event in the Plio-Pleistocene (posterior mean 1.34 Ma, 95% HPD 0.1–4.23) in the clade composed by T. inunguis and T. manatus; T. senegalensis is the sister group of this clade with higher support values (pp > 0.90). The branch-site test identified positive selection on T. inunguis in the 181st position of the ND4 amino acid gene (LRT = 6.06, p = 0.0069, BEB posterior probability = 0.96). The ND4 gene encodes one subunit of the NADH dehydrogenase complex, part of the oxidative phosphorylation machinery. In conclusion, our results provide novel insight into the evolutionary history of the Trichechidae during the Late Miocene, which was influenced by geological events, such as Amazon Basin formation.

Quantification and evolution of mitochondrial genome rearrangement in Amphibians

BACKGROUND

Rearrangement is an important topic in the research of amphibian mitochondrial genomes ("mitogenomes" hereafter), whose causes and mechanisms remain enigmatic. Globally examining mitogenome rearrangements and uncovering their characteristics can contribute to a better understanding of mitogenome evolution.

RESULTS

Here we systematically investigated mitogenome arrangements of 232 amphibians including four newly sequenced Dicroglossidae mitogenomes. The results showed that our new sequenced mitogenomes all possessed a trnM tandem duplication, which was not exclusive to Dicroglossidae. By merging the same arrangements, the mitogenomes of ~ 80% species belonged to the four major patterns, the major two of which were typical vertebrate arrangement and typical neobatrachian arrangement. Using qMGR for calculating rearrangement frequency (RF) (%), we found that the control region (CR) (RF = 45.04) and trnL2 (RF = 38.79) were the two most frequently rearranged components. Forty-seven point eight percentage of amphibians possessed rearranged mitogenomes including all neobatrachians and their distribution was significantly clustered in the phylogenetic trees (p < 0.001). In addition, we argued that the typical neobatrachian arrangement may have appeared in the Late Jurassic according to possible occurrence time estimation.

CONCLUSION

It was the first global census of amphibian mitogenome arrangements from the perspective of quantity statistics, which helped us to systematically understand the type, distribution, frequency and phylogenetic characteristics of these rearrangements.

Re-Evaluating the Internal Phylogenetic Relationships of Collembola by Means of Mitogenome Data

Collembola are an ancient and early diverging lineage of basal hexapods that occur in virtually all terrestrial habitats on Earth. Phylogenetic relationships between the different orders of Collembola are fiercely debated. Despite a range of studies and the application of both morphological and genetic approaches (singly or in combination) to assess the evolutionary relationships of major lineages in the group, no consensus has been reached. Several mitogenome sequences have been published for key taxa of the class (and their number is increasing rapidly). Here, we describe two new Antarctic Collembola mitogenomes and compare all complete or semi-complete springtail mitogenome sequences available on GenBank in terms of both gene order and DNA sequence analyses in a genome evolution and molecular phylogenetic framework. With minor exceptions, we confirm the monophyly of Poduromorpha and Symphypleona sensu stricto (the latter placed at the most basal position in the springtail phylogenetic tree), whereas monophyly of Neelipleona and Entomobryomorpha is only supported when a handful of critical taxa in these two lineages are excluded. Finally, we review gene order models observed in the class, as well as the overall mitochondrial nucleotide composition.

ACDC, a global database of amphibian cytochrome-b sequences using reproducible curation for GenBank records

Genetic data are a crucial and exponentially growing resource across all biological sciences, yet curated databases are scarce. The widespread occurrence of sequence and (meta)data errors in public repositories calls for comprehensive improvements of curation protocols leading to robust research and downstream analyses. We collated and curated all available GenBank cytochrome-b sequences for amphibians, a benchmark marker in this globally declining vertebrate clade. The Amphibia's Curated Database of Cytochrome-b (ACDC) consists of 36,514 sequences representing 2,309 species from 398 genera (median = 2 with 50% interquartile ranges of 1-7 species/genus). We updated the taxonomic identity of >4,800 sequences (ca. 13%) and found 2,359 (6%) conflicting sequences with 84% of the errors originating from taxonomic misidentifications. The database (accessible at https://doi.org/10.6084/m9.figshare.9944759 ) also includes an R script to replicate our study for other loci and taxonomic groups. We provide recommendations to improve genetic-data quality in public repositories and flag species for which there is a need for taxonomic refinement in the face of increased rate of amphibian extinctions in the Anthropocene.

Improved annotation of protein-coding genes boundaries in metazoan mitochondrial genomes

With the rapid increase of sequenced metazoan mitochondrial genomes, a detailed manual annotation is becoming more and more infeasible. While it is easy to identify the approximate location of protein-coding genes within mitogenomes, the peculiar processing of mitochondrial transcripts, however, makes the determination of precise gene boundaries a surprisingly difficult problem. We have analyzed the properties of annotated start and stop codon positions in detail, and use the inferred patterns to devise a new method for predicting gene boundaries in de novo annotations. Our method benefits from empirically observed prevalances of start/stop codons and gene lengths, and considers the dependence of these features on variations of genetic codes. Albeit not being perfect, our new approach yields a drastic improvement in the accuracy of gene boundaries and upgrades the mitochondrial genome annotation server MITOS to an even more sophisticated tool for fully automatic annotation of metazoan mitochondrial genomes.

Using multiple reference genomes to identify and resolve annotation inconsistencies

BACKGROUND

Advances in sequencing technologies have led to the release of reference genomes and annotations for multiple individuals within more well-studied systems. While each of these new genome assemblies shares significant portions of synteny between each other, the annotated structure of gene models within these regions can differ. Of particular concern are split-gene misannotations, in which a single gene is incorrectly annotated as two distinct genes or two genes are incorrectly annotated as a single gene. These misannotations can have major impacts on functional prediction, estimates of expression, and many downstream analyses.

RESULTS

We developed a high-throughput method based on pairwise comparisons of annotations that detect potential split-gene misannotations and quantifies support for whether the genes should be merged into a single gene model. We demonstrated the utility of our method using gene annotations of three reference genomes from maize (B73, PH207, and W22), a difficult system from an annotation perspective due to the size and complexity of the genome. On average, we found several hundred of these potential split-gene misannotations in each pairwise comparison, corresponding to 3-5% of gene models across annotations. To determine which state (i.e. one gene or multiple genes) is biologically supported, we utilized RNAseq data from 10 tissues throughout development along with a novel metric and simulation framework. The methods we have developed require minimal human interaction and can be applied to future assemblies to aid in annotation efforts.

CONCLUSIONS

Split-gene misannotations occur at appreciable frequency in maize annotations. We have developed a method to easily identify and correct these misannotations. Importantly, this method is generic in that it can utilize any type of short-read expression data. Failure to account for split-gene misannotations has serious consequences for biological inference, particularly for expression-based analyses.

Using multiple reference genomes to identify and resolve annotation inconsistencies

BACKGROUND

Advances in sequencing technologies have led to the release of reference genomes and annotations for multiple individuals within more well-studied systems. While each of these new genome assemblies shares significant portions of synteny between each other, the annotated structure of gene models within these regions can differ. Of particular concern are split-gene misannotations, in which a single gene is incorrectly annotated as two distinct genes or two genes are incorrectly annotated as a single gene. These misannotations can have major impacts on functional prediction, estimates of expression, and many downstream analyses.

RESULTS

We developed a high-throughput method based on pairwise comparisons of annotations that detect potential split-gene misannotations and quantifies support for whether the genes should be merged into a single gene model. We demonstrated the utility of our method using gene annotations of three reference genomes from maize (B73, PH207, and W22), a difficult system from an annotation perspective due to the size and complexity of the genome. On average, we found several hundred of these potential split-gene misannotations in each pairwise comparison, corresponding to 3-5% of gene models across annotations. To determine which state (i.e. one gene or multiple genes) is biologically supported, we utilized RNAseq data from 10 tissues throughout development along with a novel metric and simulation framework. The methods we have developed require minimal human interaction and can be applied to future assemblies to aid in annotation efforts.

CONCLUSIONS

Split-gene misannotations occur at appreciable frequency in maize annotations. We have developed a method to easily identify and correct these misannotations. Importantly, this method is generic in that it can utilize any type of short-read expression data. Failure to account for split-gene misannotations has serious consequences for biological inference, particularly for expression-based analyses.

The gene arrangement and phylogeny using mitochondrial genomes in spiders (Arachnida: Araneae)

The complete mitochondrial genome (mitogenome) of Cheiracanthium triviale was sequenced for the first time. The 14,595 bp C. triviale mitogenome contained 37 genes (13 protein coding genes, 2 ribosomal RNAs, 22 transfer RNAs) and one control region. The mitogenome of Dysdera silvatica which was available at NCBI GenBank was annotated. The mitogenome of C. triviale was compared with 43 previously sequenced spider species to observe the gene arrangements, control region and phylogeny. TreeREx analysis identified 19 mitochondrial gene rearrangements (11 transposition, 6 inversion, 2 inverse transposition) in spiders as compared with the putative ancestral gene order and lead to form new gene boundaries: trnQ-trnA, trnA-trnM for Loxosceles similis; nad3-trnS1, trnE-trnL2, trnL2-trnA, trnN-trnF for Agelena silvatica; trnN-trnE, trnE-trnA, trnR-trnF, nad4L-trnW, trnW-trnP for Carrhotus xanthogramma; trnQ-trnW, trnW-trnG, trnG-trnM for Tetragnatha nitens. Our study revealed that the gene rearrangement in spiders with putative ancestor is accelerated in Araneomorphae as compared to Mygalomorphae. Phylogenetic analysis of spiders using mitochondrial sequence data supports the monophyly of two infraorders, and sister relationship of Cheiracanthiidae with Selenopidae and Salticidae. The systematic position of the Cheiracanthium species always a controversial issue as this taxa was placed in different families by different authors.

Mechanisms of mammalian mitochondrial transcription

Numerous age-related human diseases have been associated with deficiencies in cellular energy production. Moreover, genetic alterations resulting in mitochondrial dysfunction are the cause of inheritable disorders commonly known as mitochondrial diseases. Many of these deficiencies have been directly or indirectly linked to deficits in mitochondrial gene expression. Transcription is an essential step in gene expression and elucidating the molecular mechanisms involved in this process is critical for understanding defects in energy production. For the past five decades, substantial efforts have been invested in the field of mitochondrial transcription. These efforts have led to the discovery of the main protein factors responsible for transcription as well as to a basic mechanistic understanding of the transcription process. They have also revealed various mechanisms of transcriptional regulation as well as the links that exist between the transcription process and downstream processes of RNA maturation. Here, we review the knowledge gathered in early mitochondrial transcription studies and focus on recent findings that shape our current understanding of mitochondrial transcription, posttranscriptional processing, as well as transcriptional regulation in mammalian systems.