Convergence of immune escape strategies highlights plasticity of SARS-CoV-2 spike

The global spread of the SARS-CoV-2 virus has resulted in emergence of lineages which impact the effectiveness of immunotherapies and vaccines that are based on the early Wuhan isolate. All currently approved vaccines employ the spike protein S, as it is the target for neutralizing antibodies. Here we describe two SARS-CoV-2 isolates with unusually large deletions in the N-terminal domain (NTD) of the spike. Cryo-EM structural analysis shows that the deletions result in complete reshaping of the NTD supersite, an antigenically important region of the NTD. For both spike variants the remodeling of the NTD negatively affects binding of all tested NTD-specific antibodies in and outside of the NTD supersite. For one of the variants, we observed a P9L mediated shift of the signal peptide cleavage site resulting in the loss of a disulfide-bridge; a unique escape mechanism with high antigenic impact. Although the observed deletions and disulfide mutations are rare, similar modifications have become independently established in several other lineages, indicating a possibility to become more dominant in the future. The observed plasticity of the NTD foreshadows its broad potential for immune escape with the continued spread of SARS-CoV-2.

A Therapeutic Hepatitis B Virus DNA Vaccine Induces Specific Immune Responses in Mice and Non-Human Primates

Despite the availability of an effective prophylactic vaccine for more than 30 years, nearly 300 million people worldwide are chronically infected with the hepatitis B virus (HBV), leading to 1 death every 30 s mainly from viral hepatitis-related cirrhosis and liver cancer. Chronic HBV patients exhibit weak, transient, or dysfunctional CD8⁺ T-cell responses to HBV, which contrasts with high CD8+ T-cell responses seen for resolvers of acute HBV infection. Therefore, a therapeutic DNA vaccine was designed, expressing both HBV core and polymerase proteins, and was sequence optimized to ensure high protein expression and secretion. Although the vaccine, administered intramuscularly via electroporation, had no effect on plasma viral parameters in a mouse model of persistent HBV infection, it did induce robust HBV-specific immune responses in healthy and adeno-associated hepatitis B virus (AAV-HBV) infected mice as well as in healthy non-human primates.

Tales of 1,008 small molecules: phenomic profiling through live-cell imaging in a panel of reporter cell lines

Phenomic profiles are high-dimensional sets of readouts that can comprehensively capture the biological impact of chemical and genetic perturbations in cellular assay systems. Phenomic profiling of compound libraries can be used for compound target identification or mechanism of action (MoA) prediction and other applications in drug discovery. To devise an economical set of phenomic profiling assays, we assembled a library of 1,008 approved drugs and well-characterized tool compounds manually annotated to 218 unique MoAs, and we profiled each compound at four concentrations in live-cell, high-content imaging screens against a panel of 15 reporter cell lines, which expressed a diverse set of fluorescent organelle and pathway markers in three distinct cell lineages. For 41 of 83 testable MoAs, phenomic profiles accurately ranked the reference compounds (AUC-ROC ≥ 0.9). MoAs could be better resolved by screening compounds at multiple concentrations than by including replicates at a single concentration. Screening additional cell lineages and fluorescent markers increased the number of distinguishable MoAs but this effect quickly plateaued. There remains a substantial number of MoAs that were hard to distinguish from others under the current study's conditions. We discuss ways to close this gap, which will inform the design of future phenomic profiling efforts.

Anti-HBV activity of the HBV capsid assembly modulator JNJ-56136379 across full-length genotype A–H clinical isolates and core site-directed mutants in vitro

Objectives

To characterize antiviral activity of the capsid assembly modulator (CAM-N) JNJ-56136379 against HBV genotypes and variants carrying amino acid substitutions in the core protein.

Methods

Anti-HBV activity of JNJ-56136379 was investigated against a diverse panel of 53 HBV clinical isolates (genotypes A–H). The impact of core amino acid substitutions using site-directed mutants (SDMs) was assessed in a transient replication assay.

Results

JNJ-56136379 median 50% effective concentration (EC50) values across all genotypes were 10–33 nM versus 17 nM (genotype D reference). JNJ-56136379 remained active against isolates carrying nucleos(t)ide analogue resistance mutations (median EC50 2–25 nM) or basal core promoter (BCP) ± precore (PC) mutations (median EC50 13–20 nM) or PC mutations (median EC50 11 nM), representing activity against isolates from HBeAg-positive and -negative hepatitis B patients. Core amino acid substitutions in the CAM-binding pocket, when tested as SDMs at positions 23, 25, 30, 33, 37, 106, 110, 118, 124, 127 and 128, reduced JNJ-56136379 anti-HBV activity; EC50 fold increases ranged from 3.0 (S106T) to 85 (T33N). All substitutions were rare in a public database of >7600 HBV core sequences (frequencies 0.01%–0.3%). Nucleos(t)ide analogues retained full activity against these core SDMs.

