Novel Derivatives of Quinoxaline-2-carboxylic Acid 1,4-Dioxides as Antimycobacterial Agents: Mechanistic Studies and Therapeutic Potential

The World Health Organization (WHO) reports that tuberculosis (TB) is one of the top 10 leading causes of global mortality. The increasing incidence of multidrug-resistant TB highlights the urgent need for an intensified quest to discover innovative anti-TB medications In this study, we investigated four new derivatives from the quinoxaline-2-carboxylic acid 1,4-dioxide class. New 3-methylquinoxaline 1,4-dioxides with a variation in substituents at positions 2 and 6(7) were synthesized via nucleophilic aromatic substitution with amines and assessed against a Mycobacteria spp. Compound 4 showed high antimycobacterial activity (1.25 μg/mL against M. tuberculosis) and low toxicity in vivo in mice. Selection and whole-genomic sequencing of spontaneous drug-resistant M. smegmatis mutants revealed a high number of single-nucleotide polymorphisms, confirming the predicted mode of action of the quinoxaline-2-carboxylic acid 1,4-dioxide 4 as a DNA-damaging agent. Subsequent reverse genetics methods confirmed that mutations in the genes MSMEG_4646, MSMEG_5122, and MSMEG_1380 mediate resistance to these compounds. Overall, the derivatives of quinoxaline-2-carboxylic acid 1,4-dioxide present a promising scaffold for the development of innovative antimycobacterial drugs.

Kanamycin and Ofloxacin Activate the Intrinsic Resistance to Multiple Antibiotics in Mycobacterium smegmatis

Drug resistance (DR) in Mycobacterium tuberculosis is the main problem in fighting tuberculosis (TB). This pathogenic bacterium has several types of DR implementation: acquired and intrinsic DR. Recent studies have shown that exposure to various antibiotics activates multiple genes, including genes responsible for intrinsic DR. To date, there is evidence of the acquisition of resistance at concentrations well below the standard MICs. In this study, we aimed to investigate the mechanism of intrinsic drug cross-resistance induction by subinhibitory concentrations of antibiotics. We showed that pretreatment of M. smegmatis with low doses of antibiotics (kanamycin and ofloxacin) induced drug resistance. This effect may be caused by a change in the expression of transcriptional regulators of the mycobacterial resistome, in particular the main transcriptional regulator whiB7.

Characterization of Genetic Variants Associated with Rifampicin Resistance Level in Mycobacterium tuberculosis Clinical Isolates Collected in Guangzhou Chest Hospital, China

Objective

Rifampicin (RIF)-resistance, a surrogate marker for multidrug-resistant tuberculosis (TB), is mediated by mutations in the rpoB gene. We aimed to investigate the prevalence of mutations pattern in the entire rpoB gene of Mycobacterium tuberculosis clinical isolates and their association with resistance level to RIF.

Methods

Among 465 clinical isolates collected from the Guangzhou Chest Hospital, drug-susceptibility of 175 confirmed Mtb strains was performed via the proportion method and Bactec MGIT 960 system. GeneXpert MTB/RIF and sanger sequencing facilitated in genetic characterization, whereas the MICs of RIF were determined by Alamar blue assay.

Results

We found 150/175 (85.71%) RIF-resistant strains (MIC: 4 to >64 µg/mL) of which 57 were MDR and 81 pre-XDR TB. Genetic analysis identified 17 types of mutations 146/150 (97.33%) within RRDR (codons 426–452) of rpoB, mainly at L430 (P), D435 (V, E, G, N), H445 (N, D, Y, R, L), S450 (L, F) and L452 (P). D435V 12/146 (8.2%), H445N 16/146 (10.9%), and S450L 70/146 (47.94%) were the most frequently encountered mutations. Mutations Q432K, M434V, and N437D are rarely identified in RRDR. Deletions at (1284–1289 CCAGCT), (1295–1303 AATTCATGG), and insertion at (1300–1302 TTC) were detected within RRDR of three RIF^R strains for the first time. We detected 47 types of mutations and insertions/deletions (indels) outside the RRDR. Four RIF^R strains were detected with only novel mutations/indels outside the RRDR. Two of the four had (K274Q + C897 del + I491M) and (A286V + L494P), respectively. The other two had (G1687del + P454L) and (TT1835-6 ins + I491L) individually. Compared with phenotypic characterization, diagnostic sensitivities of GeneXpert MTB/RIF and sequencing analysis were 95.33% (143/150), and 100% (150/150) respectively.

Conclusion

Our findings underscore the key role of RRDR mutations and the contribution of non-RRDR mutations in rapid molecular diagnosis of RIF^R clinical isolates. Such insights will support early detection of disease and recommend the appropriate anti-TB regimens in high-burden settings.

Synthesis and Characterization of Novel 2-Acyl-3-trifluoromethylquinoxaline 1,4-Dioxides as Potential Antimicrobial Agents

The emergence of drug resistance in pathogens leads to a loss of effectiveness of antimicrobials and complicates the treatment of bacterial infections. Quinoxaline 1,4-dioxides represent a prospective scaffold for search of new compounds with improved chemotherapeutic characteristics. Novel 2-acyl-3-trifluoromethylquinoxaline 1,4-dioxides with alteration of substituents at position 2 and 6 were synthesized via nucleophilic substitution with piperazine moiety and evaluated against a broad panel of bacteria and fungi by measuring their minimal inhibitory concentrations. Their mode of action was assessed by whole-genomic sequencing of spontaneous drug-resistant Mycobacterium smegmatis mutants, followed by comparative genomic analysis, and on an original pDualrep2 system. Most of the 2-acyl-3-trifluoromethylquinoxaline 1,4-dioxides showed high antibacterial properties against Gram-positive strains, including mycobacteria, and the introduction of a halogen atom in the position 6 of the quinoxaline ring further increased their activity, with 13c being the most active compound. The mode of action studies confirmed the DNA-damaging nature of the obtained quinoxaline 1,4-dioxides, while drug-resistance may be provided by mutations in redox homeostasis genes, encoding enzymes potentially involved in the activation of the compounds. This study extends views about the antimicrobial and antifungal activities of the quinoxaline 1,4-dioxides and can potentially lead to the discovery of new antibacterial drugs.

