Radical Amino Acid Substitutions Research Articles

Analysis of single nucleotide polymorphisms in sequences of Pina and Pinb genes of the diploid wheats Triticum monococcum and T. urartu

Puroindolines (puroindoline a and puroindoline b) are low molecular weight proteins that determine the endosperm texture of grain in the tribes Triticeae and Avenae. The aim of our study was to analyse single nucleotide polymorphisms (SNPs) of the coding sequences of the Pina and Pinb genes of the diploid wheats Triticum monococcum (AmAm) and T. urartu (AA) from the NCBI database in comparison with the reference sequences of the common wheat variety Chinese Spring (CS).For this, 62 sequences of the T. monococcum puroindoline a gene (Pina-Am1), 22 sequences of T. urartu Pina (Pina-A1), 32 sequences of the T. monococcum puroindoline b gene (Pinb-Am1) and 13 sequences of T. urartu Pinb (Pinb-A1) were retrieved from the NCBI database. The sequences of the puroindoline a gene DQ363911.1 of the variety CS (the Pina-D1a allele) and the puroindoline b gene DQ363913.1 of CS (the Pinb-D1a allele) from the NCBI database were used as the reference sequences. The sequences were aligned using MEGA 11. In total, 34 SNPs (13 synonymous and 21 nonsynonymous differences, of which 15 result in radical amino acid substitutions and 6 in conservative ones) were identified in the total sample of 84 diploid wheat Pina sequences, some of which were observed in all sequences and some were rare. Among the 45 Pinb sequences, there were 36 SNPs, but, unlike the Pina gene, synonymous substitutions prevailed (22); 7 substitutions led to radical amino acid substitutions and 7 to conservative ones. Substitutions in the sequences of the puroindoline genes relative to the CS genes can be divided into those fixed in both diploid wheat species, those fixed in T. urartu and polymorphic in T. monococcum, and species-specific ones. Significant differences in the frequencies of alternative nucleotides at certain positions (81, 318, 322 and 384 of Pina and 135 and 359 of Pinb) were found between the wild einkorn wheat T. monococcum ssp. aegilopoides and the cultivated wheat T. monococcum ssp. monococcum.

Defensins of the stable fly (Stomoxys Calcitrans) have developmental-specific regulation and evolve at different rates.

Organisms produce antimicrobial peptides (AMPs) either in response to infection (induced) or continuously (constitutively) to combat microbes encountered during normal trophic activities and/or through pathogenic infections. The expression of AMPs is tightly regulated often with specificity to particular tissues or developmental stages. As a result, AMPs face varying selective pressures based on the microbes the organism's tissue or developmental stage encounters. Here, we analyzed the evolution and developmental-specific expression of Defensins, which are ancient AMPs in insects, in the stable fly (Stomoxys calcitrans). Stable fly larvae inhibit microbe-rich environments, whereas adult flies, as blood-feeders, experience comparatively fewer encounters with diverse microbial communities. Using existing RNA-seq datasets, we identified six Defensins that were only expressed in larvae (larval Defensins) and five that were not expressed in larvae (non-larval Defensins). Each of the non-larval Defensins was expressed in at least one adult tissue sample. Half of the larval Defensins were induced by mating or feeding in adults, and all three of the induced Defensins were located downstream of canonical binding sites for an Imd transcription factor involved in the highly conserved NF-κB signaling that regulates induction of AMPs. The larval and non-larval Defensins were located in distinct genomic regions, and the amino acid sequences of the larval Defensins formed a monophyletic clade. There were more amino acid substitutions across non-larval Defensins, with multiple genes losing a highly conserved furin cleavage site thought to be required for the removal of the amino terminus from the mature Defensin domain. However, larval Defensins had a higher proportion of radical amino acid substitutions, altering amino acid size and polarity. Our results reveal insights into the developmental stage-specific regulation of AMPs, and they suggest different regulatory regimes impose unique selection pressures on AMPs, possibly as a result of variation in exposure to microbial communities across development.

Relationship between acaricide resistance and acetylcholinesterase gene polymorphisms in the cattle tick Rhipicephalus microplus.

