Duplication of a truncated paralog of the family B DNA polymerase gene Aa-polB in the Agrocybe aegerita mitochondrial genome

The 203 kbp mitochondrial genome of the phytopathogenic fungus Sclerotinia borealis reveals multiple invasions of introns and genomic duplications

Here we report the complete sequence of the mitochondrial (mt) genome of the necrotrophic phytopathogenic fungus Sclerotinia borealis, a member of the order Helotiales of Ascomycetes. The 203,051 bp long mtDNA of S. borealis represents one of the largest sequenced fungal mt genomes. The large size is mostly determined by the presence of mobile genetic elements, which include 61 introns. Introns contain a total of 125,394 bp, are scattered throughout the genome, and are found in 12 protein-coding genes and in the ribosomal RNA genes. Most introns contain complete or truncated ORFs that are related to homing endonucleases of the LAGLIDADG and GIY-YIG families. Integrations of mobile elements are also evidenced by the presence of two regions similar to fragments of inverton-like plasmids. Although duplications of some short genome regions, resulting in the appearance of truncated extra copies of genes, did occur, we found no evidences of extensive accumulation of repeat sequences accounting for mitochondrial genome size expansion in some other fungi. Comparisons of mtDNA of S. borealis with other members of the order Helotiales reveal considerable gene order conservation and a dynamic pattern of intron acquisition and loss during evolution. Our data are consistent with the hypothesis that horizontal DNA transfer has played a significant role in the evolution and size expansion of the S. borealis mt genome.

Strikingly bacteria-like and gene-rich mitochondrial genomes throughout jakobid protists

The most bacteria-like mitochondrial genome known is that of the jakobid flagellate Reclinomonas americana NZ. This genome also encodes the largest known gene set among mitochondrial DNAs (mtDNAs), including the RNA subunit of RNase P (transfer RNA processing), a reduced form of transfer-messenger RNA (translational control), and a four-subunit bacteria-like RNA polymerase, which in other eukaryotes is substituted by a nucleus-encoded, single-subunit, phage-like enzyme. Further, protein-coding genes are preceded by potential Shine-Dalgarno translation initiation motifs. Whether similarly ancestral mitochondrial characters also exist in relatives of R. americana NZ is unknown. Here, we report a comparative analysis of nine mtDNAs from five distant jakobid genera: Andalucia, Histiona, Jakoba, Reclinomonas, and Seculamonas. We find that Andalucia godoyi has an even larger mtDNA gene complement than R. americana NZ. The extra genes are rpl35 (a large subunit mitoribosomal protein) and cox15 (involved in cytochrome oxidase assembly), which are nucleus encoded throughout other eukaryotes. Andalucia cox15 is strikingly similar to its homolog in the free-living α-proteobacterium Tistrella mobilis. Similarly, a long, highly conserved gene cluster in jakobid mtDNAs, which is a clear vestige of prokaryotic operons, displays a gene order more closely resembling that in free-living α-proteobacteria than in Rickettsiales species. Although jakobid mtDNAs, overall, are characterized by bacteria-like features, they also display a few remarkably divergent characters, such as 3'-tRNA editing in Seculamonas ecuadoriensis and genome linearization in Jakoba libera. Phylogenetic analysis with mtDNA-encoded proteins strongly supports monophyly of jakobids with Andalucia as the deepest divergence. However, it remains unclear which α-proteobacterial group is the closest mitochondrial relative.

Reprever: resolving low-copy duplicated sequences using template driven assembly

Genomic sequence duplication is an important mechanism for genome evolution, often resulting in large sequence variations with implications for disease progression. Although paired-end sequencing technologies are commonly used for structural variation discovery, the discovery of novel duplicated sequences remains an unmet challenge. We analyze duplicons starting from identified high-copy number variants. Given paired-end mapped reads, and a candidate high-copy region, our tool, Reprever, identifies (a) the insertion breakpoints where the extra duplicons inserted into the donor genome and (b) the actual sequence of the duplicon. Reprever resolves ambiguous mapping signatures from existing homologs, repetitive elements and sequencing errors to identify breakpoint. At each breakpoint, Reprever reconstructs the inserted sequence using profile hidden Markov model (PHMM)-based guided assembly. In a test on 1000 artificial genomes with simulated duplication, Reprever could identify novel duplicates up to 97% of genomes within 3 bp positional and 1% sequence errors. Validation on 680 fosmid sequences identified and reconstructed eight duplicated sequences with high accuracy. We applied Reprever to reanalyzing a re-sequenced data set from the African individual NA18507 to identify >800 novel duplicates, including insertions in genes and insertions with additional variation. polymerase chain reaction followed by capillary sequencing validated both the insertion locations of the strongest predictions and their predicted sequence.

Comparison of chaetognath mitochondrial genomes and phylogenetical implications

Two complete mitochondrial genomes (mtDNAs) of chaetognaths, Spadella cephaloptera and Paraspadella gotoi, have been recently published. These genomes are highly unusual. They are the smallest metazoan mtDNAs so far known; atp6 and atp8 genes are missing; lastly, our reanalysis has evidenced that, contrarily to what has been previously published for one sequence, both contain a unique transfer RNA (tRNA(Met)) evidencing that both have the same gene content. Indeed, even if the gene order seems very different, two gene blocks are conserved. In addition, comparison of gene arrangement suggests phylogenetical relationships between chaetognaths and some lophotrochozoa like annelids and molluscs.

Molecular evolution of a mitochondrial polB gene, encoding a family B DNA polymerase, towards the elimination from Agrocybe mitochondrial genomes

Two genes (Ac-polB O1and Ac-polB O2), each encoding a family B DNA polymerase, were characterized from the mitochondrial genome of the basidiomycete Agrocybe chaxingu. These two polB genes constitute orthologs of the potentially functional Aa-polB gene and its disrupted paralog Aa-polB P1, previously described in the closely related species A. aegerita. Unlike the case in Aa-polB, both gene copies in A. chaxingu are constituted by large but disrupted ORFs, which very probably encode nonfunctional enzymes: Ac-polB O1 has a deletion of 126 bp between the segments encoding the Exo II and Pol I domains and a 78-bp insertion between the Exo II and Exo III domains, whereas Ac-polB O2 has a large deletion of 1208 bp between the Exo II and Pol III domains and a deletion of 54-bp involving the 3' end of the gene. Hence, rearrangements in the Ac-polB ORFs appear to have led to their functional erosion in the mitochondrial genome in this species. Phylogenetic analysis has shown a close relationship between the mitochondrial polB genes and homologous genes carried by fungal linear plasmids, suggesting that they may have been acquired by the integration of linear plasmids into the mitochondrial genome.