Transcriptomic Profile of Mycobacterium smegmatis in Response to an Imidazo[1,2-b][1,2,4,5]tetrazine Reveals Its Possible Impact on Iron Metabolism

Tuberculosis (TB), caused by the Mycobacterium tuberculosis complex bacteria, is one of the most pressing health problems. The development of new drugs and new therapeutic regimens effective against the pathogen is one of the greatest challenges in the way of tuberculosis control. Imidazo[1,2-b][1,2,4,5]tetrazines have shown promising activity against M. tuberculosis and M. smegmatis strains. Mutations in MSMEG_1380 lead to mmpS5–mmpL5 operon overexpression, which provides M. smegmatis with efflux-mediated resistance to imidazo[1,2-b][1,2,4,5]tetrazines, but the exact mechanism of action of these compounds remains unknown. To assess the mode of action of imidazo[1,2-b][1,2,4,5]tetrazines, we analyzed the transcriptomic response of M. smegmatis to three different concentrations of 3a compound: 1/8×, 1/4×, and 1/2× MIC. Six groups of genes responsible for siderophore synthesis and transport were upregulated in a dose-dependent manner, while virtual docking revealed proteins involved in siderophore synthesis as possible targets for 3a. Thus, we suggest that imidazo[1,2-b][1,2,4,5]tetrazines may affect mycobacterial iron metabolism.

Transcriptomic dataset of Mycolicibacterium smegmatis exposed to an imidazo[1,2-b][1,2,4,5]tetrazine

Deciphering the mechanism of action of novel anti-tuberculosis compounds is a key step in the drug development process. We have previously described a number of imidazo[1,2-b][1,2,4,5]tetrazines with a promising activity on Mycobacterium tuberculosis[1]. These compounds had predicted activity as serine‑threonine protein kinase inhibitors, however spontaneous drug resistant Mycolicibacterium smegmatis mc ² 155 (formerly Mycobacterium smegmatis) revealed only the mycobacterial mechanism of resistance to imidazo[1,2-b][1,2,4,5]tetrazines: mutations in MSMEG_1380 gene lead to overexpression of the mmpS5-mmpL5 operon in M. smegmatis, thus providing resistance to imidazo[1,2-b][1,2,4,5]tetrazines via enhanced efflux [2]. Here we report the RNA sequencing data of M. smegmatis mc ²  155 culture treated with one of the imidazo[1,2-b][1,2,4,5]tetrazines for 1.5 h and the untreated culture as a control. The mapped reads showed that a total of 1386 genes are differentially expressed in this experiment. A further analysis of these data can shed light of the mechanism of action of imidazo[1,2-b][1,2,4,5]tetrazines. The data generated by RNA-seq (raw reads) have been deposited to NCBI sequence read archive (SRA) and have been assigned a BioProject accession number PRJNA615922.

Transcriptional regulation of drug resistance mechanisms in Salmonella: where we stand and what we need to know

Salmonellae have evolved a wide range of molecular mechanisms to neutralize the effect of antibiotics and evade the host immune system response. These mechanisms are exquisitely controlled by global and local regulators and enable the pathogens to use its energy as per need and hence allow the pathogen to economize the consumption of energy by its cellular machinery. Several families that regulate the expression of different drug resistance genes are known; some of these are: the TetR family (which affects tetracycline resistance genes), the AraC/XylS family (regulators that can act as both transcriptional activators and repressors), two-component signal transduction systems (e.g. PhoPQ, a key regulator for virulence), mercury resistance Mer-R and multiple antibiotic resistance Mar-R regulators, LysR-type global regulators (e.g. LeuO) and histone-like protein regulators (involved in the repression of newly transferred resistance genes). This minireview focuses on the role of different regulators harbored by the Salmonella genome and characterized for mediating the drug resistance mechanisms particularly via efflux and influx systems. Understanding of such transcriptional regulation mechanisms is imperative to address drug resistance issues in Salmonella and other bacterial pathogens.

MmpS5-MmpL5 Transporters Provide Mycobacterium smegmatis Resistance to imidazo[1,2-b][1,2,4,5]tetrazines

The emergence and spread of drug-resistant Mycobacterium tuberculosis strains (including MDR, XDR, and TDR) force scientists worldwide to search for new anti-tuberculosis drugs. We have previously reported a number of imidazo[1,2-b][1,2,4,5]tetrazines - putative inhibitors of mycobacterial eukaryotic-type serine-threonine protein-kinases, active against M. tuberculosis. Whole genomic sequences of spontaneous drug-resistant M. smegmatis mutants revealed four genes possibly involved in imidazo[1,2-b][1,2,4,5]tetrazines resistance; however, the exact mechanism of resistance remain unknown. We used different approaches (construction of targeted mutants, overexpression of the wild-type (w.t.) and mutant genes, and gene-expression studies) to assess the role of the previously identified mutations. We show that mutations in MSMEG_1380 gene lead to overexpression of the mmpS5-mmpL5 operon in M. smegmatis, thus providing resistance to imidazo[1,2-b][1,2,4,5]tetrazines by increased efflux through the MmpS5-MmpL5 system, similarly to the mechanisms of resistance described for M. tuberculosis and M. abscessus. Mycobacterial MmpS5-MmpL5 transporters should be considered as an MDR-efflux system and they should be taken into account at early stages of anti-tuberculosis drug development.

Draft Genome Sequences of Mycobacterium tuberculosis Clinical Isolates from the Ural Region of Russia That Carry the pks15/1 Gene

