LOW METABOLIC RATES IN SALAMANDERS ARE CORRELATED WITH WEAK SELECTIVE CONSTRAINTS ON MITOCHONDRIAL GENES

Mitochondria are the site for the citric acid cycle and oxidative phosphorylation (OXPHOS), the final steps of ATP synthesis via cellular respiration. Each mitochondrion contains its own genome; in vertebrates, this is a small, circular DNA molecule that encodes 13 subunits of the multiprotein OXPHOS electron transport complexes. Vertebrate lineages vary dramatically in metabolic rates; thus, functional constraints on mitochondrial‐encoded proteins likely differ, potentially impacting mitochondrial genome evolution. Here, we examine mitochondrial genome evolution in salamanders, which have the lowest metabolic requirements among tetrapods. We show that salamanders experience weaker purifying selection on protein‐coding sequences than do frogs, a comparable amphibian clade with higher metabolic rates. In contrast, we find no evidence for weaker selection against mitochondrial genome expansion in salamanders. Together, these results suggest that different aspects of mitochondrial genome evolution (i.e., nucleotide substitution, accumulation of noncoding sequences) are differently affected by metabolic variation across tetrapod lineages.     

Genetic drift and mutational hazard in the evolution of salamander genomic gigantism

Salamanders have the largest nuclear genomes among tetrapods and, excepting lungfishes, among vertebrates as a whole. Lynch and Conery (2003) have proposed the mutational‐hazard hypothesis to explain variation in genome size and complexity. Under this hypothesis, noncoding DNA imposes a selective cost by increasing the target for degenerative mutations (i.e., the mutational hazard). Expansion of noncoding DNA, and thus genome size, is driven by increased levels of genetic drift and/or decreased mutation rates; the former determines the efficiency with which purifying selection can remove excess DNA, whereas the latter determines the level of mutational hazard. Here, we test the hypothesis that salamanders have experienced stronger long‐term, persistent genetic drift than frogs, a related clade with more typically sized vertebrate genomes. To test this hypothesis, we compared dN/dS and Kr/Kc values of protein‐coding genes between these clades. Our results do not support this hypothesis; we find that salamanders have not experienced stronger genetic drift than frogs. Additionally, we find evidence consistent with a lower nucleotide substitution rate in salamanders. This result, along with previous work showing lower rates of small deletion and ectopic recombination in salamanders, suggests that a lower mutational hazard may contribute to genomic gigantism in this clade.     

Incorporation of transition to transversion ratio and nonsense mutations,improves the estimation of the number of synonymous and non-synonymous sites in codons

A common approach to estimate the strength and direction of selection acting on protein coding sequences is to calculate the dN/dS ratio. The method to calculate dN/dS has been widely used by many researchers and many critical reviews have been made on its application after the proposition by Nei and Gojobori in 1986. However, the method is still evolving considering the non-uniform substitution rates and pretermination codons. In our study of SNPs in 586 genes across 156 Escherichia coli strains, synonymous polymorphism in 2-fold degenerate codons were higher in comparison to that in 4-fold degenerate codons, which could be attributed to the difference between transition (Ti) and transversion (Tv) substitution rates where the average rate of a transition is four times more than that of a transversion in general. We considered both the Ti/Tv ratio, and nonsense mutation in pretermination codons, to improve estimates of synonymous (S) and non-synonymous (NS) sites. The accuracy of estimating dN/dS has been improved by considering the Ti/Tv ratio and nonsense substitutions in pretermination codons. We showed that applying the modified approach based on Ti/Tv ratio and pretermination codons results in higher values of dN/dS in 29 common genes of equal reading-frames between E. coli and Salmonella enterica. This study emphasizes the robustness of amino acid composition with varying codon degeneracy, as well as the pretermination codons when calculating dN/dS values.     

Patterns of evolutionary selection pressure in the immune signaling protein TRAF3IP2 in mammals

TRAF3 interacting protein 2 (TRAF3IP2) is important for immune responses to pathogens, inflammatory signals and autoimmunity in mammals. In the present study, we collected 19 mammalian TRAF3IP2 sequences and investigated the various types of selection pressure acting on them. Maximum likelihood estimations of nonsynonymous (dN) to synonymous (dS) substitution (dN/dS) ratios for the aligned coding sequences indicated that, as a whole, TRAF3IP2 has been subject to purifying selection. However, the N-terminus of the protein has been subject to higher selection pressure than the C-terminal domain. While eight amino acid residues within the N-terminus appear to have evolved under positive selection, no evidence for such selection was found in the C-terminus. The positively selected residues, which fall outside the currently known functional sites within TRAF3IP2, may have novel functions. The different selection pressures acting on the N- and C-terminal regions are consistent with their protein structures: the C-terminal structure is an ordered structure, whereas the N-terminus is disordered. Taken together with the results of previous studies, it is plausible that positive selection on the N-terminus of TRAF3IP2 may have occurred by competitive coevolution between mammalian hosts and viruses.     

Whole-Genome Sequencing of Theileria parva Strains Provides Insight into Parasite Migration and Diversification in the African Continent