In this study, we aimed to develop a comprehensive methodology for identifying amino acid polymorphisms in acetylcholinesterase transcript 2 (AChE2) in acaricide-resistant Rhipicephalus microplus ticks. This included assessing AChE2 expression levels through qPCR and conducting 3D modeling to evaluate the interaction between acaricides and AChE2 using docking techniques. The study produced significant results, demonstrating that acaricide-resistant R. microplus ticks exhibit significantly higher levels of AChE expression than susceptible reference ticks. In terms of amino acid sequence, we identified 9 radical amino acid substitutions in AChE2 from acaricide-resistant ticks, when compared to the gene sequence of the susceptible reference strain. To further understand the implications of these substitutions, we utilized 3D acaricide-AChE2 docking modeling to examine the interaction between the acaricide and the AChE2 catalytic site. Our models suggest that these amino acid polymorphisms alter the configuration of the binding pocket, thereby contributing to differences in acaricide interactions and ultimately providing insights into the acaricide-resistance phenomenon in R. microplus.

Single Amino Acid Polymorphisms in the Fasciola hepatica Carboxylesterase Type B Gene and Their Potential Role in Anthelmintic Resistance.

The expression of the Fasciola hepatica carboxylesterase type B (CestB) gene is known to be induced upon exposure to the anthelmintic triclabendazole (TCBZ), leading to a substantial rise in enzyme-specific activity. Furthermore, the nucleotide sequence of the CestB gene displays variations that can potentially result in radical amino acid substitutions at the ligand binding site. These substitutions hold the potential to impact both the ligand-protein interaction and the catalytic properties of the enzyme. Thus, the objective of our study was to identify novel CestB polymorphisms in TCBZ-resistant parasites and field isolates obtained from a highly endemic region in Central Mexico. Additionally, we aimed to assess these amino acid polymorphisms using 3D modeling against the metabolically oxidized form of the anthelmintic TCBZSOX. Our goal was to observe the formation of TCBZSOX-specific binding pockets that might provide insights into the role of CestB in the mechanism of anthelmintic resistance. We identified polymorphisms in TCBZ-resistant parasites that exhibited three radical amino acid substitutions at positions 147, 215, and 263. These substitutions resulted in the formation of a TCBZSOX-affinity pocket with the potential to bind the anthelmintic drug. Furthermore, our 3D modeling analysis revealed that these amino acid substitutions also influenced the configuration of the CestB catalytic site, leading to alterations in the enzyme's interaction with chromogenic carboxylic ester substrates and potentially affecting its catalytic properties. However, it is important to note that the TCBZSOX-binding pocket, while significant for drug binding, was located separate from the enzyme's catalytic site, rendering enzymatic hydrolysis of TCBZSOX impossible. Nonetheless, the observed increased affinity for the anthelmintic may provide an explanation for a drug sequestration type of anthelmintic resistance. These findings lay the groundwork for the future development of a molecular diagnostic tool to identify anthelmintic resistance in F. hepatica.

A Single Nucleotide Polymorphism Translates into a Radical Amino Acid Substitution at the Ligand-Binding Site in Fasciola hepatica Carboxylesterase B.

Fasciola hepatica anthelmintic resistance may be associated with the catalytic activity of xenobiotic metabolizing enzymes. The gene expression of one of these enzymes, identified as carboxylesterase B (CestB), was previously described as inducible in adult parasites under anthelmintic treatment and exhibited a single nucleotide polymorphism at position 643 that translates into a radical amino acid substitution at position 215 from Glutamic acid to Lysine. Alphafold 3D models of both allelic sequences exhibited a significant affinity pocket rearrangement and different ligand-docking modeling results. Further bioinformatics analysis confirmed that the radical amino acid substitution is located at the ligand affinity site of the enzyme, affecting its affinity to serine hydrolase inhibitors and preferences for ester ligands. A field genotyping survey from parasite samples obtained from two developmental stages isolated from different host species from Argentina and Mexico exhibited a 37% allele distribution for 215E and a 29% allele distribution for 215K as well as a 34% E/K heterozygous distribution. No linkage to host species or geographic origin was found in any of the allele variants.

Overexpression of glutamate-gated chloride channel in the integument is mainly responsible for emamectin benzoate resistance in the western flower thrips Frankliniella occidentalis.