Conclusions

JNJ-56136379, a potent HBV CAM-N, currently in Phase 2 clinical development, was generally fully active against an extensive panel of genotype A–H clinical isolates, regardless of the presence of nucleos(t)ide analogue resistance or BCP/PC mutations. JNJ-56136379 activity was reduced by some core amino acid substitutions in the CAM-binding pocket.

The Need of Industry to Go FAIR

The industry sector is a very large producer and consumer of data, and many companies traditionally focused on production or manufacturing are now relying on the analysis of large amounts of data to develop new products and services. As many of the data sources needed are distributed and outside the company, FAIR data will have a major impact, both by reducing the existing internal data silos and by enabling the efficient integration with external (public and commercial) data. Many companies are still in the early phases of internal data “FAIRification”, providing opportunities for SMEs and academics to apply and develop their expertise on FAIR data in collaborations and public-private partnerships. For a global Internet of FAIR Data & Services to thrive, also involving industry, professional tools and services are essential. FAIR metrics and certifications on individuals, data, organizations, and software, must ensure that data producers and consumers have independent quality metrics on their data. In this opinion article we reflect on some industry specific challenges of FAIR implementation to be dealt with when choices are made regarding “Industry GOing FAIR”.

Protocols for the Design of Kinase-focused Compound Libraries

Protocols for the design of kinase-focused compound libraries are presented. Kinase-focused compound libraries can be differentiated based on the design goal. Depending on whether the library should be a discovery library specific for one particular kinase, a general discovery library for multiple distinct kinase projects, or even phenotypic screening, there exists today a variety of in silico methods to design candidate compound libraries. We address the following scenarios: 1) Datamining of SAR databases and kinase focused vendor catalogues; 2) Predictions and virtual screening; 3) Structure-based design of combinatorial kinase inhibitors; 4) Design of covalent kinase inhibitors; 5) Design of macrocyclic kinase inhibitors; and 6) Design of allosteric kinase inhibitors and activators.

Extending kinome coverage by analysis of kinase inhibitor broad profiling data

The explored kinome was extended with broad profiling using the DiscoveRx and Millipore assay panels. The analysis of the profiling of 3368 selected inhibitors on 456 kinases in the DiscoveRx format delivered several insights. First, the coverage depended on the threshold of the selectivity parameter. Second, the relation between hit confirmation rates and inhibitor selectivity showed unexpectedly that higher selectivity can increase the likelihood of false positives. Third, comparing the coverage of a focused to that of a random library showed that the design based on a maximum number of scaffolds was superior to a limited number of scaffolds. Therefore, selective compounds can be used in target validation, enable the jumpstarting of new kinase drug discovery projects, and chart new biological space via phenotypic screening.

The application of the open pharmacological concepts triple store (open PHACTS) to support drug discovery research

Integration of open access, curated, high-quality information from multiple disciplines in the Life and Biomedical Sciences provides a holistic understanding of the domain. Additionally, the effective linking of diverse data sources can unearth hidden relationships and guide potential research strategies. However, given the lack of consistency between descriptors and identifiers used in different resources and the absence of a simple mechanism to link them, gathering and combining relevant, comprehensive information from diverse databases remains a challenge. The Open Pharmacological Concepts Triple Store (Open PHACTS) is an Innovative Medicines Initiative project that uses semantic web technology approaches to enable scientists to easily access and process data from multiple sources to solve real-world drug discovery problems. The project draws together sources of publicly-available pharmacological, physicochemical and biomolecular data, represents it in a stable infrastructure and provides well-defined information exploration and retrieval methods. Here, we highlight the utility of this platform in conjunction with workflow tools to solve pharmacological research questions that require interoperability between target, compound, and pathway data. Use cases presented herein cover 1) the comprehensive identification of chemical matter for a dopamine receptor drug discovery program 2) the identification of compounds active against all targets in the Epidermal growth factor receptor (ErbB) signaling pathway that have a relevance to disease and 3) the evaluation of established targets in the Vitamin D metabolism pathway to aid novel Vitamin D analogue design. The example workflows presented illustrate how the Open PHACTS Discovery Platform can be used to exploit existing knowledge and generate new hypotheses in the process of drug discovery.