Here, we report the draft genome sequences of 15 Mycobacterium tuberculosis isolates of the Beijing-B0/W-148 sublineage that carry a 7-bp insertion within the pks15 gene, which leads to the synthesis of Pks15/1 fusion protein. Pks15/1 is involved in phenolglycolipid synthesis and biofilm formation, thus potentially contributing to the B0/W-148 lineage's enhanced virulence and drug resistance.

Sequencing and Analysis of Three Mycobacterium tuberculosis Genomes of the B0/N-90 Sublineage

We report the draft genome sequences of three Mycobacterium tuberculosis isolates belonging to the B0/N-90 sublineage, EKB34, EKB53, and EKB79. The B0/N-90 sublineage belongs to the prevalent (in Russia) and highly virulent Beijing-B0/W148 sublineage. Isolates EKB34 and EKB79 were obtained from people with immune deficiency.

Synthesis and antimycobacterial activity of imidazo[1,2-b][1,2,4,5]tetrazines

Tuberculosis (TB) has recently become the leading killer among infectious diseases. Multidrug and extensively drug-resistant Mycobacterium tuberculosis strains urge the need to develop anti-TB drugs with a novel mechanism of action. We describe synthesis of 22 novel imidazo[1,2-b][1,2,4,5]tetrazine derivatives with different substituents at C(3) and C(6) positions, and their antimycobacterial activity in vitro. 8 compounds show activity as potential serine/threonine protein kinase (STPK) inhibitors in M. smegmatis aphVIII+ test-system, which is characteristic for this class. 3 compounds out of 5 most active STPK inhibitors have a prominent minimal inhibitory concentration on M. tuberculosis H37Rv of 1 μg/ml. We were able to obtain M. smegmatis mc2 155 mutants resistant to 4 compounds and show that they do not have cross resistance with other drugs, but have a common mechanism of resistance among these 4 imidazo[1,2-b][1,2,4,5]tetrazines. Compound 3h seems the most promising, combining a predicted STPK inhibitor activity, the lowest MIC on M. tuberculosis and a low frequency of drug resistant mutants' emergence.

Draft Genome Sequences of 12 Mycolicibacterium smegmatis Strains Resistant to Imidazo[1,2-b][1,2,4,5]Tetrazines

Here, we report 12 draft genome sequences of mutant Mycolicibacterium smegmatis strains resistant to imidazo[1,2-b][1,2,4,5]tetrazines, which are antituberculosis drug candidates. We have identified 7 different mutations in the MSMEG_1380 gene, which encodes the AcrR/TetR_N transcriptional repressor, which may activate efflux-mediated resistance.

Mutations in Efflux Pump Rv1258c (Tap) Cause Resistance to Pyrazinamide, Isoniazid, and Streptomycin in M. tuberculosis

Although drug resistance in Mycobacterium tuberculosis is mainly caused by mutations in drug activating enzymes or drug targets, there is increasing interest in the possible role of efflux in causing drug resistance. Previously, efflux genes have been shown to be upregulated upon drug exposure or implicated in drug resistance in overexpression studies, but the role of mutations in efflux pumps identified in clinical isolates in causing drug resistance is unknown. Here we investigated the role of mutations in efflux pump Rv1258c (Tap) from clinical isolates in causing drug resistance in M. tuberculosis. We constructed point mutations V219A and S292L in Rv1258c in the chromosome of M. tuberculosis and the point mutations were confirmed by DNA sequencing. The susceptibility of the constructed M. tuberculosis Rv1258c mutants to different tuberculosis drugs was assessed using conventional drug susceptibility testing in 7H11 agar in the presence and absence of efflux pump inhibitor piperine. A C14-labeled PZA uptake experiment was performed to demonstrate higher efflux activity in the M. tuberculosis Rv1258c mutants. Interestingly, the V219A and S292L point mutations caused clinically relevant drug resistance to pyrazinamide (PZA), isoniazid (INH), and streptomycin (SM), but not to other drugs in M. tuberculosis. While V219A point mutation conferred low-level drug resistance, the S292L mutation caused a higher level of resistance. Efflux inhibitor piperine inhibited INH and PZA resistance in the S292L mutant but not in the V219A mutant. The S292L mutant had higher efflux activity for pyrazinoic acid (the active form of PZA) than the parent strain. We conclude that point mutations in the efflux pump Rv1258c in clinical isolates can confer clinically relevant drug resistance, including PZA resistance, and could explain some previously unaccounted drug resistance in clinical strains. Future studies need to take efflux mutations into consideration for improved detection of drug resistance in M. tuberculosis and address their role in affecting treatment outcome in vivo.

Resistance to pyrazinamide in Russian Mycobacterium tuberculosis isolates: pncA sequencing versus Bactec MGIT 960

Resistance to pyrazinamide (PZA) may impact clinical outcome of anti-tuberculosis chemotherapy. PZA susceptibility testing using MGIT 960 is not reliable and little information is available on the prevalence of PZA resistance in Russia. A collection of 64 clinical isolates of Mycobacterium tuberculosis, including 35 multidrug resistant and extensively drug-resistant (MDR/XDR), was analyzed for PZA resistance using MGIT 960, Wayne test, and sequencing of PZA resistance genes pncA, rpsA and panD. In addition, we analyzed 519 MDR-TB strains for susceptibility to PZA by MGIT 960. Sequencing of pncA revealed 17 of 25 (68%) MDR strains and all 10 XDR strains harboring pncA mutations. A correlation of φ = 0.81 between MGIT 960 and pncA sequencing was observed. Mutations in rpsA and panD not associated with PZA resistance as defined by MGIT 960 were identified. We found 1 PZA-resistant strain without mutations in known PZA resistance genes. Almost 73% of MDR-TB strains isolated in Moscow, Russia, were PZA-resistant by MGIT 960 testing of 519 MDR-TB clinical isolates. Further studies are needed to determine the role of rpsA and panD mutations in possible low-level PZA resistance and to identify the molecular basis of new PZA resistance in the isolate without known PZA resistance mutations.