Western flower thrips Frankliniella occidentalis is a serious polyphagous pest worldwide. In this study, we investigated the potential mechanisms of resistance including enhanced metabolism and target site insensitivity in an emamectin benzoate (EB)-resistant (EB-R) strain. The EB-R strain of F. occidentalis showed 356-fold increased resistance compared to a susceptible RDA strain. Analysis of cross-resistance to four other insecticides confirmed that EB resistance is highly specific to the contact toxicity of EB. Synergistic bioassay and quantitative PCR of cytochrome P450 monooxygenase (CYP) genes revealed that three overexpressed Cyps were likely involved in resistance. Among three putative glutamate-gated chloride channel (GluCl) genes identified, FoGluClc showed four radical amino acid substitutions and 3.8-fold and 31-fold transcription level in the head and integument in the EB-R strain when compared to the RDA strain. Backcrossing analysis and RNA interference confirmed that both amino acid substitution and overexpression of FoGluClc are responsible for EB resistance. In situ hybridization revealed that FoGluClc is mainly distributed in the integument in the EB-R strain. Cross-comparison of known genomes and transcriptomes of thrips species revealed that FoGluClc is unique to the Frankliniella genus. While mutations and overexpression of FoGluClc play major roles in EB resistance, the overexpressed Cyps are partially involved as metabolic factors. Higher expression of FoGluClc in the integument may suggest its role in the first-line defense against EB in the EB-R strain. Unique distribution of FoGluClc in the Frankliniella genus but not in other thrips species further suggests that FoGluClc may be a surplus channel not having an essential endogenous function and is thus recruited as a defense barrier against xenobiotics. © 2022 Society of Chemical Industry.

Evolutionary changes between pre- and post-vaccine South African group A G2P[4] rotavirus strains, 2003-2017.

The transient upsurge of G2P[4] group A rotavirus (RVA) after Rotarix vaccine introduction in several countries has been a matter of concern. To gain insight into the diversity and evolution of G2P[4] strains in South Africa pre- and post-RVA vaccination introduction, whole-genome sequencing was performed for RVA positive faecal specimens collected between 2003 and 2017 and samples previously sequenced were obtained from GenBank (n=103; 56 pre- and 47 post-vaccine). Pre-vaccine G2 sequences predominantly clustered within sub-lineage IVa-1. In contrast, post-vaccine G2 sequences clustered mainly within sub-lineage IVa-3, whereby a radical amino acid (AA) substitution, S15F, was observed between the two sub-lineages. Pre-vaccine P[4] sequences predominantly segregated within sub-lineage IVa while post-vaccine sequences clustered mostly within sub-lineage IVb, with a radical AA substitution R162G. Both S15F and R162G occurred outside recognised antigenic sites. The AA residue at position 15 is found within the signal sequence domain of Viral Protein 7 (VP7) involved in translocation of VP7 into endoplasmic reticulum during infection process. The 162 AA residue lies within the hemagglutination domain of Viral Protein 4 (VP4) engaged in interaction with sialic acid-containing structure during attachment to the target cell. Free energy change analysis on VP7 indicated accumulation of stable point mutations in both antigenic and non-antigenic regions. The segregation of South African G2P[4] strains into pre- and post-vaccination sub-lineages is likely due to erstwhile hypothesized stepwise lineage/sub-lineage evolution of G2P[4] strains rather than RVA vaccine introduction. Our findings reinforce the need for continuous whole-genome RVA surveillance and investigation of contribution of AA substitutions in understanding the dynamic G2P[4] epidemiology.

Meta-analysis of tumor- and T cell-intrinsic mechanisms of sensitization to checkpoint inhibition

SummaryCheckpoint inhibitors (CPIs) augment adaptive immunity. Systematic pan-tumor analyses may reveal the relative importance of tumor-cell-intrinsic and microenvironmental features underpinning CPI sensitization. Here, we collated whole-exome and transcriptomic data for >1,000 CPI-treated patients across seven tumor types, utilizing standardized bioinformatics workflows and clinical outcome criteria to validate multivariable predictors of CPI sensitization. Clonal tumor mutation burden (TMB) was the strongest predictor of CPI response, followed by total TMB and CXCL9 expression. Subclonal TMB, somatic copy alteration burden, and histocompatibility leukocyte antigen (HLA) evolutionary divergence failed to attain pan-cancer significance. Dinucleotide variants were identified as a source of immunogenic epitopes associated with radical amino acid substitutions and enhanced peptide hydrophobicity/immunogenicity. Copy-number analysis revealed two additional determinants of CPI outcome supported by prior functional evidence: 9q34 (TRAF2) loss associated with response and CCND1 amplification associated with resistance. Finally, single-cell RNA sequencing (RNA-seq) of clonal neoantigen-reactive CD8 tumor-infiltrating lymphocytes (TILs), combined with bulk RNA-seq analysis of CPI-responding tumors, identified CCR5 and CXCL13 as T-cell-intrinsic markers of CPI sensitivity.