Delayed bactericidal response of Mycobacterium tuberculosis to bedaquiline involves remodelling of bacterial metabolism

Bedaquiline (BDQ), an ATP synthase inhibitor, is the first drug to be approved for treatment of multidrug-resistant tuberculosis in decades. Though BDQ has shown excellent efficacy in clinical trials, its early bactericidal activity during the first week of chemotherapy is minimal. Here, using microfluidic devices and time-lapse microscopy of Mycobacterium tuberculosis, we confirm the absence of significant bacteriolytic activity during the first 3-4 days of exposure to BDQ. BDQ-induced inhibition of ATP synthesis leads to bacteriostasis within hours after drug addition. Transcriptional and proteomic analyses reveal that M. tuberculosis responds to BDQ by induction of the dormancy regulon and activation of ATP-generating pathways, thereby maintaining bacterial viability during initial drug exposure. BDQ-induced bacterial killing is significantly enhanced when the mycobacteria are grown on non-fermentable energy sources such as lipids (impeding ATP synthesis via glycolysis). Our results show that BDQ exposure triggers a metabolic remodelling in mycobacteria, thereby enabling transient bacterial survival.

Nociceptive and behavioural sensitisation by protein kinase Cepsilon signalling in the CNS

Despite the apparent homology in the protein kinase C (PKC) family, it has become clear that slight structural differences are sufficient to have unique signalling properties for each individual isoform. For PKCepsilon in depth investigation of these aspects revealed unique actions in the CNS and lead to development of specific modulators with clinical perspective. In this review, we describe to which extent PKCepsilon is distinct from other isoforms on the level of tissue expression and protein structure. As this kinase is highly expressed in the brain, we outline three main aspects of PKCepsilon signalling in the CNS. First, its ability to alter the permeability of N-type Ca2+ channels in dorsal root ganglia has been shown to enhance nociception. Secondly, PKCepsilon increases anxiety by diminishing GABA(A)R-induced inhibitory post-synaptic currents in the prefrontal cortex. Another important aspect of the latter inhibition is the reduced sensitivity of GABA(A) receptors to ethanol, a mechanism potentially contributing to abuse. A third signalling cascade improves cognitive functions by facilitating cholinergic signalling in the hippocampus. Collectively, these findings point to a physical and behavioural sensitising role for this kinase.

Downregulation of two novel genes in Sl/Sld and W(LacZ)/Wv mouse jejunum

Interstitial cells of Cajal (ICC) are the so-called pacemaker cells of the gut. W(LacZ)/Wv and Sl/Sld mice lack ICC surrounding the myenteric plexus (MP) in the jejunum. We compared the gene expression profile of wild type (WT) and W(LacZ)/Wv and Sl/Sld mice using suppression subtractive hybridization (SSH), generating a cDNA library of 1303 clones from which 48 unique sequences were differentially expressed with Southern blot. Among them, we identified heme oxygenase2, TROY, and phospholamban in ICC using immunohistochemistry. Using RT-qPCR, c-Kit and two new transcripts Dithp and prenylcysteine oxidase1 were significantly lower expressed in Sl/Sld and W(LacZ)/Wv versus WT. Prenylcysteine oxidase1 appeared cytotoxic for COS-7 cells and was highly expressed in liver while Dithp was mainly expressed in small intestine. The combination of SSH, Southern blot, RT-qPCR, and immunohistochemistry turned out to be a useful approach to identify rarely expressed genes and genes with small differences in expression.

A diarylquinoline drug active on the ATP synthase of Mycobacterium tuberculosis

The incidence of tuberculosis has been increasing substantially on a worldwide basis over the past decade, but no tuberculosis-specific drugs have been discovered in 40 years. We identified a diarylquinoline, R207910, that potently inhibits both drug-sensitive and drug-resistant Mycobacterium tuberculosis in vitro (minimum inhibitory concentration 0.06 mug/ml). In mice, R207910 exceeded the bactericidal activities of isoniazid and rifampin by at least 1 log unit. Substitution of drugs included in the World Health Organization's first-line tuberculosis treatment regimen (rifampin, isoniazid, and pyrazinamide) with R207910 accelerated bactericidal activity, leading to complete culture conversion after 2 months of treatment in some combinations. A single dose of R207910 inhibited mycobacterial growth for 1 week. Plasma levels associated with efficacy in mice were well tolerated in healthy human volunteers. Mutants selected in vitro suggest that the drug targets the proton pump of adenosine triphosphate (ATP) synthase.

