Chromosome-level Genome Assembly of the Yeast Lodderomyces beijingensis Reveals the Genetic Nature of Metabolic Adaptations and Identifies Subtelomeres as Hotspots for Amplification of Mating Type Loci

Lodderomyces beijingensis is an ascosporic ascomycetous yeast. In contrast to related species Lodderomyces elongisporus, which is a recently emerging human pathogen, L. beijingensis is associated with insects. To provide an insight into its genetic makeup, we investigated the genome of its type strain, CBS 14171. We demonstrate that this yeast is diploid and describe the high contiguity nuclear genome assembly consisting of eight chromosome-sized contigs with a total size of about 15.1 Mbp. We find that the genome sequence contains multiple copies of the mating type loci and codes for essential components of the mating pheromone response pathway, however, the missing orthologs of several genes involved in the meiotic program raise questions about the mode of sexual reproduction. We also show that L. beijingensis genome codes for the 3-oxoadipate pathway enzymes, which allow the assimilation of protocatechuate. In contrast, the GAL gene cluster underwent a decay resulting in an inability of L. beijingensis to utilize galactose. Moreover, we find that the 56.5 kbp long mitochondrial DNA is structurally similar to known linear mitochondrial genomes terminating on both sides with covalently closed single-stranded hairpins. Finally, we discovered a new double-stranded RNA mycovirus from the Totiviridae family and characterize its genome sequence.

y-mtPTM: Yeast Mitochondrial Post-Translational Modification Database