Differences in the sucrose synthase gene SUS1 expression pattern between Solanum lycopersicum and wild tomato species

Plant sucrose synthase (SUS) reversibly hydrolyzes sucrose, incorporating sucrose into different metabolic pathways. SUS1 homologous genes were characterized in 12 accessions of cultivated tomato Solanum lycopersicum and wild tomato species (Solanum section Lycopersicon). In total, 34 non-synonymous SNPs were found; among them, 17 caused radical amino acid substitutions, mostly in the wild species. The SUS1-based phylogenetic tree topology was represented by two clusters of red-fruited and green-fruited species, where S. pennellii and S. habrochaites showed the highest divergence. SUS1 transcription was detected in source and sink organs of one cultivated and four wild tomato species; it was the highest in buds and flowers, but lower in mature green fruit of the final size. In leaves and roots, the level of SUS1 expression corresponded to species evolutionary history, whereas no correlation was found between SUS1 levels and fruit size and sugar composition. Our findings suggest that SUS1-mediated reversible sucrose hydrolysis is important for maintaining the balance between sucrose and its monomers not only in fruits, but also in the other tomato organs.

Adaptive evolution of the ACSL gene family in Carnivora.

Carnivores exhibit various fat contents and energy reserves to adapt to their environments. However, the molecular mechanisms underlying lipid metabolic differences among carnivores have not been well explored. Long-chain acyl-CoA synthetases (ACSLs) catalyze the initial step in lipid metabolism by activating fatty acids (FAs), and they drive acyl-CoAs toward anabolic lipid synthesis or catabolic β-oxidation. We identified the sequences of the genes of the ACSL family (ACSL1, ACSL3, ACSL4, ACSL5 and ACSL6) in the sable (Martes zibellina) via transcriptome sequencing. The ACSL gene sequences of 13 other carnivores were obtained from NCBI. Phylogenetic results showed that unlike the widely accepted carnivore phylogeny, Canidae and Felidae tend to group together based on ACSL4 and ACSL6. The evolutionary analyses identified a series of positively selected amino acid residues in ACSL1, ACSL4 and ACSL5. Two radical amino acid substitutions detected in sable suggested potential insights into the molecular mechanism underlying the relatively low fat content in this animal. This is the first study to investigate the molecular mechanisms underlying the adaptive evolution of fat metabolism in carnivores. Overall, the ACSL genes were under different evolutionary forces in carnivores, and some genes have undergone adaptive evolution in lipid metabolism.

Comprehensive phylogenomic analysis reveals a novel cluster of simian endogenous retroviral sequences in Colobinae monkeys.

Simian retrovirus (SRV) is a type-D betaretrovirus infectious to the Old World monkeys causing a variety of symptoms. SRVs are also present in the Old World monkey genomes as endogenous forms, which are referred to as Simian endogenous retroviruses (SERVs). Although many SERV sequences have been identified in Cercopithecinae genomes, with potential of encoding all functional genes, the distribution of SERVs in genomes and evolutionary relationship between exogeneous SRVs and SERVs remains unclear. In this study, we comprehensively investigated seven draft genome sequences of the Old World monkeys, and identified a novel cluster of SERVs in the two Rhinopithecus (R. roxellana and R. bieti) genomes, which belong to the Colobinae subfamily. The Rhinopithecus genomes harbored higher copy numbers of SERVs than the Cercopithecinae genomes. A reconstructed phylogenetic tree showed that the Colobinae SERVs formed a distinct cluster from SRVs and Cercopithecinae SERVs, and more closely related to exogenous SRVs than Cercopithecinae SERVs. Three radical amino acid substitutions specific to Cercopithecinae SERVs, which potentially affect the infectious ability of SERVs, were also identified in the proviral envelope protein. In addition, we found many integration events of SERVs were genus- or species-specific, suggesting the recent activity of SERVs in the Old World monkey genomes. The results suggest that SERVs in Cercopithecinae and Colobinae monkeys were endogenized after the divergence of subfamilies and do not transmit across subfamilies. Our findings also support the hypothesis that Colobinae SERVs are direct ancestors of SRV-6, which has a different origin from the other exogenous SRVs. These findings shed novel insight into the evolutionary history of SERVs, and may help to understand the process of endogenization of SRVs.