Hevin is down-regulated in many cancers and is a negative regulator of cell growth and proliferation

We have cloned a human Hevin cDNA from omental adipose tissue of different patients by reverse transcription polymerase chain reaction and shown a sequence variation due to a possible polymorphism at amino acid position 161 (E/G). Hevin protein expressed in vitro showed molecular weights of approximately 75 kDa and 150 kDa, suggesting that Hevin may form a homodimer in vitro. Using Northern blots and a human expressed sequence tAg database analysis, Hevin was shown to be widely expressed in human normal or non-neoplastic diseased tissues with various levels. In contrast to this, its expression was strongly down-regulated in most neoplastic cells or tissues tested. However, neither the mechanism nor the physiological meaning of this down-regulation is known. As an initial step towards investigating the functional role of Hevin in cell growth and differentiation, we transiently or stably expressed this gene in cancer cells (HeLa 3S) that are devoid of endogenous Hevin and measured DNA synthesis (cell proliferation) by 5-bromo-2'-deoxyuridine incorporation. Hevin was shown to be a negative regulator of cell proliferation. Furthermore, we have shown that Hevin can inhibit progression of cells from G1 to S phase or prolong G1 phase. This is the first report which describes the function of Hevin in cell growth and proliferation. Through database analysis, Hevin was found to be located on chromosome 4 which contains loss of heterozygosity of many tumour suppressor genes. Taken together, these results suggest that Hevin may be a candidate for a tumour suppressor gene and a potential target for cancer diagnosis/therapy.

Isolation and expression of novel human glutamate carboxypeptidases with N-acetylated alpha-linked acidic dipeptidase and dipeptidyl peptidase IV activity

Hydrolysis of the neuropeptide N-acetyl-L-aspartyl-L-glutamate (NAAG) by N-acetylated alpha-linked acidic dipeptidase (NAALADase) to release glutamate may be important in a number of neurodegenerative disorders in which excitotoxic mechanisms are implicated. The gene coding for human prostate-specific membrane antigen, a marker of prostatic carcinomas, and its rat homologue glutamate carboxypeptidase II have recently been shown to possess such NAALADase activity. In contrast, a closely related member of this gene family, rat ileal 100-kDa protein, possesses a dipeptidyl peptidase IV activity. Here, we describe the cloning of human ileal 100-kDa protein, which we have called a NAALADase- "like" (NAALADase L) peptidase based on its sequence similarity to other members of this gene family, and its inability to hydrolyze NAAG in transient transfection experiments. Furthermore, we describe the cloning of a third novel member of this gene family, NAALADase II, which codes for a type II integral membrane protein and which we have localized to chromosome 11 by fluorescent in situ hybridization analysis. Transient transfection of NAALADase II cDNA confers both NAALADase and dipeptidyl peptidase IV activity to COS cells. Expression studies using reverse transcription-polymerase chain reaction and Northern blot hybridization show that NAALADase II is highly expressed in ovary and testis as well as within discrete brain areas.

Structure of the human histamine H1 receptor gene

Histamine H1 receptor expression has been reported to change in disorders such as allergic rhinitis, autoimmune myocarditis, rheumatoid arthritis and atherosclerosis. Here we report the isolation and characterization of genomic clones containing the 5' flanking (regulatory) region of the human histamine H1 receptor gene. An intron of approx. 5.8 kb was identified in the 5' untranslated region, which suggests that an entire subfamily of G-protein-coupled receptors may contain an intron immediately upstream of the start codon. The transcription initiation site was mapped by 5' rapid amplification of cDNA ends to a region 6.2 kb upstream of the start codon. Immediately upstream of the transcription start site a fragment of 1.85 kb was identified that showed promoter activity when placed upstream of a luciferase reporter gene and transiently transfected into cells expressing the histamine H1 receptor. The promoter sequence shares a number of characteristics with the promoter sequences of other G-protein-coupled receptor encoding genes, including binding sites for several transcription factors, and the absence of TATA and CAAT sequences at the appropriate locations. The promoter sequence described here differs from that reported previously [Fukui, Fujimoto, Mizuguchi, Sakamoto, Horio, Takai, Yamada and Ito (1994) Biochem. Biophys. Res. Commun. 201, 894-901] because the reported genomic clone was chimaeric. Furthermore our study provides evidence that the 3' untranslated region of the H1 receptor mRNA is much longer than previously accepted. Together, these findings provide a complete view of the structure of the human histamine H1 receptor gene. Both the coding region of the H1 receptor gene and its promoter region were independently mapped to chromosome 3p25.