Leptomonas costaricensis sp. n. (Kinetoplastea: Trypanosomatidae), a member of the novel phylogenetic group of insect trypanosomatids closely related to the genus Leishmania

A flagellate isolated from the intestinal tract of a reduviid bug Ricolla simillima (Heteroptera) in Costa Rica was found to represent a new trypanosomatid species by the phylogenetic analysis of small subunit ribosomal RNA (SSU rRNA), glyceraldehyde phosphate dehydrogenase (GAPDH) and large subunit of RNA polymerase II (RPOIILS) genes. The phylogenetic position of this trypanosomatid, together with its typical promastigote morphology and the host identity, allowed its classification as a species that belongs to the polyphyletic genus Leptomonas. Interestingly, the new species was revealed as a member of the novel phylogenetic clade representing the closest known relative of Leishmania. With the new species used as an outgroup to root the Leishmania RPOIILS phylogenetic tree, the lineage of the Neotropical species L. enriettii was found to branch off early, and was followed by a deep split between the Old World and the remaining New World species. This tree topology supports the hypothesis that the initial transition to dixenous parasitism in this group pre-dated the continental split and that afterwards the Neotropical and the Old World groups evolved largely independently.

Trypanosomatid biodiversity in Costa Rica: genotyping of parasites from Heteroptera using the spliced leader RNA gene

The biodiversity of insect trypanosomes is largely unknown, resulting in significant gaps in the understanding of pathogen evolution. A culture-independent preliminary survey of trypanosomatid fauna was conducted for the parasites of Heteroptera (Hemiptera) from several localities in Costa Rica. Trypanosomatid infections were detected by light microscopy of smeared gut contents. Out of 257 insects representing 6 families, infections were found in 62 cases; cultures were obtained for 29 new isolates. Gut material from infected hosts was preserved in the field using an SDS–EDTA buffer solution for subsequent DNA extraction in the laboratory. PCR amplification of the trypanosomatid-specific spliced leader (SL) RNA gene repeats was successful for 60 field samples. Eighteen distinct SL RNA typing units were identified in a set of 28 samples analysed in detail. Cluster analysis indicated that these typing units were unique and thus could represent new species and, in some cases, new genera. This study reveals only a minor fraction of the trypanosomatid biodiversity, which is anticipated to be high.

Diplonema spp. possess spliced leader RNA genes similar to the Kinetoplastida

The phylogenetic placement of the genus Diplonema in relation to fellow phylum members Euglena and Trypanosoma has been uncertain. The spliced leader RNA gene, present in the euglenids and kinetoplastids in distinct forms, was a potential target for resolving this question. The first indication supporting a closer relationship to the kinetoplastids was the recognition of potential spliced leader RNA exon sequences in the genomic DNA of two Diplonema isolates. Examination of total cell RNA revealed transcripts in the anticipated size range at approximately 120 and 130 nt. Specific PCR amplification of a spliced leader RNA gene repeat was performed. The hallmark features of the kinetoplastid-type spliced leader RNA, specifically the 39-nt exon, splice-donor site, Sm-binding site and poly-T tract and the potential to form the requisite stem-loop structures, were found. Diplonema spp. are different from the kinetoplastids by virtue of C residues at positions 4 and 18 in the exon. While the intergenic spacer regions varied in size, each contained the complete sequence or remnants of a 5S ribosomal RNA gene. Possession of a functional spliced leader RNA gene of the kinetoplastid variety in Diplonema supports a closer evolutionary relationship with this group than with the euglenids.

Phylogeny of the kinetoplastida: taxonomic problems and insights into the evolution of parasitism

To further investigate phylogeny of kinetoplastid protozoa, the sequences of small subunit (18S) ribosomal RNA of nine bodonid isolates and ten isolates of insect trypanosomatids have been determined. The root of the kinetoplastid tree was attached to the branch of Bodo designis and/or Cruzella marina. The suborder Trypanosomatina appeared as a monophyletic group, while the suborder Bodonina was paraphyletic. Among bodonid lineages, parasitic organisms were intermingled with free-living ones, implying multiple transitions to parasitism and supporting the 'vertebrate-first hypothesis'. The tree indicated that the genera Cryptobia and Bodo are artificial taxa. Separation of fish cryptobias and Trypanoplasma borreli as different genera was not supported. In trypanosomatids, the genera Leptomonas and Blastocrithidia were polyphyletic, similar to the genera Herpetomonas and Crithidia and in contrast to the monophyletic genera Trypanosoma and Phytomonas. This analysis has shown that the morphological classification of kinetoplastids does not in general reflect their genetic affinities and needs a revision.

Diversity and phylogeny of insect trypanosomatids based on small subunit rRNA genes: polyphyly of Leptomonas and Blastocrithidia

With the aim of further investigating phylogenetic relationships in insect trypanosomatids, we have determined the sequences of small subunit rRNA genes from ten isolates, which were originally classified as Leptomonas, Blastocrithidia, and Wallaceina based on their morphology in the hosts. The inferred maximum likelihood, parsimony, and distance trees indicate that the Leptomonas and Blastocrithidia are polyphyletic, and confirm the polyphyly of Herpetomonas and Crithidia. Blastocrithidia triatoma and Leptomonas collosoma were among the earliest branching lineages among the insect trypanosomatids, while most other isolates were found within a closely related terminal clade, which also included Crithidia fasciculata. This analysis has clearly demonstrated that the morphological classification system of insect trypanosomatids does not always reflect their genetic affinities warranting its revision in the future.

Leishmania tarentolae: a parallel isolation of cytochrome bc(1) and cytochrome c oxidase