Substitution of Aromatic Residues with Polar Residues in the Active Site Pocket of epi-Isozizaene Synthase Leads to the Generation of New Cyclic Sesquiterpenes.

The sesquiterpene cyclase epi-isozizaene synthase (EIZS) catalyzes the cyclization of farnesyl diphosphate to form the tricyclic hydrocarbon precursor of the antibiotic albaflavenone. The hydrophobic active site pocket of EIZS serves as a template as it binds and chaperones the flexible substrate and carbocation intermediates through the conformations required for a multistep reaction sequence. We previously demonstrated that the substitution of hydrophobic residues with other hydrophobic residues remolds the template and expands product chemodiversity [Li, R., Chou, W. K. W., Himmelberger, J. A., Litwin, K. M., Harris, G. G., Cane, D. E., and Christianson, D. W. (2014) Biochemistry 53, 1155-1168]. Here, we show that the substitution of hydrophobic residues-specifically, Y69, F95, F96, and W203-with polar side chains also yields functional enzyme catalysts that expand product chemodiversity. Fourteen new EIZS mutants are reported that generate product arrays in which eight new sesquiterpene products have been identified. Of note, some mutants generate acyclic and cyclic hydroxylated products, suggesting that the introduction of polarity in the hydrophobic pocket facilitates the binding of water capable of quenching carbocation intermediates. Furthermore, the substitution of polar residues for F96 yields high-fidelity sesquisabinene synthases. Crystal structures of selected mutants reveal that residues defining the three-dimensional contour of the hydrophobic pocket can be substituted without triggering significant structural changes elsewhere in the active site. Thus, more radical nonpolar-polar amino acid substitutions should be considered when terpenoid cyclase active sites are remolded by mutagenesis with the goal of exploring and expanding product chemodiversity.

Patterns of Gene Conversion in Duplicated Yeast Histones Suggest Strong Selection on a Coadapted Macromolecular Complex.

We find evidence for interlocus gene conversion in five duplicated histone genes from six yeast species. The sequences of these duplicated genes, surviving from the ancient genome duplication, show phylogenetic patterns inconsistent with the well-resolved orthology relationships inferred from a likelihood model of gene loss after the genome duplication. Instead, these paralogous genes are more closely related to each other than any is to its nearest ortholog. In addition to simulations supporting gene conversion, we also present evidence for elevated rates of radical amino acid substitutions along the branches implicated in the conversion events. As these patterns are similar to those seen in ribosomal proteins that have undergone gene conversion, we speculate that in cases where duplicated genes code for proteins that are a part of tightly interacting complexes, selection may favor the fixation of gene conversion events in order to maintain high protein identities between duplicated copies.

Development of a Homologous Expression System for and Systematic Site-Directed Mutagenesis Analysis of Thurincin H, a Bacteriocin Produced by Bacillus thuringiensis SF361

Thurincin H is an antimicrobial peptide produced by Bacillus thuringiensis SF361. With a helical back bone, the 31 amino acids of thurincin H form a hairpin structure maintained by four pairs of very unique sulfur-to-α-carbon thioether bonds. The production of thurincin H depends on a putative gene cluster containing 10 open reading frames. The gene cluster includes three tandem structural genes (thnA1, thnA2, and thnA3) encoding three identical 40-amino-acid thurincin H prepeptides and seven other genes putatively responsible for prepeptide processing, regulation, modification, exportation, and self-immunity. A homologous thurincin H expression system was developed by transforming a thurincin H-deficient host with a novel expression vector, pGW133. The host, designated B. thuringiensis SF361 ΔthnA1 ΔthnA2 ΔthnA3, was constructed by deletion of the three tandem structural genes from the chromosome of the natural thurincin H producer. The thurincin H expression vector pGW133 was constructed by cloning the thurincin H native promoter, thnA1, and a Cry protein terminator into the Escherichia coli-B. thuringiensis shuttle vector pHT315. Thirty-three different pGW133 variants, each containing a different point mutation in the thnA1 gene, were generated and separately transformed into B. thuringiensis SF361 ΔthnA1 ΔthnA2 ΔthnA3. Those site-directed mutants contained either a single radical or conservative amino acid substitution on the thioether linkage-forming positions or a radical substitution on all other nonalanine amino acids. The bacteriocin activities of B. thuringiensis SF361 ΔthnA1 ΔthnA2 ΔthnA3 carrying different pGW133 variants against three different indicator strains were subsequently compared.