Evolutionary affiliation of the marine nitrogen-fixing cyanobacterium Trichodesmium sp. strain NIBB 1067, derived by 16S ribosomal RNA sequence analysis

The 16S rRNA sequence of Trichodesmium sp. strain NIBB 1067 was determined and used for the construction of a distance tree and bootstrap analysis. The tree shows that, among the available cyanobacterial 16S rRNA sequences, Trichodesmium NIBB 1067 has Oscillatoria PCC 7515 as its closest relative, presenting 94.9% of sequence similarity with the latter strain. This is in contrast to a difference of 9 mol% G+C in mean genomic DNA base composition between the two organisms. Nevertheless, the genotypic heterogeneity presented by a number of strains assigned to the genus Oscillatoria hinders a taxonomic decision on the separate existence of the genera Trichodesmium and Oscillatoria. The sequence of the internal transcribed spacer (ITS) between the 16S and 23S rRNA genes was also determined, as a possible marker to study inter- and intraspecific variability. The ITS contains the genes coding for tRNA(Ile) and tRNA(Ala) and its total length is 547 nucleotides. In six out of eight sequenced clones, there is a duplication of 29 nucleotides, surrounding the 5' end of the tRNA(Ile).

Reconstructing evolution from eukaryotic small-ribosomal-subunit RNA sequences: calibration of the molecular clock

The detailed descriptions now available for the secondary structure of small-ribosomal-subunit RNA, including areas of highly variable primary structure, facilitate the alignment of nucleotide sequences. However, for optimal exploitation of the information contained in the alignment, a method must be available that takes into account the local sequence variability in the computation of evolutionary distance. A quantitative definition for the variability of an alignment position is proposed in this study. It is a parameter in an equation which expresses the probability that the alignment position contains a different nucleotide in two sequences, as a function of the distance separating these sequences, i.e., the number of substitutions per nucleotide that occurred during their divergence. This parameter can be estimated from the distance matrix resulting from the conversion of pairwise sequence dissimilarities into pairwise distances. Alignment positions can then be subdivided into a number of sets of matching variability, and the average variability of each set can be derived. Next, the conversion of dissimilarity into distance can be recalculated for each set of alignment positions separately, using a modified version of the equation that corrects for multiple substitutions and changing for each set the parameter that reflects its average variability. The distances computed for each set are finally averaged, giving a more precise distance estimation. Trees constructed by the algorithm based on variability calibration have a topology markedly different from that of trees constructed from the same alignments in the absence of calibration. This is illustrated by means of trees constructed from small-ribosomal-subunit RNA sequences of Metazoa. A reconstruction of vertebrate evolution based on calibrated alignments matches the consensus view of paleontologists, contrary to trees based on uncalibrated alignments. In trees derived from sequences covering several metazoan phyla, artefacts in topology that are probably due to a high clock rate in certain lineages are avoided.

Compilation of small ribosomal subunit RNA structures

The database on small ribosomal subunit RNA structure contained 1804 nucleotide sequences on April 23, 1993. This number comprises 365 eukaryotic, 65 archaeal, 1260 bacterial, 30 plastidial, and 84 mitochondrial sequences. These are stored in the form of an alignment in order to facilitate the use of the database as input for comparative studies on higher-order structure and for reconstruction of phylogenetic trees. The elements of the postulated secondary structure for each molecule are indicated by special symbols. The database is available on-line directly from the authors by ftp and can also be obtained from the EMBL nucleotide sequence library by electronic mail, ftp, and on CD ROM disk.

The gene coding for small ribosomal subunit RNA in the basidiomycete Ustilago maydis contains a group I intron

The nucleotide sequence of the gene coding for small ribosomal subunit RNA in the basidiomycete Ustilago maydis was determined. It revealed the presence of a group I intron with a length of 411 nucleotides. This is the third occurrence of such an intron discovered in a small subunit rRNA gene encoded by a eukaryotic nuclear genome. The other two occurrences are in Pneumocystis carinii, a fungus of uncertain taxonomic status, and Ankistrodesmus stipitatus, a green alga. The nucleotides of the conserved core structure of 101 group I intron sequences present in different genes and genome types were aligned and their evolutionary relatedness was examined. This revealed a cluster including all group I introns hitherto found in eukaryotic nuclear genes coding for small and large subunit rRNAs. A secondary structure model was designed for the area of the Ustilago maydis small ribosomal subunit RNA precursor where the intron is situated. It shows that the internal guide sequence pairing with the intron boundaries fits between two helices of the small subunit rRNA, and that minimal rearrangement of base pairs suffices to achieve the definitive secondary structure of the 18S rRNA upon splicing.