A rapid and simple method which allowed for a parallel isolation of cytochrome c reductase (cytochrome bc(1) ) and cytochrome c oxidase from kinetoplast-mitochondria of Leishmania tarentolae was developed. The method involved the lysis of kinetoplasts with dodecyl maltoside in the presence of 260 mM NaCl, followed by purification of bc(1) complexes on DEAE-sepharose CL-6B. The oxidase which was found in the flow-through fractions of the first chromatographic step was diluted and then repurified on a similar DEAE-sepharose column. The investigated properties of the isolated cytochrome c oxidase and reductase, such as their absolute and difference spectrum absorption maxima, heme content, specific activity, and subunit composition, confirm the usefulness of this method for obtaining highly active preparations of the enzymes.

Phylogeny of the bodonid flagellates (Kinetoplastida) based on small-subunit rRNA gene sequences

The phylogeny of kinetoplastid flagellates was investigated by determining the sequences of the small-subunit (18S) rRNA from Bodo designis, Bodo saltans K, Bodo saltans P, Bodo sorokini, Bodo sp. (cf. uncinatus), Cruzella marina, Cryptobia helicis, Dimastigella mimosa and Parabodo nitrophilus and analysing these data together with several previously obtained sequences. The root of the kinetoplastid tree was tentatively determined to be attached to the branch of B. designis and/or Cruzella marina. Within this topology, the suborder Trypanosomatina appears as a late-emerging monophyletic group, while the suborder Bodonina is paraphyletic. Within the bodonid subtree, the branches of parasitic organisms were intermingled with free-living ones, implying multiple transitions to parasitism. The tree indicates that the genera Cryptobia and Bodo are artificial taxa. In addition, the separation of the fish cryptobias and Trypanoplasma borreli as different genera was not supported.

The absence of genes for cytochrome c oxidase and reductase subunits in maxicircle kinetoplast DNA of the respiration-deficient plant trypanosomatid Phytomonas serpens

By completing the sequencing of the maxicircle conserved region in the kinetoplast DNA of Phytomonas serpens, we showed that the genes for subunits I and II (COI and COII) of cytochrome c oxidase in this organism were missing. We had previously shown that the genes for cytochrome c oxidase subunit III and apocytochrome b were also missing. These deletions did not affect the structure or expression of the remaining genes. Partial editing of the mRNA for NADH dehydrogenase subunit 8, previously found in strain IG from insects, was complete in two other strains isolated from plants. The appearance of a novel maxicircle gene for MURF2 block I gRNA, which substitutes for the gene missing due to the COII gene deletion, may illustrate a general mechanism for the origin of gRNAs.

Evolution of RNA editing in trypanosome mitochondria

Two different RNA editing systems have been described in the kinetoplast-mitochondrion of trypanosomatid protists. The first involves the precise insertion and deletion of U residues mostly within the coding regions of maxicircle-encoded mRNAs to produce open reading frames. This editing is mediated by short overlapping complementary guide RNAs encoded in both the maxicircle and the minicircle molecules and involves a series of enzymatic cleavage-ligation steps. The second editing system is a C(34) to U(34) modification in the anticodon of the imported tRNA(Trp), thereby permitting the decoding of the UGA stop codon as tryptophan. U-insertion editing probably originated in an ancestor of the kinetoplastid lineage and appears to have evolved in some cases by the replacement of the original pan-edited cryptogene with a partially edited cDNA. The driving force for the evolutionary fixation of these retroposition events was postulated to be the stochastic loss of entire minicircle sequence classes and their encoded guide RNAs upon segregation of the single kinetoplast DNA network into daughter cells at cell division. A large plasticity in the relative abundance of minicircle sequence classes has been observed during cell culture in the laboratory. Computer simulations provide theoretical evidence for this plasticity if a random distribution and segregation model of minicircles is assumed. The possible evolutionary relationship of the C to U and U-insertion editing systems is discussed.

Detection of the mitochondrially encoded cytochrome c oxidase subunit I in the trypanosomatid protozoan Leishmania tarentolae. Evidence for translation of unedited mRNA in the kinetoplast

With the aim of identification of kinetoplast-encoded proteins we investigated the subunit composition of cytochrome c oxidase (respiratory complex IV) from kinetoplast mitochondria of the trypanosomatid protozoan Leishmania tarentolae. Eleven stoichiometric subunits were visible in Coomassie-stained, two-dimensional Blue Native/Tricine-SDS electrophoretic gels. Their partial amino acid sequences indicated that these polypeptides are nuclear-encoded. The mitochondrial subunit I was detected with the polyclonal antibodies against an internal region of this polypeptide. In two-dimensional (9 versus 14%) polyacrylamide glycine-SDS gels this subunit is found as a series of spots located off the main diagonal, a property that can be explained by abnormal electrophoretic migration and aggregation. In gels loaded with high amounts of the purified, enzymatically active oxidase, the subunit I spots could be visualized by staining. The determined N-terminal amino acid sequence of the putative monomeric subunit I (MFXLCLVCLSVS) matched with the predicted sequence, thus indicating that the corresponding kinetoplast unedited mRNA is translated into a functional protein.

Translation of the edited mRNA for cytochrome b in trypanosome mitochondria

The type of RNA editing found in the kinetoplast-mitochondria of trypanosomes and related protozoa, involving uridylate insertions and deletions, creates translatable messenger RNAs (mRNAs) out of nonsense pre-edited RNAs by correcting encoded defects that vary from simple frameshifts to large "cryptic" regions. However, any evidence for translation of these mRNAs in the kinetoplast has been missing for decades. We identified a kinetoplast-encoded protein, apocytochrome b, whose mRNA is edited in the 5' region. The determined amino-terminal sequence of the protein coincides with the predicted sequence derived from the edited region, demonstrating that the cognate apocytochrome b mRNA is translated into a functional protein. This finding represents the first direct evidence for a functional translation system in the kinetoplasts.

Unexpectedly high variability of the histone H4 gene in Leishmania