Functional Divergence for Every Paralog

Because genes can be constrained by selection at more than one phenotypic level, the relaxation of constraints following gene duplication allows for functional divergence (FD) along multiple phenotypic axes. Many studies have generated individual measures of FD, but the profile of FD between paralogs across levels of phenotypic space remains largely uncharted. We evaluate paralog pairs that originated via the yeast whole-genome duplication (ohnolog pairs) at three distinct phenotypic levels (properties of proteins, gene expression, and overall organismal growth) using eight complementary measures of FD (protein: evolutionary rates, radical amino acid substitutions, and domain architecture; gene expression: expression differences in a single species and condition, across species in a single condition, and in a single species across conditions; and organismal: genetic interaction profiles and growth profiles in multiple conditions). We find that the majority of ohnolog pairs show FD by multiple phenotypic measures. Within each phenotypic level, measures of FD are strongly correlated but are generally weakly correlated between levels, suggesting that functional constraints exerted on genes from distinct phenotypic levels are largely decoupled. Our results suggest that redundancy is a rare functional fate for retained paralogs and that FD cannot be fully captured by measures at any single phenotypic level.

Selective processes and adaptive evolution of the cytochrome b gene in salamanders of the genus Salamandrella

Sequence analysis of the cytochrome b gene fragment in the salamanders of the genus Salamandrella, Siberian salamander and Schrenk salamander was performed with the purpose to elucidate the effect of natural selection on the evolution of mitochondrial DNA (mtDNA) in these species. It was demonstrated that despite of notable influence of negative selection (expressed as very low dN/dS values), speciation and intraspecific divergence in salamanders was accompanied by the appearance of radical amino acid substitutions, caused by the influence of positive (directional) selection. To examine the evolutionary pattern of synonymous mtDNA sites, distribution of conservative and non-conservative substitutions was analyzed. The rates of conservative and non-conservative substitutions were nearly equal, pointing to neutrality of mutation process at synonymous mtDNA sites of salamanders. Analysis of conservative and non-conservative synonymous substitution distributions in different parts of phylogenetic trees showed that the differences between the synonymous groups compared were statistically significant only in one phylogenetic group of Siberian salamander (haplogroup C) (P = 0.02). In the group of single substitutions, located at terminal phylogenetic branches of Siberian salamanders from group C, increased rate of conservative substitutions was observed. Based on these findings, it was suggested that selective processes could have an influence on the formation of the synonymous substitution profile in the Siberian salamander mtDNA fragment examined.

Reduced Selective Constraint in Endosymbionts: Elevation in Radical Amino Acid Replacements Occurs Genome-Wide

As predicted by the nearly neutral model of evolution, numerous studies have shown that reduced Ne accelerates the accumulation of slightly deleterious changes under genetic drift. While such studies have mostly focused on eukaryotes, bacteria also offer excellent models to explore the effects of Ne. Most notably, the genomes of host-dependent bacteria with small Ne show signatures of genetic drift, including elevated Ka/Ks. Here, I explore the utility of an alternative measure of selective constraint: the per-site rate of radical and conservative amino acid substitutions (Dr/Dc). I test the hypothesis that purifying selection against radical amino acid changes is less effective in two insect endosymbiont groups (Blochmannia of ants and Buchnera of aphids), compared to related gamma-Proteobacteria. Genome comparisons demonstrate a significant elevation in Dr/Dc in endosymbionts that affects the majority (66–79%) of shared orthologs examined. The elevation of Dr/Dc in endosymbionts affects all functional categories examined. Simulations indicate that Dr/Dc estimates are sensitive to codon frequencies and mutational parameters; however, estimation biases occur in the opposite direction as the patterns observed in genome comparisons, thereby making the inference of elevated Dr/Dc more conservative. Increased Dr/Dc and other signatures of genome degradation in endosymbionts are consistent with strong effects of genetic drift in their small populations, as well as linkage to selected sites in these asexual bacteria. While relaxed selection against radical substitutions may contribute, genome-wide processes such as genetic drift and linkage best explain the pervasive elevation in Dr/Dc across diverse functional categories that include basic cellular processes. Although the current study focuses on a few bacterial lineages, it suggests Dr/Dc is a useful gauge of selective constraint and may provide a valuable alternative to Ka/Ks when high sequence divergences preclude estimates of Ks. Broader application of Dr/Dc will benefit from approaches less prone to estimation biases.