Phylogenetic analysis of five medically important Candida species as deduced on the basis of small ribosomal subunit RNA sequences

The classification of species belonging to the genus Candida Berkhout is problematic. Therefore, we have determined the small ribosomal subunit RNA (srRNA) sequences of the type strains of three human pathogenic Candida species; Candida krusei, C. lusitaniae and C. tropicalis. The srRNA sequences were aligned with published eukaryotic srRNA sequences and evolutionary trees were inferred using a matrix optimization method. An evolutionary tree comprising all available eukaryotic srRNA sequences, including two other pathogenic Candida species, C. albicans and C. glabrata, showed that the yeasts diverge rather late in the course of eukaryote evolution, namely at the same depth as green plants, ciliates and some smaller taxa. The cluster of the higher fungi consists of 10 ascomycetes and ascomycete-like species with the first branches leading to Neurospora crassa, Pneumocystis carinii, Candida lusitaniae and C. krusei, in that order. Next there is a dichotomous divergence leading to a group consisting of Torulaspora delbrueckii, Saccharomyces cerevisiae, C. glabrata and Kluyveromyces lactis and a smaller group comprising C. tropicalis and C. albicans. The divergence pattern obtained on the basis of srRNA sequence data is also compared to various other chemotaxonomic data.

A proposal for the secondary structure of a variable area of eukaryotic small ribosomal subunit RNA involving the existence of a pseudoknot

Eukaryotic small ribosomal subunit RNAs contain an area of variable structure, V4, which comprises about 250 nucleotides in most species, whereas the corresponding area in bacterial small ribosomal subunit RNAs consists of about 64 nucleotides folded into a single hairpin. There is no consensus on the secondary structure of area V4 in eukaryotes, about 10 different models having been proposed. The prediction of a model on a comparative basis poses special problems because, due to the variability of the area in length as well as sequence, a dependable alignment is very difficult to achieve. A new model was derived by systematic examination of all combinations of helices that have been hitherto proposed, plus some new ones. The following properties of the helices were examined: transposability to all presently known sequences, presence of compensating substitutions, and thermodynamic stability. A model was selected by ranking all possible combinations of transposable helices according to the number of compensating substitutions scored. The optimal model comprises a pseudoknot and four hairpin structures. Certain species contain additional hairpins inserted between these structural elements, while in others the structure is partially or entirely deleted.

The 18S ribosomal RNA sequence of the sea anemone Anemonia sulcata and its evolutionary position among other eukaryotes

Evolutionary trees based on partial small ribosomal subunit RNA sequences of 22 metazoa species have been published [(1988) Science 239, 748-753]. In these trees, cnidarians (Radiata) seemed to have evolved independently from the Bilateria, which is in contradiction with the general evolutionary view. In order to further investigate this problem, the complete srRNA sequence of the sea anemone Anemonia sulcata was determined and evolutionary trees were constructed using a matrix optimization method. In the tree thus obtained the sea anemone and Bilateria together form a monophyletic cluster, with the sea anemone forming the first line of the metazoan group.

Small ribosomal subunit RNA sequences, evolutionary relationships among different life forms, and mitochondrial origins

A tree was constructed from a structurally conserved area in an alignment of 83 small ribosomal subunit sequences of eukaryotic, archaebacterial, eubacterial, plastidial, and mitochondrial origin. The algorithm involved computation and optimization of a dissimilarity matrix. According to the tree, only plant mitochondria belong to the eubacterial primary kingdom, whereas animal, fungal, algal, and ciliate mitochondria branch off from an internal node situated between the tree primary kingdoms. This result is at variance with a parsimony tree of similar size published by Cedergren et al. (J Mol Evol 28:98-112, 1988), which postulates the mitochondria to be monophyletic and to belong to the eubacterial primary kingdom. The discrepancy does not follow from the use of conflicting sequence alignments, hence it must be due to the use of different treeing algorithms. We tested our algorithm on a set of sequences resulting from a simulated evolution and found it capable of faithfully reconstructing a branching topology that involved very unequal evolutionary rates. The use of more limited or more extended areas of the complete sequence alignment, comprising only very conserved or also more variable portions of the small ribosomal subunit structure, does have some influence on the tree topology. In all cases, however, the nonplant mitochondria seem to branch off before the emergence of eubacteria, and the differences are limited to the branching pattern among different types of mitochondria.