The disease caused by the apicomplexan protozoan parasite Theileria parva, known as East Coast fever or Corridor disease, is one of the most serious cattle diseases in Eastern, Central, and Southern Africa. We performed whole-genome sequencing of nine T. parva strains, including one of the vaccine strains (Kiambu 5), field isolates from Zambia, Uganda, Tanzania, or Rwanda, and two buffalo-derived strains. Comparison with the reference Muguga genome sequence revealed 34 814–121 545 single nucleotide polymorphisms (SNPs) that were more abundant in buffalo-derived strains. High-resolution phylogenetic trees were constructed with selected informative SNPs that allowed the investigation of possible complex recombination events among ancestors of the extant strains. We further analysed the dN/dS ratio (non-synonymous substitutions per non-synonymous site divided by synonymous substitutions per synonymous site) for 4011 coding genes to estimate potential selective pressure. Genes under possible positive selection were identified that may, in turn, assist in the identification of immunogenic proteins or vaccine candidates. This study elucidated the phylogeny of T. parva strains based on genome-wide SNPs analysis with prediction of possible past recombination events, providing insight into the migration, diversification, and evolution of this parasite species in the African continent.     

Structure-function analysis of a bacterial deoxyadenosine kinase reveals the basis for substrate specificity

Deoxyribonucleoside kinases (dNKs) catalyze the transfer of a phosphoryl group from ATP to a deoxyribonucleoside (dN), a key step in DNA precursor synthesis. Recently structural information concerning dNKs has been obtained, but no structure of a bacterial dCK/dGK enzyme is known. Here we report the structure of such an enzyme, represented by deoxyadenosine kinase from Mycoplasma mycoides subsp. mycoides small colony type (Mm-dAK). Superposition of Mm-dAK with its human counterpart's deoxyguanosine kinase (dGK) and deoxycytidine kinase (dCK) reveals that the overall structures are very similar with a few amino acid alterations in the proximity of the active site. To investigate the substrate specificity, Mm-dAK has been crystallized in complex with dATP and dCTP, as well as the products dCMP and dCDP. Both dATP and dCTP bind to the enzyme in a feedback-inhibitory manner with the dN part in the deoxyribonucleoside binding site and the triphosphates in the P-loop. Substrate specificity studies with clinically important nucleoside analogs as well as several phosphate donors were performed. Thus, in this study we combine structural and kinetic data to gain a better understanding of the substrate specificity of the dCK/dGK family of enzymes. The structure of Mm-dAK provides a starting point for making new anti bacterial agents against pathogenic bacteria.     

Dictyostelium discoideum salvages purine deoxyribonucleosides by highly specific bacterial-like deoxyribonucleoside kinases

The salvage of deoxyribonucleosides in the social amoeba Dictyostelium discoideum, which has an extremely A+T-rich genome, was investigated. All native deoxyribonucleosides were phosphorylated by D. discoideum cell extracts and we subcloned three deoxyribonucleoside kinase (dNK) encoding genes. D. discoideum thymidine kinase was similar to the human thymidine kinase 1 and was specific for thymidine with a K(m) of 5.1 microM. The other two cloned kinases were phylogenetically closer to bacterial deoxyribonucleoside kinases than to the eukaryotic enzymes. D. discoideum deoxyadenosine kinase (DddAK) had a K(m) for deoxyadenosine of 22.7 microM and a k(cat) of 3.7 s(-1) and could not efficiently phosphorylate any other native deoxyribonucleoside. D. discoideum deoxyguanosine kinase was also a purine-specific kinase and phosphorylated significantly only deoxyguanosine, with a K(m) of 1.4 microM and a k(cat) of 3 s(-1). The two purine-specific deoxyribonucleoside kinases could represent ancient enzymes present in the common ancestor of bacteria and eukaryotes but remaining only in a few eukaryote lineages. The narrow substrate specificity of the D. discoideum dNKs reflects the biased genome composition and we attempted to explain the strict preference of DddAK for deoxyadenosine by modeling the active center with different substrates. Apart from its native substrate, deoxyadenosine, DddAK efficiently phosphorylated fludarabine. Hence, DddAK could be used in the enzymatic production of fludarabine monophosphate, a drug used in the treatment of chronic lymphocytic leukemia.     

BASE: A novel workflow to integrate nonubiquitous genes in comparative genomics analyses for selection

Inferring the selective forces that orthologous genes underwent across different lineages can help us understand the evolutionary processes that have shaped their extant diversity and the phenotypes they underlie. The most widespread metric to estimate the selection regimes of coding genes—across sites and phylogenies—is the ratio of nonsynonymous to synonymous substitutions (dN/dS, also known as ω). Nowadays, modern sequencing technologies and the large amount of already available sequence data allow the retrieval of thousands of orthologous genes across large numbers of species. Nonetheless, the tools available to explore selection regimes are not designed to automatically process all genes, and their practical usage is often restricted to the single‐copy ones which are found across all species considered (i.e., ubiquitous genes). This approach limits the scale of the analysis to a fraction of single‐copy genes, which can be as low as an order of magnitude in respect to those which are not consistently found in all species considered (i.e., nonubiquitous genes). Here, we present a workflow named BASE that—leveraging the CodeML framework—eases the inference and interpretation of gene selection regimes in the context of comparative genomics. Although a number of bioinformatics tools have already been developed to facilitate this kind of analyses, BASE is the first to be specifically designed to allow the integration of nonubiquitous genes in a straightforward and reproducible manner. The workflow—along with all relevant documentation—is available at github.com/for‐giobbe/BASE.