Adaptive intraspecific divergence: An example using the animal cytochrome b gene

The topologies of phylogenetic trees characterized by a high level of intraspecific divergence between the phylogenetic DNA groups (clades) are often explained in terms of the theory of Pleistocene refugia. To elucidate the issue of the adaptive role of intraspecific divergence, the changes in the physicochemical properties of amino acids in the course ofcladogenesis (MM01 model of the TreeSAAP 3.2 package) were analyzed in this work using as an example the nucleotide sequences of the cytochrome b gene in some species of northern animals (lemmings, redbacked voles, chipmunks, flying squirrels, ermines). It was shown that the process of intraspecific divergence was rarely accompanied by radical amino acid substitutions in cytochrome b caused by adaptation (directional selection). In connection with this, the hypothesis is discussed according to which the adaptive variants formed in the species at the peak of cold were lost with climatic warming due to the drift or selection against individuals adapted to cold.

Adaptive Introgression of Anticoagulant Rodent Poison Resistance by Hybridization between Old World Mice

Polymorphisms in the vitamin K 2,3-epoxide reductase subcomponent 1 (vkorc1) of house mice (Mus musculus domesticus) can cause resistance to anticoagulant rodenticides such as warfarin [1-3]. Here we show that resistant house mice can also originate from selection on vkorc1 polymorphisms acquired from the Algerian mouse (M. spretus) through introgressive hybridization. We report on a polymorphic introgressed genomic region in European M. m. domesticus that stems from M. spretus, spans >10 Mb on chromosome 7, and includes the molecular target of anticoagulants vkorc1 [1-4]. We show that in the laboratory, the homozygous complete vkorc1 allele of M. spretus confers resistance when introgressed into M. m. domesticus. Consistent with selection on the introgressed allele after the introduction of rodenticides in the 1950s, we found signatures of selection in patterns of variation in M. m. domesticus. Furthermore, we detected adaptive protein evolution of vkorc1 in M. spretus (Ka/Ks = 1.54-1.93) resulting in radical amino acid substitutions that apparently cause anticoagulant tolerance in M. spretus as a pleiotropic effect. Thus, positive selection produced an adaptive, divergent, and pleiotropic vkorc1 allele in the donor species, M. spretus, which crossed a species barrier and produced an adaptive polymorphic trait in the recipient species, M. m. domesticus.

The linker domain of basal transcription factor TFIIB controls distinct recruitment and transcription stimulation functions

RNA polymerases (RNAPs) require basal transcription factors to assist them during transcription initiation. One of these factors, TFIIB, combines promoter recognition, recruitment of RNAP, promoter melting, start site selection and various post-initiation functions. The ability of 381 site-directed mutants in the TFIIB ‘linker domain’ to stimulate abortive transcription was systematically quantitated using promoter-independent dinucleotide extension assays. The results revealed two distinct clusters (mjTFIIB E78-R80 and mjTFIIB R90-G94, respectively) that were particularly sensitive to substitutions. In contrast, a short sequence (mjTFIIB A81-K89) between these two clusters tolerated radical single amino acid substitutions; short deletions in that region even caused a marked increase in the ability of TFIIB to stimulate abortive transcription (‘superstimulation’). The superstimulating activity did, however, not correlate with increased recruitment of the TFIIB/RNAP complex because substitutions in a particular residue (mjTFIIB K87) increased recruitment by more than 5-fold without affecting the rate of abortive transcript stimulation. Our work demonstrates that highly localized changes within the TFIIB linker have profound, yet surprisingly disconnected, effects on RNAP recruitment, TFIIB/RNAP complex stability and the rate of transcription initiation. The identification of superstimulating TFIIB variants reveals the existence of a previously unknown rate-limiting step acting on the earliest stages of gene expression.