The sequence of the cDNA of the histone H4 gene of Leishmania tarentolae is reported herein. The predicted 100-amino-acid-long protein has the highest degree of identity with the histone H4 gene of L. infantum and shares with it a 5' region that shows a very low degree of identity with the corresponding region of histone H4 genes from other organisms. However, between these two genes is a 7.7% nucleotide difference that results in seven different amino acids, located in the 5', central, and 3' regions of the coding sequence. Such a divergence in the H4 gene, which is considered to be one of the most highly conserved genes, between closely related members of the genus Leishmania is unexpected and may reflect some unusual features of these important proteins in kinetoplastid flagellates.

Partial kinetoplast-mitochondrial gene organization and expression in the respiratory deficient plant trypanosomatid Phytomonas serpens

In plant-dwelling trypanosomatids from the genus Phytomonas, mitochondrial functions, such as cytochrome mediated respiration, ATP production and Krebs cycle, are missing, and cell energetics is based on the glycolysis. Using Blue Native/Tricine-SDS two-dimensional gel electrophoretic analysis, we observed that mitochondrial respiratory Complexes III (cytochrome bc1) and IV (cytochrome c oxidase) were absent in Phytomonas serpens; however, Complex V (ATPase) was present. A deletion of the genes for cytochrome c oxidase subunit III (COIII) and apocytochrome b (Cyb) was identified within the 6234 bp sequenced region of the 31 kb maxicircle kinetoplast DNA. Genes, found in this region, include 12S and 9S ribosomal RNAs, subunits 7, 8 and 9 of NADH dehydrogenase (ND7, ND8 and ND9) and subunit 6 of ATPase (A6 or MURF4), as well as the genes (MURF1, MURF5 and G3) with unknown function. Most genes are actively transcribed and some mRNAs are edited. Fully edited mRNAs for A6 and G3 were abundant, while edited ND7 transcripts were rare, and only partially edited and pre-edited transcripts for ND8 were detected. The data show that the mitochondrial genome of P. serpens is functional, although its functions may be limited to expressing the ATPase and, possibly, NADH dehydrogenase complexes.

Phylogenetic affinities of Diplonema within the Euglenozoa as inferred from the SSU rRNA gene and partial COI protein sequences

In order to shed light on the phylogenetic position of diplonemids within the phylum Euglenozoa, we have sequenced small subunit rRNA (SSU rRNA) genes from Diplonema (syn. Isonema) papillatum and Diplonema sp. We have also analyzed a partial sequence of the mitochondrial gene for cytochrome c oxidase subunit I from D. papillatum. With both markers, the maximum likelihood method favored a closer grouping of diplonemids with kinetoplastids, while the parsimony and distance suggested a closer relationship of diplonemids with euglenoids. In each case, the differences between the best tree and the alternative trees were small. The frequency of codon usage in the partial D. papillatum COI was different from both related groups; however, as is the case in kinetoplastids but not in Euglena, both the non-canonical UGA codon and the canonical UGG codon were used to encode tryptophan in Diplonema.

Demonstration of mRNA editing and localization of guide RNA genes in kinetoplast-mitochondria of the plant trypanosomatid Phytomonas serpens

Maxicircle molecules of kDNA in several isolates of Phytomonas were detected by hybridization with the 12S rRNA gene probe from Leishmania tarentolae. The estimated size of maxicircles is isolate-specific and varies from 27 to 36 kb. Fully edited and polyadenylated mRNA for kinetoplast-encoded ribosomal protein S12 (RPS12) was found in the steady-state kinetoplast RNA isolated from Phytomonas serpens strain 1G. Two minicircles (1.45 kb) from this strain were also sequenced. Each minicircle contains two 120 bp conserved regions positioned 180 degrees apart, a region enriched with G and T bases and a variable region. One minicircle encodes a gRNA for the first block of editing of RPSl2 mRNA, and the other encodes a gRNA with unknown function. A gRNA gene for the second block of RPSl2 was found on a minicircle sequenced previously. On each minicircle, a gRNA gene is located in the variable region in a similar position and orientation with respect to the conserved regions.

Monophyly of endosymbiont containing trypanosomatids: phylogeny versus taxonomy

To obtain additional information on the phylogenetic relationships within the family Trypanosomatidae (order Kinetoplastida), we have sequenced the small subunit ribosomal RNA genes from the endosymbiont containing species Herpetomonas roitmani TCC080, Herpetomonas sp. TCC263, Crithidia oncopelti ATCC 12982 and a partial large subunit rRNA gene from H. roitmani. The small subunit sequences in the two isolates of Herpetomonas are very similar but not identical, and so are their restriction digest profiles of kinetoplast DNA. The size of minicircles in both isolates is 4.2 kilobases. The inferred ribosomal RNA phylogenetic trees shows the genera Herpetomonas and Crithidia as polyphyletic. Endosymbiont-bearing herpetomonads cluster with the endosymbiont-bearing crithidias and a blastocrithidia to form a monophyletic clade, whereas the endosymbiont-free members of these genera are found elsewhere in the tree. These data support the hypothesis of a monophyletic origin of endosymbiosis in trypanosomatid evolution and also suggest that a taxonomic revision is needed in order to better describe the natural affinities in this family.

U-insertion/deletion Edited Sequence Database

Uridine insertion/deletion RNA editing is a post-transcriptional RNA modification occurring in the mitochondria of kinetoplastid protozoa. The U-insertion/deletion Edited Sequence Database is a compilation of mitochondrial genes and edited mRNAs from five kinetoplastid species. It contains separate files with the DNA, mRNA (both unedited and edited) and predicted protein sequences, as well as alignments of the Leishmania tarentolae and Trypanosoma brucei protein sequences from edited and unedited genes. The sequence files are in GCG format. A 'map' sequence file showing the location of U-deletions, U-insertions and the translated amino acid sequences is also provided for each gene. Genomic maps for each species are also provided with clickable genes, including maxicircle-encoded gRNAs. Sets of aligned nuclear rRNA sequences from kinetoplastid protozoa are also provided, which were used for phylogenetic reconstructions in an analysis of the origin of RNA editing. The database is available through the World Wide Web as an HTML document at the URLhttp://www.lifesci.ucla.edu/RNA/trypanosome/ database.html