One powerful strategy of how to increase the complexity of cellular proteomes is through post-translational modifications (PTMs) of proteins. Currently there are approximately 400 types of PTMs, the different combinations of which yield a large variety of protein isoforms with distinct biochemical properties. Although mitochondrial proteins undergoing PTMs were identified nearly six decades ago, studies on the roles and extent of PTMs on mitochondrial functions lagged behind the other cellular compartments. The application of mass spectrometry for the characterization of the mitochondrial proteome as well as for the detection of various PTMs resulted in the identification of thousands of amino acid positions that can be modified by different chemical groups. However, the data on mitochondrial PTMs are scattered in several datasets, and the available databases do not contain a complete list of modified residues. To integrate information on PTMs of the mitochondrial proteome of the yeast Saccharomyces cerevisiae, we built the y-mtPTM (Yeast Mitochondrial Post-Translational Modification) Database (http://compbio.fmph.uniba.sk/y-mtptm/). It lists nearly 20,000 positions on mitochondrial proteins affected by approximately twenty various PTMs, with phosphorylated, succinylated, acetylated and ubiquitylated sites being the most abundant. A simple search of a protein of interest reveals the modified amino acid residues, their position within the primary sequence as well as on its 3D structure and links to the source reference(s). The database will serve yeast mitochondrial researchers as a comprehensive platform to investigate the functional significance of the PTMs of mitochondrial proteins.

Transcriptome and proteome profiling reveals complex adaptations of Candida parapsilosis cells assimilating hydroxyaromatic carbon sources

Many fungal species utilize hydroxyderivatives of benzene and benzoic acid as carbon sources. The yeast Candida parapsilosis metabolizes these compounds via the 3-oxoadipate and gentisate pathways, whose components are encoded by two metabolic gene clusters. In this study, we determine the chromosome level assembly of the C. parapsilosis strain CLIB214 and use it for transcriptomic and proteomic investigation of cells cultivated on hydroxyaromatic substrates. We demonstrate that the genes coding for enzymes and plasma membrane transporters involved in the 3-oxoadipate and gentisate pathways are highly upregulated and their expression is controlled in a substrate-specific manner. However, regulatory proteins involved in this process are not known. Using the knockout mutants, we show that putative transcriptional factors encoded by the genes OTF1 and GTF1 located within these gene clusters function as transcriptional activators of the 3-oxoadipate and gentisate pathway, respectively. We also show that the activation of both pathways is accompanied by upregulation of genes for the enzymes involved in β-oxidation of fatty acids, glyoxylate cycle, amino acid metabolism, and peroxisome biogenesis. Transcriptome and proteome profiles of the cells grown on 4-hydroxybenzoate and 3-hydroxybenzoate, which are metabolized via the 3-oxoadipate and gentisate pathway, respectively, reflect their different connection to central metabolism. Yet we find that the expression profiles differ also in the cells assimilating 4-hydroxybenzoate and hydroquinone, which are both metabolized in the same pathway. This finding is consistent with the phenotype of the Otf1p-lacking mutant, which exhibits impaired growth on hydroxybenzoates, but still utilizes hydroxybenzenes, thus indicating that additional, yet unidentified transcription factor could be involved in the 3-oxoadipate pathway regulation. Moreover, we propose that bicarbonate ions resulting from decarboxylation of hydroxybenzoates also contribute to differences in the cell responses to hydroxybenzoates and hydroxybenzenes. Finally, our phylogenetic analysis highlights evolutionary paths leading to metabolic adaptations of yeast cells assimilating hydroxyaromatic substrates.

Genome analysis of five recently described species of the CUG-Ser clade uncovers Candida theae as a new hybrid lineage with pathogenic potential in the Candida parapsilosis species complex

Candida parapsilosis species complex comprises three important pathogenic species: Candida parapsilosis sensu stricto, Candida orthopsilosis and Candida metapsilosis. The majority of C. orthopsilosis and all C. metapsilosis isolates sequenced thus far are hybrids, and most of the parental lineages remain unidentified. This led to the hypothesis that hybrids with pathogenic potential were formed by the hybridization of non-pathogenic lineages that thrive in the environment. In a search for the missing hybrid parentals, and aiming to get a better understanding of the evolution of the species complex, we sequenced, assembled and analysed the genome of five close relatives isolated from the environment: Candida jiufengensis, Candida pseudojiufengensis, Candida oxycetoniae, Candida margitis and Candida theae. We found that the linear conformation of mitochondrial genomes in Candida species emerged multiple times independently. Furthermore, our analyses discarded the possible involvement of these species in the mentioned hybridizations, but identified C. theae as an additional hybrid in the species complex. Importantly, C. theae was recently associated with a case of infection, and we also uncovered the hybrid nature of this clinical isolate. Altogether, our results reinforce the hypothesis that hybridization is widespread among Candida species, and potentially contributes to the emergence of lineages with opportunistic pathogenic behaviour.

Monoclonal antibodies targeting two immunodominant epitopes on the Spike protein neutralize emerging SARS-CoV-2 variants of concern

BACKGROUND

The emergence of new SARS-CoV-2 variants of concern B.1.1.7 (Alpha), B.1.351 (Beta), P.1 (Gamma) and B.1.617.2 (Delta) that harbor mutations in the viral S protein raised concern about activity of current vaccines and therapeutic antibodies. Independent studies have shown that mutant variants are partially or completely resistant against some of the therapeutic antibodies authorized for emergency use.

METHODS

We employed hybridoma technology, ELISA-based and cell-based S-ACE2 interaction assays combined with authentic virus neutralization assays to develop second-generation antibodies, which were specifically selected for their ability to neutralize the new variants of SARS-CoV-2.

FINDINGS

AX290 and AX677, two monoclonal antibodies with non-overlapping epitopes, exhibit subnanomolar or nanomolar affinities to the receptor binding domain of the viral Spike protein carrying amino acid substitutions N501Y, N439K, E484K, K417N, and a combination N501Y/E484K/K417N found in the circulating virus variants. The antibodies showed excellent neutralization of an authentic SARS-CoV-2 virus representing strains circulating in Europe in spring 2020 and also the variants of concern B.1.1.7 (Alpha), B.1.351 (Beta) and B.1.617.2 (Delta). In addition, AX677 is able to bind Omicron Spike protein just like the wild type Spike. The combination of the two antibodies prevented the appearance of escape mutations of the authentic SARS-CoV-2 virus. Prophylactic administration of AX290 and AX677, either individually or in combination, effectively reduced viral burden and inflammation in the lungs, and prevented disease in a mouse model of SARS-CoV-2 infection.

INTERPRETATION

The virus-neutralizing properties were fully reproduced in chimeric mouse-human versions of the antibodies, which may represent a promising tool for COVID-19 therapy.

FUNDING

The study was funded by AXON Neuroscience SE and AXON COVIDAX a.s.

Nanopore sequencing of SARS-CoV-2: Comparison of short and long PCR-tiling amplicon protocols

Surveillance of the SARS-CoV-2 variants including the quickly spreading mutants by rapid and near real-time sequencing of the viral genome provides an important tool for effective health policy decision making in the ongoing COVID-19 pandemic. Here we evaluated PCR-tiling of short (~400-bp) and long (~2 and ~2.5-kb) amplicons combined with nanopore sequencing on a MinION device for analysis of the SARS-CoV-2 genome sequences. Analysis of several sequencing runs demonstrated that using the long amplicon schemes outperforms the original protocol based on the 400-bp amplicons. It also illustrated common artefacts and problems associated with PCR-tiling approach, such as uneven genome coverage, variable fraction of discarded sequencing reads, including human and bacterial contamination, as well as the presence of reads derived from the viral sub-genomic RNAs.

Surveillance of SARS-CoV-2 lineage B.1.1.7 in Slovakia using a novel, multiplexed RT-qPCR assay

The emergence of a novel SARS-CoV-2 B.1.1.7 variant sparked global alarm due to increased transmissibility, mortality, and uncertainty about vaccine efficacy, thus accelerating efforts to detect and track the variant. Current approaches to detect B.1.1.7 include sequencing and RT-qPCR tests containing a target assay that fails or results in reduced sensitivity towards the B.1.1.7 variant. Since many countries lack genomic surveillance programs and failed assays detect unrelated variants containing similar mutations as B.1.1.7, we used allele-specific PCR, and judicious placement of LNA-modified nucleotides to develop an RT-qPCR test that accurately and rapidly differentiates B.1.1.7 from other SARS-CoV-2 variants. We validated the test on 106 clinical samples with lineage status confirmed by sequencing and conducted a country-wide surveillance study of B.1.1.7 prevalence in Slovakia. Our multiplexed RT-qPCR test showed 97% clinical sensitivity and retesting 6,886 SARS-CoV-2 positive samples obtained during three campaigns performed within one month, revealed pervasive spread of B.1.1.7 with an average prevalence of 82%. Labs can easily implement this test to rapidly scale B.1.1.7 surveillance efforts and it is particularly useful in countries with high prevalence of variants possessing only the ΔH69/ΔV70 deletion because current strategies using target failure assays incorrectly identify these as putative B.1.1.7 variants.

A SARS-CoV-2 mutant from B.1.258 lineage with ∆H69/∆V70 deletion in the Spike protein circulating in Central Europe in the fall 2020

SARS-CoV-2 mutants carrying the ∆H69/∆V70 deletion in the amino-terminal domain of the Spike protein emerged independently in at least six lineages of the virus (namely, B.1.1.7, B.1.1.298, B.1.160, B.1.177, B.1.258, B.1.375). We analyzed SARS-CoV-2 samples collected from various regions of Slovakia between November and December 2020 that were presumed to contain B.1.1.7 variant due to drop-out of the Spike gene target in an RT-qPCR test caused by this deletion. Sequencing of these samples revealed that although in some cases the samples were indeed confirmed as B.1.1.7, a substantial fraction of samples contained another ∆H69/∆V70 carrying mutant belonging to the lineage B.1.258, which has been circulating in Central Europe since August 2020, long before the import of B.1.1.7. Phylogenetic analysis shows that the early sublineage of B.1.258 acquired the N439K substitution in the receptor-binding domain (RBD) of the Spike protein and, later on, also the deletion ∆H69/∆V70 in the Spike N-terminal domain (NTD). This variant was particularly common in several European countries including the Czech Republic and Slovakia but has been quickly replaced by B.1.1.7 early in 2021.

Genome analysis of Candida subhashii reveals its hybrid nature and dual mitochondrial genome conformations

Candida subhashii belongs to the CUG-Ser clade, a group of phylogenetically closely related yeast species that includes some human opportunistic pathogens, such as Candida albicans. Despite being present in the environment, C. subhashii was initially described as the causative agent of a case of peritonitis. Considering the relevance of whole-genome sequencing and analysis for our understanding of genome evolution and pathogenicity, we sequenced, assembled and annotated the genome of C. subhashii type strain. Our results show that C. subhashii presents a highly heterozygous genome and other signatures that point to a hybrid ancestry. The presence of functional pathways for assimilation of hydroxyaromatic compounds goes in line with the affiliation of this yeast with soil microbial communities involved in lignin decomposition. Furthermore, we observed that different clones of this strain may present circular or linear mitochondrial DNA. Re-sequencing and comparison of strains with differential mitochondrial genome topology revealed five candidate genes potentially associated with this conformational change: MSK1, SSZ1, ALG5, MRPL9 and OYE32.

OCT1 - a yeast mitochondrial thiolase involved in the 3-oxoadipate pathway

The 3-oxoacyl-CoA thiolases catalyze the last step of the fatty acid β-oxidation pathway. In yeasts and plants, this pathway takes place exclusively in peroxisomes, whereas in animals it occurs in both peroxisomes and mitochondria. In contrast to baker's yeast Saccharomyces cerevisiae, yeast species from the Debaryomycetaceae family also encode a thiolase with predicted mitochondrial localization. These yeasts are able to utilize a range of hydroxyaromatic compounds via the 3-oxoadipate pathway the last step of which is catalyzed by 3-oxoadipyl-CoA thiolase and presumably occurs in mitochondria. In this work, we studied Oct1p, an ortholog of this enzyme from Candida parapsilosis. We found that the cells grown on a 3-oxoadipate pathway substrate exhibit increased levels of the OCT1 mRNA. Deletion of both OCT1 alleles impairs the growth of C. parapsilosis cells on 3-oxoadipate pathway substrates and this defect can be rescued by expression of the OCT1 gene from a plasmid vector. Subcellular localization experiments and LC-MS/MS analysis of enriched organellar fraction-proteins confirmed the presence of Oct1p in mitochondria. Phylogenetic profiling of Oct1p revealed an intricate evolutionary pattern indicating multiple horizontal gene transfers among different fungal groups.

Shared evolutionary footprints suggest mitochondrial oxidative damage underlies multiple complex I losses in fungi

Oxidative phosphorylation is among the most conserved mitochondrial pathways. However, one of the cornerstones of this pathway, the multi-protein complex NADH : ubiquinone oxidoreductase (complex I) has been lost multiple independent times in diverse eukaryotic lineages. The causes and consequences of these convergent losses remain poorly understood. Here, we used a comparative genomics approach to reconstruct evolutionary paths leading to complex I loss and infer possible evolutionary scenarios. By mining available mitochondrial and nuclear genomes, we identified eight independent events of mitochondrial complex I loss across eukaryotes, of which six occurred in fungal lineages. We focused on three recent loss events that affect closely related fungal species, and inferred genomic changes convergently associated with complex I loss. Based on these results, we predict novel complex I functional partners and relate the loss of complex I with the presence of increased mitochondrial antioxidants, higher fermentative capabilities, duplications of alternative dehydrogenases, loss of alternative oxidases and adaptation to antifungal compounds. To explain these findings, we hypothesize that a combination of previously acquired compensatory mechanisms and exposure to environmental triggers of oxidative stress (such as hypoxia and/or toxic chemicals) induced complex I loss in fungi.

Design and Synthesis of Pyrano[3,2-b]indolones Showing Antimycobacterial Activity

Latent Mycobacterium tuberculosis infection presents one of the largest challenges for tuberculosis control and novel antimycobacterial drug development. A series of pyrano[3,2-b]indolone-based compounds was designed and synthesized via an original eight-step scheme. The synthesized compounds were evaluated for their in vitro activity against M. tuberculosis strains H37Rv and streptomycin-starved 18b (SS18b), representing models for replicating and nonreplicating mycobacteria, respectively. Compound 10a exhibited good activity with MIC₉₉ values of 0.3 and 0.4 μg/mL against H37Rv and SS18b, respectively, as well as low toxicity, acceptable intracellular activity, and satisfactory metabolic stability and was selected as the lead compound for further studies. An analysis of 10a-resistant M. bovis mutants disclosed a cross-resistance with pretomanid and altered relative amounts of different forms of cofactor F₄₂₀ in these strains. Complementation experiments showed that F₄₂₀-dependent glucose-6-phosphate dehydrogenase and the synthesis of mature F₄₂₀ were important for 10a activity. Overall these studies revealed 10a to be a prodrug that is activated by an unknown F₄₂₀-dependent enzyme in mycobacteria.

Mitochondrial protein phosphorylation in yeast revisited

Protein phosphorylation is one of the best-known post-translational modifications occurring in all domains of life. In eukaryotes, protein phosphorylation affects all cellular compartments including mitochondria. High-throughput techniques of mass spectrometry combined with cell fractionation and biochemical methods yielded thousands of phospho-sites on hundreds of mitochondrial proteins. We have compiled the information on mitochondrial protein kinases and phosphatases and their substrates in Saccharomyces cerevisiae and provide the current state-of-the-art overview of mitochondrial protein phosphorylation in this model eukaryote. Using several examples, we describe emerging features of the yeast mitochondrial phosphoproteome and present challenges lying ahead in this exciting field.

Step-by-Step Evolution of Telomeres: Lessons from Yeasts

In virtually every eukaryotic species, the ends of nuclear chromosomes are protected by telomeres, nucleoprotein structures counteracting the end-replication problem and suppressing recombination and undue DNA repair. Although in most cases, the primary structure of telomeric DNA is conserved, there are several exceptions to this rule. One is represented by the telomeric repeats of ascomycetous yeasts, which encompass a great variety of sequences, whose evolutionary origin has been puzzling for several decades. At present, the key questions concerning the driving force behind their rapid evolution and the means of co-evolution of telomeric repeats and telomere-binding proteins remain largely unanswered. Previously published studies addressed mostly the general concepts of the evolutionary origin of telomeres, key properties of telomeric proteins as well as the molecular mechanisms of telomere maintenance; however, the evolutionary process itself has not been analyzed thoroughly. Here, we aimed to inspect the evolution of telomeres in ascomycetous yeasts from the subphyla Saccharomycotina and Taphrinomycotina, with special focus on the evolutionary origin of species-specific telomeric repeats. We analyzed the sequences of telomeric repeats from 204 yeast species classified into 20 families and as a result, we propose a step-by-step model, which integrates the diversity of telomeric repeats, telomerase RNAs, telomere-binding protein complexes and explains a propensity of certain species to generate the repeat heterogeneity within a single telomeric array.

Mitochondrial HMG-Box Containing Proteins: From Biochemical Properties to the Roles in Human Diseases

Mitochondrial DNA (mtDNA) molecules are packaged into compact nucleo-protein structures called mitochondrial nucleoids (mt-nucleoids). Their compaction is mediated in part by high-mobility group (HMG)-box containing proteins (mtHMG proteins), whose additional roles include the protection of mtDNA against damage, the regulation of gene expression and the segregation of mtDNA into daughter organelles. The molecular mechanisms underlying these functions have been identified through extensive biochemical, genetic, and structural studies, particularly on yeast (Abf2) and mammalian mitochondrial transcription factor A (TFAM) mtHMG proteins. The aim of this paper is to provide a comprehensive overview of the biochemical properties of mtHMG proteins, the structural basis of their interaction with DNA, their roles in various mtDNA transactions, and the evolutionary trajectories leading to their rapid diversification. We also describe how defects in the maintenance of mtDNA in cells with dysfunctional mtHMG proteins lead to different pathologies at the cellular and organismal level.

Role of folding kinetics of secondary structures in telomeric G-overhangs in the regulation of telomere maintenance in Saccharomyces cerevisiae

The ends of eukaryotic chromosomes typically contain a 3' ssDNA G-rich protrusion (G-overhang). This overhang must be protected against detrimental activities of nucleases and of the DNA damage response machinery and participates in the regulation of telomerase, a ribonucleoprotein complex that maintains telomere integrity. These functions are mediated by DNA-binding proteins, such as Cdc13 in Saccharomyces cerevisiae, and the propensity of G-rich sequences to form various non-B DNA structures. Using CD and NMR spectroscopies, we show here that G-overhangs of S. cerevisiae form distinct Hoogsteen pairing-based secondary structures, depending on their length. Whereas short telomeric oligonucleotides form a G-hairpin, their longer counterparts form parallel and/or antiparallel G-quadruplexes (G4s). Regardless of their topologies, non-B DNA structures exhibited impaired binding to Cdc13 in vitro as demonstrated by electrophoretic mobility shift assays. Importantly, whereas G4 structures formed relatively quickly, G-hairpins folded extremely slowly, indicating that short G-overhangs, which are typical for most of the cell cycle, are present predominantly as single-stranded oligonucleotides and are suitable substrates for Cdc13. Using ChIP, we show that the occurrence of G4 structures peaks at the late S phase, thus correlating with the accumulation of long G-overhangs. We present a model of how time- and length-dependent formation of non-B DNA structures at chromosomal termini participates in telomere maintenance.

Co-evolution in the Jungle: From Leafcutter Ant Colonies to Chromosomal Ends

Biological entities are multicomponent systems where each part is directly or indirectly dependent on the others. In effect, a change in a single component might have a consequence on the functioning of its partners, thus affecting the fitness of the entire system. In this article, we provide a few examples of such complex biological systems, ranging from ant colonies to a population of amino acids within a single-polypeptide chain. Based on these examples, we discuss one of the central and still challenging questions in biology: how do such multicomponent consortia co-evolve? More specifically, we ask how telomeres, nucleo-protein complexes protecting the integrity of linear DNA chromosomes, originated from the ancestral organisms having circular genomes and thus not dealing with end-replication and end-protection problems. Using the examples of rapidly evolving topologies of mitochondrial genomes in eukaryotic microorganisms, we show what means of co-evolution were employed to accommodate various types of telomere-maintenance mechanisms in mitochondria. We also describe an unprecedented runaway evolution of telomeric repeats in nuclei of ascomycetous yeasts accompanied by co-evolution of telomere-associated proteins. We propose several scenarios derived from research on telomeres and supported by other studies from various fields of biology, while emphasizing that the relevant answers are still not in sight. It is this uncertainty and a lack of a detailed roadmap that makes the journey through the jungle of biological systems still exciting and worth undertaking.

Identification of telomerase RNAs in species of the Yarrowia clade provides insights into the co-evolution of telomerase, telomeric repeats and telomere-binding proteins

Telomeric repeats in fungi of the subphylum Saccharomycotina exhibit great inter- and intra-species variability in length and sequence. Such variations challenged telomeric DNA-binding proteins that co-evolved to maintain their functions at telomeres. Here, we compare the extent of co-variations in telomeric repeats, encoded in the telomerase RNAs (TERs), and the repeat-binding proteins from 13 species belonging to the Yarrowia clade. We identified putative TER loci, analyzed their sequence and secondary structure conservation, and predicted functional elements. Moreover, in vivo complementation assays with mutant TERs showed the functional importance of four novel TER substructures. The TER-derived telomeric repeat unit of all species, except for one, is 10 bp long and can be represented as 5′-TTNNNNAGGG-3′, with repeat sequence variations occuring primarily outside the vertebrate telomeric motif 5′-TTAGGG-3′. All species possess a homologue of the Yarrowia lipolytica Tay1 protein, YlTay1p. In vitro, YlTay1p displays comparable DNA-binding affinity to all repeat variants, suggesting a conserved role among these species. Taken together, these results add significant insights into the co-evolution of TERs, telomeric repeats and telomere-binding proteins in yeasts.

A New View of the T-Loop Junction: Implications for Self-Primed Telomere Extension, Expansion of Disease-Related Nucleotide Repeat Blocks, and Telomere Evolution

Telomere loops (t-loops) are formed at the ends of chromosomes in species ranging from humans to worms, plants, and with genetic manipulation, some yeast. Recent in vitro studies demonstrated that transcription of telomeric DNA leads to highly efficient t-loop formation. It was also shown that both DNA termini are inserted into the preceding DNA to generate a highly stable t-loop junction. Furthermore, some telomeric RNA remains present at the junction, potentially acting as a plug to further protect and stabilize the t-loop. Modeling the loop junction reveals two mechanisms by which the canonical chromosomal replication factors could extend the telomere in the absence of telomerase. One mechanism would utilize the annealed 3’ terminus as a de novo replication origin. In vitro evidence for the ability of the t-loop to prime telomere extension using the T7 replication factors is presented. A second mechanism would involve resolution of the Holliday junction present in the t-loop bubble by factors such as GEN1 to generate a rolling circle template at the extreme terminus of the telomere. This could lead to large expansions of the telomeric tract. Here, we propose that telomeres evolved as terminal elements containing long arrays of short nucleotide repeats due to the ability of such arrays to fold back into loops and self-prime their replicative extension. In this view, telomerase may have evolved later to provide a more precise mechanism of telomere maintenance. Both pathways have direct relevance to the alternative lengthening of telomeres (ALT) pathway. This view also provides a possible mechanism for the very large repeat expansions observed in nucleotide repeat diseases such as Fragile X syndrome, myotonic dystrophy, familial amyotrophic lateral sclerosis (ALS), and frontotemporal dementia (FTD). The evolution of telomeres is discussed in the framework of these models.

Rare causes of acute appendicitis

INTRODUCTION

A clinical diagnosis of acute appendicitis does not necessarily imply only inflammatory changes found upon histopathological examination. Even though appendiceal tumors are very rare, accounting for less than 4% of all gastrointestinal tumors‚ they need to be taken into account and require the right choices regarding subsequent onco-surgical management.

METHODS

1036 appendectomies were performed and 38 (3.7%) rare appendiceal tumor cases were identified at the department of Surgery in Pilsen from 1 January 2008 to 1 October 2018. The following text presents the difficult nature of the initial diagnosis, current therapy and principles of follow-up care in this relatively rare group of patients.

CONCLUSION

Appendiceal tumors are very rare. Depending on the histopathological type of tumor, it is necessary to choose adequate therapy and also provide proper follow up care.

Genome sequence of the opportunistic human pathogen Magnusiomyces capitatus

The yeast Magnusiomyces capitatus is an opportunistic human pathogen causing rare yet severe infections, especially in patients with hematological malignancies. Here, we report the 20.2 megabase genome sequence of an environmental strain of this species as well as the genome sequences of eight additional isolates from human and animal sources providing an insight into intraspecies variation. The distribution of single-nucleotide variants is indicative of genetic recombination events, supporting evidence for sexual reproduction in this heterothallic yeast. Using RNAseq-aided annotation, we identified genes for 6518 proteins including several expanded families such as kexin proteases and Hsp70 molecular chaperones. Several of these families are potentially associated with the ability of M. capitatus to infect and colonize humans. For the purpose of comparative analysis, we also determined the genome sequence of a closely related yeast, Magnusiomyces ingens. The genome sequences of M. capitatus and M. ingens exhibit many distinct features and represent a basis for further comparative and functional studies.

Eukaryotic transporters for hydroxyderivatives of benzoic acid

Several yeast species catabolize hydroxyderivatives of benzoic acid. However, the nature of carriers responsible for transport of these compounds across the plasma membrane is currently unknown. In this study, we analyzed a family of genes coding for permeases belonging to the major facilitator superfamily (MFS) in the pathogenic yeast Candida parapsilosis. Our results revealed that these transporters are functionally equivalent to bacterial aromatic acid: H⁺ symporters (AAHS) such as GenK, MhbT and PcaK. We demonstrate that the genes HBT1 and HBT2 encoding putative transporters are highly upregulated in C. parapsilosis cells assimilating hydroxybenzoate substrates and the corresponding proteins reside in the plasma membrane. Phenotypic analyses of knockout mutants and hydroxybenzoate uptake assays provide compelling evidence that the permeases Hbt1 and Hbt2 transport the substrates that are metabolized via the gentisate (3-hydroxybenzoate, gentisate) and 3-oxoadipate pathway (4-hydroxybenzoate, 2,4-dihydroxybenzoate and protocatechuate), respectively. Our data support the hypothesis that the carriers belong to the AAHS family of MFS transporters. Phylogenetic analyses revealed that the orthologs of Hbt permeases are widespread in the subphylum Pezizomycotina, but have a sparse distribution among Saccharomycotina lineages. Moreover, these analyses shed additional light on the evolution of biochemical pathways involved in the catabolic degradation of hydroxyaromatic compounds.

Double-stranded telomeric DNA binding proteins: Diversity matters

Telomeric sequences constitute only a small fraction of the whole genome yet they are crucial for ensuring genomic stability. This function is in large part mediated by protein complexes recruited to telomeric sequences by specific telomere-binding proteins (TBPs). Although the principal tasks of nuclear telomeres are the same in all eukaryotes, TBPs in various taxa exhibit a surprising diversity indicating their distinct evolutionary origin. This diversity is especially pronounced in ascomycetous yeasts where they must have co-evolved with rapidly diversifying sequences of telomeric repeats. In this article we (i) provide a historical overview of the discoveries leading to the current list of TBPs binding to double-stranded (ds) regions of telomeres, (ii) describe examples of dsTBPs highlighting their diversity in even closely related species, and (iii) speculate about possible evolutionary trajectories leading to a long list of various dsTBPs fulfilling the same general role(s) in their own unique ways.

Mitochondrial Carriers Link the Catabolism of Hydroxyaromatic Compounds to the Central Metabolism in Candida parapsilosis

The pathogenic yeast Candida parapsilosis metabolizes hydroxyderivatives of benzene and benzoic acid to compounds channeled into central metabolism, including the mitochondrially localized tricarboxylic acid cycle, via the 3-oxoadipate and gentisate pathways. The orchestration of both catabolic pathways with mitochondrial metabolism as well as their evolutionary origin is not fully understood. Our results show that the enzymes involved in these two pathways operate in the cytoplasm with the exception of the mitochondrially targeted 3-oxoadipate CoA-transferase (Osc1p) and 3-oxoadipyl-CoA thiolase (Oct1p) catalyzing the last two reactions of the 3-oxoadipate pathway. The cellular localization of the enzymes indicates that degradation of hydroxyaromatic compounds requires a shuttling of intermediates, cofactors, and products of the corresponding biochemical reactions between cytosol and mitochondria. Indeed, we found that yeast cells assimilating hydroxybenzoates increase the expression of genes SFC1, LEU5, YHM2, and MPC1 coding for succinate/fumarate carrier, coenzyme A carrier, oxoglutarate/citrate carrier, and the subunit of pyruvate carrier, respectively. A phylogenetic analysis uncovered distinct evolutionary trajectories for sparsely distributed gene clusters coding for enzymes of both pathways. Whereas the 3-oxoadipate pathway appears to have evolved by vertical descent combined with multiple losses, the gentisate pathway shows a striking pattern suggestive of horizontal gene transfer to the evolutionarily distant Mucorales.

Mgm101: A double-duty Rad52-like protein

Mgm101 has well-characterized activity for the repair and replication of the mitochondrial genome. Recent work has demonstrated a further role for Mgm101 in nuclear DNA metabolism, contributing to an S-phase specific DNA interstrand cross-link repair pathway that acts redundantly with a pathway controlled by Pso2 exonuclease. Due to involvement of FANCM, FANCJ and FANCP homologues (Mph1, Chl1 and Slx4), this pathway has been described as a Fanconi anemia-like pathway. In this pathway, Mgm101 physically interacts with the DNA helicase Mph1 and the MutSα (Msh2/Msh6) heterodimer, but its precise role is yet to be elucidated. Data presented here suggests that Mgm101 functionally overlaps with Rad52, supporting previous suggestions that, based on protein structure and biochemical properties, Mgm101 and Rad52 belong to a family of proteins with similar function. In addition, our data shows that this overlap extends to the function of both proteins at telomeres, where Mgm101 is required for telomere elongation during chromosome replication in rad52 defective cells. We hypothesize that Mgm101 could, in Rad52-like manner, preferentially bind single-stranded DNAs (such as at stalled replication forks, broken chromosomes and natural chromosome ends), stabilize them and mediate single-strand annealing-like homologous recombination event to prevent them from converting into toxic structures.

Correction: The Genomic Aftermath of Hybridization in the Opportunistic Pathogen Candida metapsilosis

[This corrects the article DOI: 10.1371/journal.pgen.1005626.].

Evolution of Telomeres in Schizosaccharomyces pombe and Its Possible Relationship to the Diversification of Telomere Binding Proteins

Telomeres of nuclear chromosomes are usually composed of an array of tandemly repeated sequences that are recognized by specific Myb domain containing DNA-binding proteins (telomere-binding proteins, TBPs). Whereas in many eukaryotes the length and sequence of the telomeric repeat is relatively conserved, telomeric sequences in various yeasts are highly variable. Schizosaccharomyces pombe provides an excellent model for investigation of co-evolution of telomeres and TBPs. First, telomeric repeats of S. pombe differ from the canonical mammalian type TTAGGG sequence. Second, S. pombe telomeres exhibit a high degree of intratelomeric heterogeneity. Third, S. pombe contains all types of known TBPs (Rap1p [a version unable to bind DNA], Tay1p/Teb1p, and Taz1p) that are employed by various yeast species to protect their telomeres. With the aim of reconstructing evolutionary paths leading to a separation of roles between Teb1p and Taz1p, we performed a comparative analysis of the DNA-binding properties of both proteins using combined qualitative and quantitative biochemical approaches. Visualization of DNA-protein complexes by electron microscopy revealed qualitative differences of binding of Teb1p and Taz1p to mammalian type and fission yeast telomeres. Fluorescence anisotropy analysis quantified the binding affinity of Teb1p and Taz1p to three different DNA substrates. Additionally, we carried out electrophoretic mobility shift assays using mammalian type telomeres and native substrates (telomeric repeats, histone-box sequences) as well as their mutated versions. We observed relative DNA sequence binding flexibility of Taz1p and higher binding stringency of Teb1p when both proteins were compared directly to each other. These properties may have driven replacement of Teb1p by Taz1p as the TBP in fission yeast.

The structure and DNA-binding properties of Mgm101 from a yeast with a linear mitochondrial genome

To study the mechanisms involved in the maintenance of a linear mitochondrial genome we investigated the biochemical properties of the recombination protein Mgm101 from Candida parapsilosis. We show that CpMgm101 complements defects associated with the Saccharomyces cerevisiae mgm101-1(ts) mutation and that it is present in both the nucleus and mitochondrial nucleoids of C. parapsilosis. Unlike its S. cerevisiae counterpart, CpMgm101 is associated with the entire nucleoid population and is able to bind to a broad range of DNA substrates in a non-sequence specific manner. CpMgm101 is also able to catalyze strand annealing and D-loop formation. CpMgm101 forms a roughly C-shaped trimer in solution according to SAXS. Electron microscopy of a complex of CpMgm101 with a model mitochondrial telomere revealed homogeneous, ring-shaped structures at the telomeric single-stranded overhangs. The DNA-binding properties of CpMgm101, together with its DNA recombination properties, suggest that it can play a number of possible roles in the replication of the mitochondrial genome and the maintenance of its telomeres.

[Coincidence of colonic lymphoma and gallstone ileus - case report]

INTRODUCTION

Primary colonic lymphoma is a very rare malignant disease of the gastrointestinal tract, accounting for 14% of all malignant diseases in this location. It is classified in the group of extranodal lymphomas; its long-term asymptomatic progression makes it different from common colorectal carcinomas making its diagnosis very difficult, more often accidental. Gallstone ileus is quite an uncommon complication of cholecystolithiasis diagnosed with difficulty. Up to 50% of cases are diagnosed during surgery. The obturated location depends on the size of the stone, location of the conjunction between the biliary and gastrointestinal tracts, and also on any preexisting stenosis due to another unknown pathology.

CASE REPORT

We present a case of an 86-year-old man treated for acute diverticulitis with typical clinical symptoms. Following further examination (colonoscopy, computed tomography) revealed a tumour-like infiltration in the sigmoid colon wall and a voluminous polyp was suspected according to the colonoscopy. Computed tomography described an obstruction by a biliary stone tumbling through the cholecystocolonic fistula. Subsequent biopsy supported the suspected malignant etiology. The patient underwent resection of the sigmoid colon sec. Hartmann; an infiltration was found in the subhepatic space, which corresponded to the described fistulisation between the biliary tract and the colon. A large 40 mm gallstone was found in the resected sigmoid colon over the stenosis and the bowel wall showed diffuse thickening with several polyps; final histopathological assessment confirmed malignant lymphoma of the plasmocytoma type. No serious complications occurred in the postoperative period; after healing, the patient was transferred to hematooncology care.

CONCLUSION

The article describes the presence of two rare diseases - colonic lymphoma and gallstone ileus. Clearly, without the biliary stone obstruction in the preexisting tumorous stenosis in the sigmoid colon, the malignant hematooncology disease would not have been diagnosed.

KEY WORDS

primary colonic lymphoma - gallstone ileus - complication of the cholecystolithiasis - extranodal lymphoma - acute diverticulitis.

The Genomic Aftermath of Hybridization in the Opportunistic Pathogen Candida metapsilosis

Candida metapsilosis is a rarely-isolated, opportunistic pathogen that belongs to a clade of pathogenic yeasts known as the C. parapsilosis sensu lato species complex. To gain insight into the recent evolution of C. metapsilosis and the genetic basis of its virulence, we sequenced the genome of 11 clinical isolates from various locations, which we compared to each other and to the available genomes of the two remaining members of the complex: C. orthopsilosis and C. parapsilosis. Unexpectedly, we found compelling genomic evidence that C. metapsilosis is a highly heterozygous hybrid species, with all sequenced clinical strains resulting from the same past hybridization event involving two parental lineages that were approximately 4.5% divergent in sequence. This result indicates that the parental species are non-pathogenic, but that hybridization between them formed a new opportunistic pathogen, C. metapsilosis, that has achieved a worldwide distribution. We show that these hybrids are diploid and we identified strains carrying loci for both alternative mating types, which supports mating as the initial mechanism for hybrid formation. We trace the aftermath of this hybridization at the genomic level, and reconstruct the evolutionary relationships among the different strains. Recombination and introgression -resulting in loss of heterozygosis- between the two subgenomes have been rampant, and includes the partial overwriting of the MTLa mating locus in all strains. Collectively, our results shed light on the recent genomic evolution within the C. parapsilosis sensu lato complex, and argue for a re-definition of species within this clade, with at least five distinct homozygous lineages, some of which having the ability to form hybrids.

[Intestinal volvulus caused by the ingestion of magnet balls: unexpected risk in children]

INTRODUCTION

The occurrence of swallowed foreign bodies in the digestive system is a common problem in children with the highest incidence in children aged six months to five years. Most swallowed objects leave the human body per vias naturales while 10-20% of swallowed foreign bodies need to be removed with an endoscope. Serious and life-threatening situations are caused by the ingestion of foreign bodies in about 1% of all cases.

CASE REPORT

The authors present a case of a two-year-old girl diagnosed with acute abdomen for which she was operated on. A small bowel volvulus and several intestinal fistulas were found intraoperatively. The cause of this finding was the ingestion of magnetic balls and a swallowed metal body drawn to them by magnetic force.

CONCLUSION

If more than one magnetic body is ingested, it is necessary to admit the patient to hospital and to remove these foreign bodies using an endoscope. The position of the magnets which is not changing in a location inaccessible for an endoscope during 2448 hours is an indication for urgent operation.

Redesigning services for the management of vitreomacular traction and macular hole

Vitreomacular traction (VMT) and VMT with macular hole (MH) are serious conditions, being associated with visual disturbance, for example, metamorphopsia, and diminished visual acuity (VA). Pars plana vitrectomy is the routine treatment for symptomatic VMT and VMT+MH. However, ocriplasmin has demonstrated favourable efficacy and safety in specific patient groups with VMT/MH and is now recommended as a treatment option for certain patients by the National Institute of Health and Care Excellence. This means that services for managing patients with VMT/MH may need to be revised, as patients can now potentially receive treatment earlier in the course of the disease. VMT triage clinics could provide a more efficient way of managing VMT/MH patients. Patient assessment should always include high-definition optical coherence tomography, as this is the most accurate means of assessing abnormalities in the vitreoretinal (VR) interface, and an accurate measurement of best-corrected VA. It has been proposed that patients with VMT+MH be managed as a routine 6-week referral, with the complete patient journey—from initial referral to treatment—taking no longer than 6 months. It is important that patients are entered onto VR surgical lists so that there is no delay if ocriplasmin treatment is unsuccessful. Patients will need appropriate counselling about the expected outcomes and possible side effects of ocriplasmin treatment. One-year follow-up data should be collected by treatment centres in order to evaluate the new VMT service.

Saccharomyces cerevisiae as a model for the study of extranuclear functions of mammalian telomerase

The experimental evidence from the last decade made telomerase a prominent member of a family of moonlighting proteins performing different functions at various cellular loci. However, the study of extratelomeric functions of the catalytic subunit of mammalian telomerase (TERT) is often complicated by the fact that it is sometimes difficult to distinguish them from its role(s) at the chromosomal ends. Here, we present an experimental model for studying the extranuclear function(s) of mammalian telomerase in the yeast Saccharomyces cerevisiae. We demonstrate that the catalytic subunit of mammalian telomerase protects the yeast cells against oxidative stress and affects the stability of the mitochondrial genome. The advantage of using S. cerevisiae to study of mammalian telomerase is that (1) mammalian TERT does not interfere with its yeast counterpart in the maintenance of telomeres, (2) yeast telomerase is not localized in mitochondria and (3) it does not seem to be involved in the protection of cells against oxidative stress and stabilization of mtDNA. Thus, yeast cells can be used as a 'test tube' for reconstitution of mammalian TERT extranuclear function(s).

Replication intermediates of the linear mitochondrial DNA of Candida parapsilosis suggest a common recombination based mechanism for yeast mitochondria

Variation in the topology of mitochondrial DNA (mtDNA) in eukaryotes evokes the question if differently structured DNAs are replicated by a common mechanism. RNA-primed DNA synthesis has been established as a mechanism for replicating the circular animal/mammalian mtDNA. In yeasts, circular mtDNA molecules were assumed to be templates for rolling circle DNA-replication. We recently showed that in Candida albicans, which has circular mapping mtDNA, recombination driven replication is a major mechanism for replicating a complex branched mtDNA network. Careful analyses of C. albicans-mtDNA did not reveal detectable amounts of circular DNA molecules. In the present study we addressed the question of how the unit sized linear mtDNA of Candida parapsilosis terminating at both ends with arrays of tandem repeats (mitochondrial telomeres) is replicated. Originally, we expected to find replication intermediates diagnostic of canonical bi-directional replication initiation at the centrally located bi-directional promoter region. However, we found that the linear mtDNA of Candida parapsilosis also employs recombination for replication initiation. The most striking findings were that the mitochondrial telomeres appear to be hot spots for recombination driven replication, and that stable RNA:DNA hybrids, with a potential role in mtDNA replication, are also present in the mtDNA preparations.

Mathematical model of alternative mechanism of telomere length maintenance

Biopolymer length regulation is a complex process that involves a large number of biological, chemical, and physical subprocesses acting simultaneously across multiple spatial and temporal scales. An illustrative example important for genomic stability is the length regulation of telomeres-nucleoprotein structures at the ends of linear chromosomes consisting of tandemly repeated DNA sequences and a specialized set of proteins. Maintenance of telomeres is often facilitated by the enzyme telomerase but, particularly in telomerase-free systems, the maintenance of chromosomal termini depends on alternative lengthening of telomeres (ALT) mechanisms mediated by recombination. Various linear and circular DNA structures were identified to participate in ALT, however, dynamics of the whole process is still poorly understood. We propose a chemical kinetics model of ALT with kinetic rates systematically derived from the biophysics of DNA diffusion and looping. The reaction system is reduced to a coagulation-fragmentation system by quasi-steady-state approximation. The detailed treatment of kinetic rates yields explicit formulas for expected size distributions of telomeres that demonstrate the key role played by the J factor, a quantitative measure of bending of polymers. The results are in agreement with experimental data and point out interesting phenomena: an appearance of very long telomeric circles if the total telomere density exceeds a critical value (excess mass) and a nonlinear response of the telomere size distributions to the amount of telomeric DNA in the system. The results can be of general importance for understanding dynamics of telomeres in telomerase-independent systems as this mode of telomere maintenance is similar to the situation in tumor cells lacking telomerase activity. Furthermore, due to its universality, the model may also serve as a prototype of an interaction between linear and circular DNA structures in various settings.

[Chronic abdominal pain with a rare cause]

The authors use a case report of a female patient with chronic abdominal pain progressing postprandially as an example to demonstrate the complexity of the diagnostic process in patients with these complaints, as well as a surprising peroperative finding in which a small neuroendocrine tumour of the small intestine was the cause of the difficulties. The ensuing discussion also deals with other possible causes of chronic abdominal pain associated with eating disorders.

Mitochondrial genome of the basidiomycetous yeast Jaminaea angkorensis

Jaminaea angkorensis is an anamorphic basidiomycetous yeast species originally isolated from decaying leaves in Cambodia. Taxonomically, J. angkorensis is affiliated with Microstromatales (Exobasidiomycetes, Ustilaginomycotina, Basidiomycota) and represents a basal phylogenetic lineage of this fungal order. To perform a comparative analysis of J. angkorensis with other basidiomycetes, we determined and analyzed its complete mitochondrial DNA sequence. The mitochondrial genome is represented by 29,999 base pairs long, circular DNA containing 32 % guanine and cytosine residues. Its genetic organization is relatively compact and comprises typical genes for 15 conserved proteins involved in oxidative phosphorylation (atp6, 8, and 9; cob; cox1, 2, and 3; and nad1, 2, 3, 4, 4L, 5, and 6) and translation (rps3), two ribosomal RNAs (rnl and rns) and twenty-two transfer RNAs (trnA-Y). Although the gene content is similar to other basidiomycetes, the gene orders in the examined species exhibit only a limited synteny, reflecting their phylogenetic distances and extensive genome rearrangements. In addition, a comparative analysis of basidiomycete mitochondrial genomes indicates that stop-to-tryptophan reassignment of the UGA codon was accompanied by structural alterations of tRNA-Trp(CCA). These results provide an insight into the evolution of the genetic code in fungal mitochondria.

Synergism of the two Myb domains of Tay1 protein results in high affinity binding to telomeres

Double-stranded regions of the telomeres are recognized by proteins containing Myb-like domains conferring specificity toward telomeric repeats. Although biochemical and structural studies revealed basic molecular principles involved in DNA binding, relatively little is known about evolutionary pathways leading to various types of Myb domain-containing proteins in divergent species of eukaryotes. Recently we identified a novel type of telomere-binding protein YlTay1p from the yeast Yarrowia lipolytica containing two Myb domains (Myb1, Myb2) very similar to the Myb domain of mammalian TRF1 and TRF2. In this study we prepared mutant versions of YlTay1p lacking Myb1, Myb2, or both Myb domains and found that YlTay1p carrying either Myb domain exhibits preferential affinity to both Y. lipolytica (GGGTTAGTCA)(n) and human (TTAGGG)(n) telomeric sequences. Quantitative measurements of the protein binding to telomeric DNA revealed that the presence of both Myb domains is required for a high affinity of YlTay1p to either telomeric repeat. Additionally, we performed detailed thermodynamic analysis of the YlTay1p interaction with its cognate telomeric DNA, which is to our knowledge the first energetic description of a full-length telomeric-protein binding to DNA. Interestingly, when compared with human TRF1 and TRF2 proteins, YlTay1p exhibited higher affinity not only for Y. lipolytica telomeres but also for human telomeric sequences. The duplication of the Myb domain region in YlTay1p thus produces a synergistic effect on its affinity toward the cognate telomeric sequence, alleviating the need for homodimerization observed in TRF-like proteins possessing a single Myb domain.

Gentisate and 3-oxoadipate pathways in the yeast Candida parapsilosis: identification and functional analysis of the genes coding for 3-hydroxybenzoate 6-hydroxylase and 4-hydroxybenzoate 1-hydroxylase

The pathogenic yeast Candida parapsilosis degrades various hydroxy derivatives of benzenes and benzoates by the gentisate and 3-oxoadipate pathways. We identified the genes MNX1, MNX2, MNX3, GDX1, HDX1 and FPH1 that code for enzymes involved in these pathways in the complete genome sequence of C. parapsilosis. Next, we demonstrated that MNX1, MNX2, MNX3 and GDX1 are inducible and transcriptionally controlled by hydroxyaromatic substrates present in cultivation media. Our results indicate that MNX1 and MNX2 code for flavoprotein monooxygenases catalysing the first steps in the 3-oxoadipate and gentisate pathways, respectively (i.e. 4-hydroxybenzoate 1-hydroxylase and 3-hydroxybenzoate 6-hydroxylase). Moreover, we found that the two pathways differ by their intracellular localization. The enzymes of the 3-oxoadipate pathway, Mnx1p and Mnx3p, localize predominantly in the cytosol. In contrast, intracellular localization of the components of the gentisate pathway, Mnx2p and Gdx1p, depends on the substrate in the cultivation medium. In cells growing on glucose these proteins localize in the cytosol, whereas in media containing hydroxyaromatic compounds they associate with mitochondria. Finally, we showed that the overexpression of MNX1 or MNX2 increases the tolerance of C. parapsilosis cells to the antifungal drug terbinafine.

Evolution of linear chromosomes and multipartite genomes in yeast mitochondria

Mitochondrial genome diversity in closely related species provides an excellent platform for investigation of chromosome architecture and its evolution by means of comparative genomics. In this study, we determined the complete mitochondrial DNA sequences of eight Candida species and analyzed their molecular architectures. Our survey revealed a puzzling variability of genome architecture, including circular- and linear-mapping and multipartite linear forms. We propose that the arrangement of large inverted repeats identified in these genomes plays a crucial role in alterations of their molecular architectures. In specific arrangements, the inverted repeats appear to function as resolution elements, allowing genome conversion among different topologies, eventually leading to genome fragmentation into multiple linear DNA molecules. We suggest that molecular transactions generating linear mitochondrial DNA molecules with defined telomeric structures may parallel the evolutionary emergence of linear chromosomes and multipartite genomes in general and may provide clues for the origin of telomeres and pathways implicated in their maintenance.