Analysis of ribosomal RNA genes suggests that trypanosomes are monophyletic

To further investigate the phylogeny of protozoa from the order Kinetoplastida we have sequenced the small subunit (SSU) and a portion of the large subunit (LSU) nuclear rRNA genes. The SSU and LSU sequences were determined from a lizard trypanosome, Trypanosoma scelopori and a bodonid, Rhynchobodo sp., and the LSU sequences were determined from an insect trypanosomatid, Crithidia oncopelti, and a bodonid, Dimastigella trypaniformis. Contrary to previous results, in which trypanosomes were found to be paraphyletic, with Trypanosoma brucei representing the earliest-diverging lineage, we have now found evidence for the monophyly of trypanosomes. Addition of new taxa which subdivide long branches (such as that of T. brucei) have helped to identify homoplasies responsible for the paraphyletic trees in previous studies. Although the monophyly of the trypanosome clade is supported in the bootstrap analyses for maximum likelihood at 97% and maximum parsimony at 92%, there is only a small difference in ln-likelihood value or tree length between the most optimal monophyletic tree and the best suboptimal paraphyletic tree. Within the trypanosomatid subtree, the clade of trypanosomes is a sister group to the monophyletic clade of the nontrypanosome genera. Different groups of trypanosomes group on the tree according to their mode of transmission. This suggests that the adaptation to invertebrate vectors plays a more important role in the trypanosome evolution than the adaptation to vertebrate hosts.

Detection and identification of human pathogenic Leishmania and Trypanosoma species by hybridization of PCR-amplified mini-exon repeats

A single pair of PCR primers within a conserved region of the mini-exon repeat was used to amplify the repeats from 10 species of pathogenic Leishmania belonging to four major clinical groups and also from three species of Trypanosoma. Oligonucleotide hybridization probes for the detection and identification of the PCR-amplified repeats were constructed from alignments of mini-exon intron and intergenic sequences. The probes generated from mini-exon intergenic regions of the L. (V.) braziliensis, L. (L.) donovani, and L. (L.) mexicana species hybridized specifically to their cognate groups without discriminating between the species within the groups. The probes for L. (L.) major and L. (L.) aethiopica were species-specific, while the L. (L.) tropica probe also hybridized with the L. (L.) aethiopica mini-exon repeat. The mini-exon intron-derived probes for T. cruzi, T. rangeli, and T. brucei were species-specific. This method involving the detection of specific PCR-amplified products produced using a single primer set represents a novel sensitive and specific assay for multiple trypanosomatid species and groups.

Phylogeny of trypanosomes as inferred from the small and large subunit rRNAs: implications for the evolution of parasitism in the trypanosomatid protozoa

Sequences of the small rRNA genes and partial sequences of the large rRNA genes were obtained by PCR amplification from a variety of vertebrate trypanosomes. The trypanosome species and hosts included Trypanosoma avium from a bird, T. rotatorium from an amphibian, T. boissoni from an elasmobranch, T. triglae from a marine teleost and T. carassii from a freshwater teleost. Phylogenetic relationships among these species and other representatives of the family Trypanosomatidae were inferred using maximum likelihood, maximum parsimony and evolutionary parsimony. The trypanosomatid tree was rooted using rRNA sequences from two species from the suborder Bodonina. All methods showed that the mammalian parasite, Trypanosoma brucei, constitutes the earliest divergent branch. The remaining trypanosomes formed a monophyletic group. Within this group, the bird trypanosome was grouped with T. cruzi, while the elasmobranch trypanosome and the two fish trypanosome species formed a group with an affinity to T. rotatorium. Our results provide no evidence for co-evolution of trypanosomatids and their hosts, either vertebrate or invertebrate. This suggests that evolution of trypanosomatids was accompanied by secondary acquisitions of hosts and habitats.

Disruption of RNA editing in Leishmania tarentolae by the loss of minicircle-encoded guide RNA genes

RNA editing in kinetoplastids appears to be a labile genetic trait that is affected by prolonged cell culture. The transcripts of the G1-G5 cryptogenes are pan-edited in the recently isolated LEM125 strain of Leishmania tarentolae, but not in the UC strain which has been in culture for 55 years. At least 32 minicircle-encoded guide RNAs (gRNAs) for the editing of G1-G5 transcripts are present in LEM125 and absent in UC. We hypothesize that specific minicircle sequence classes encoding gRNAs for the editing of these transcripts were lost during the long culture history of the UC strain. The protein products, which include components of complex I of the respiratory chain, are probably not required during the culture stage of the Leishmania life cycle.

Ancient origin of RNA editing in kinetoplastid protozoa

A rooted phylogenetic tree of the kinetoplastid protozoa has been constructed that, together with a comparative analysis of editing of several genes, leads to the surprising conclusion that extensive or pan-editing with multiple overlapping guide RNAs is more ancient than 5'-editing. The mechanism of editing is still uncertain, but multiple ribonucleoprotein complexes have been identified that contain components of the enzymatic machinery.

RNA editing and mitochondrial genomic organization in the cryptobiid kinetoplastid protozoan Trypanoplasma borreli

The bodonids and cryptobiids represent an early diverged sister group to the trypanosomatids among the kinetoplastid protozoa. The trypanosome type of uridine insertion-deletion RNA editing was found to occur in the cryptobiid fish parasite Trypanoplasma borreli. A pan-edited ribosomal protein, S12, and a novel 3'- and 5'-edited cytochrome b, in addition to an unedited cytochrome oxidase III gene and an apparently unedited 12S rRNA gene, were found in a 6-kb fragment of the 80- to 90-kb mitochondrial genome. The gene order differs from that in trypanosomatids, as does the organization of putative guide RNA genes; guide RNA-like molecules are transcribed from tandemly repeated 1-kb sequences organized in 200- and 170-kb molecules instead of minicircles. The presence of pan-editing in this lineage is consistent with an ancient evolutionary origin of this process.

Monophyletic origin of beta-division proteobacterial endosymbionts and their coevolution with insect trypanosomatid protozoa Blastocrithidia culicis and Crithidia spp

Some trypanosomatid protozoa (order Kinetoplastida) are well known to harbor bacterial endosymbionts. Their phylogenetic positions and evolutionary relationships with the hosts were deduced by comparing the rRNA gene sequences. Earlier, we observed that these symbionts from three Crithidia spp. are identical and are closely related to Bordetella bronchiseptica. We have now sequenced the genes of another endosymbiont and the host protozoan Blastocrithidia culicis. The 16S rRNA genes of the Blastocrithidia and Crithidia symbionts share approximately 97% identity and form a distinct group, branching off the B. bronchiseptica lineage in the beta-division of Proteobacteria. Comparison of their secondary structures in the stem regions suggests compensatory mutations of the symbiont sequences, contributing to their biased base transitions from G to A and C to T. Two putative genes encoding tRNA(Ile) and tRNA(Ala) are highly conserved in the otherwise variable internal transcribed spacer region. Comparisons of the host rRNA gene sequences suggest that the symbiont-containing Crithidia and Blastocrithidia are more akin to each other than to other trypanosomatids. The evidence suggests that Blastocrithidia and Crithidia symbionts descend from a common ancestor, which had presumably entered an ancestral host and thence coevolved with it into different species. We therefore propose naming the symbionts Kinetoplastibacterium blastocrithidii and Kinetoplastibacterium crithidii.

Evolution of RNA editing in kinetoplastid protozoa

The editing of RNA in trypanosomatid mitochondria involves the insertion and occasional deletion of uridine residues within coding regions of maxicircle messenger RNA transcripts. The extent to which the transcripts of homologous genes undergo editing differs in different species. In some, entire genes are edited (pan-editing), whereas in others, editing is limited to the 5' termini of editing domains (5' editing). Here we investigate which type of editing is ancestral and which is derived, by analysing RNA editing in the different lineages, using a kinetoplastid phylogeny reconstructed from nuclear small subunit ribosomal RNA sequences. We conclude that the ancestral cryptogenes were pan-edited, and we hypothesize that the 5'-edited homologues were generated by several independent events from partially edited RNAs, in which case editing may be a more primitive mechanism than previously thought.

Organization of mini-exon and 5S rRNA genes in the kinetoplastid Trypanoplasma borreli

Mini-exon gene repeats from Trypanoplasma borreli, which belongs to the Cryptobiidae family of the Bodonina suborder of the Kinetoplastida, were isolated by PCR amplification and cloning. The presence of kinetoplastid-like mini-exon genes in T. borreli is consistent with the taxonomic status of this organism as a kinetoplastid protozoan. Two families of repeats were found: 597 nt (T1) and 794 nt (T2), each of which encodes an approximately 95-nt medRNA transcript. The T1 repeats also contain a complete 5S rRNA gene on the complementary strand. The T2 repeats contain a defective copy of a 5S gene, in which the 5' portion is absent. The intergenic regions between the 5'-ends of the mini-exon genes and the 5S rRNA genes in the T1 and T2 repeats are highly diverged. All or most mini-exon genes and 5S genes are located within either the T1 or the T2 repeats. The T1 repeats were localized to a megabase-size chromosome, while the T2 repeats were localized within at least 4 large chromosomes.

Generation of unexpected editing patterns in Leishmania tarentolae mitochondrial mRNAs: misediting produced by misguiding

We have analyzed the generation of unexpected patterns of RNA editing, i.e., those not following a strict 3' to 5' progression, which occur in junction regions between fully edited and preedited sequences. Evidence is presented that these patterns are generated by misediting due to specific events of misguiding. Misediting can occur through the interaction of inappropriate gRNAs with mRNAs or appropriate gRNAs in an incorrect fashion. Four possible mechanisms for the generation of misedited sequences are presented. Chimeric molecules have been detected in steady-state mitochondrial RNAs that are composed of misguiding gRNAs covalently linked to mRNAs at misediting sites by the 3' oligo(U) tail. We propose that misediting within junction regions can be corrected by appropriately acting gRNAs.

The polarity of editing within a multiple gRNA-mediated domain is due to formation of anchors for upstream gRNAs by downstream editing

Seventeen kinetoplast minicircle-encoded and nine maxicircle-encoded gRNA genes have been identified. Six overlapping minicircle-encoded gRNAs mediate editing for the 5'-pan-edited MURF4 gene and two for the 5'-edited COIII gene. The pan-edited RPS12 mRNA is edited by seven minicircle-encoded gRNAs and one maxicircle-encoded gRNA. The 3'-most gRNA in each domain forms an anchor with unedited mRNA, whereas upstream gRNAs form anchors only with edited mRNA, thereby explaining the observed 3' to 5' polarity of editing within an editing domain. We suggest that a role of G-U base pairs is to allow breathing of the edited mRNA-gRNA hybrid and formation of the upstream anchor hybrid.

An intergenic G-rich region in Leishmania tarentolae kinetoplast maxicircle DNA is a pan-edited cryptogene encoding ribosomal protein S12

Six short G-rich intergenic regions in the maxicircle of Leishmania tarentolae are conserved in location and polarity in two other kinetoplastid species. We show here that G-rich region 6 (G6) represents a pan-edited cryptogene which contains at least two domains edited independently in a 3'-to-5' manner connected by short unedited regions. In the completely edited RNA, 117 uridines are added at 49 sites and 32 uridines are deleted at 13 sites, creating a translated 85-amino-acid polypeptide. Similar polypeptides are probably encoded by pan-edited G6 transcripts in two other species. The G6 polypeptide has significant sequence similarity to the family of S12 ribosomal proteins. A minicircle-encoded gRNA overlaps 12 editing sites in G6 mRNA, and chimeric gRNA/mRNA molecules were shown to exist, in agreement with the transesterification model for editing.