Codon Sites Research Articles

Mitochondrial genomes of twelve species of hyperdiverse Trigonopterus weevils.

Mitochondrial genomes of twelve species of Trigonopterus weevils are presented, ten of them complete. We describe their gene order and molecular features and test their potential for reconstructing the phylogeny of this hyperdiverse genus comprising > 1,000 species. The complete mitochondrial genomes examined herein ranged from 16,501 bp to 21,007 bp in length, with an average AT content of 64.2% to 69.7%. Composition frequencies and skews were generally lower across species for atp6, cox1-3, and cob genes, while atp8 and genes coded on the minus strand showed much higher divergence at both nucleotide and amino acid levels. Most variation within genes was found at the codon level with high variation at third codon sites across species, and with lesser degree at the coding strand level. Two large non-coding regions were found, CR1 (between rrnS and trnI genes) and CR2 (between trnI and trnQ), but both with large variability in length; this peculiar structure of the non-coding region may be a derived character of Curculionoidea. The nad1 and cob genes exhibited an unusually high interspecific length variation of up to 24 bp near the 3′ end. This pattern was probably caused by a single evolutionary event since both genes are only separated by trnS2 and length variation is extremely rare in mitochondrial protein coding genes. We inferred phylogenetic trees using protein coding gene sequences implementing both maximum likelihood and Bayesian approaches, each for both nucleotide and amino acid sequences. While some clades could be retrieved from all reconstructions with high confidence, there were also a number of differences and relatively low support for some basal nodes. The best partition scheme of the 13 protein coding sequences obtained by IQTREE suggested that phylogenetic signal is more accurate by splitting sequence variation at the codon site level as well as coding strand, rather than at the gene level. This result corroborated the different patterns found in Trigonopterus regarding to A+T frequencies and AT and GC skews that also greatly diverge at the codon site and coding strand levels.

Variation around the dominant viral genome sequence contributes to viral load and outcome in patients with Ebola virus disease

BackgroundViral load is a major contributor to outcome in patients with Ebola virus disease (EVD), with high values leading to a fatal outcome. Evidence from the 2013–2016 Ebola virus (EBOV) outbreak indicated that different genotypes of the virus can have different phenotypes in patients. Additionally, due to the error-prone nature of viral RNA synthesis in an individual patient, the EBOV genome exists around a dominant viral genome sequence. The minor variants within a patient may contribute to the overall phenotype in terms of viral protein function. To investigate the effects of these minor variants, blood samples from patients with acute EVD were deeply sequenced.ResultsWe examine the minor variant frequency between patients with acute EVD who survived infection with those who died. Non-synonymous differences in viral proteins were identified that have implications for viral protein function. The greatest frequency of substitution was identified at three codon sites in the L gene—which encodes the viral RNA-dependent RNA polymerase (RdRp). Recapitulating this in an assay for virus replication, these substitutions result in aberrant viral RNA synthesis and correlate with patient outcome.ConclusionsTogether, these findings support the notion that in patients who survived EVD, in some cases, the genetic variability of the virus resulted in deleterious mutations that affected viral protein function, leading to reduced viral load. Such mutations may also lead to persistent strains of the virus and be associated with recrudescent infections.

Differences of cytotoxic T-lymphocyte pressure and divergent evolution of several CRF07_BC clusters circulating in men who have sex with men in China

The HIV-1 CRF07_BC strain is now the major recombinant form in China among the population of men who have sex with men (MSM), and has critically contributed to the HIV-1 epidemic in recent years. The phylodynamic and virological differences among CRF07_BC clusters circulating in MSM, and the factors that could be driving their evolution, remains unclear. Using the HIV-1 CRF07_BC strains obtained from the Los Alamos HIV database, we undertook a large-scale phylogenetic analysis using the maximum likelihood method of partial gag, pol, and env segments to infer their evolutionary relationships. The demographic histories of clusters were determined using the Bayesian Markov chain Monte Carlo (MCMC) method. For four pol clusters we analyzed the non-synonymous (dN) to synonymous (dS) substitution rates and performed site to site analysis to identify positive selection sites and cytotoxic T-lymphocyte (CTL) escape mutation positions. MSM was found to be the predominant risk factor for all four of the CRF07_BC epidemic pol segment clusters with the largest number of infections. Two of those clusters had higher growth in the effective number of infections, and two clusters had slower growth. Analysis of all four clusters showed no significant differences in the mean substitution rates and dN/dS selection pressure ratios. However, a site to site codon analysis found thirteen significant positive selection sites. Ten of these sites are CTL escape mutations. The two clusters with the higher growth in infections had seven and eight pol segment CTL escape mutation sites respectively, while the two with slower growth had only one or two. Our findings demonstrated differences in the CTL escape mutation and divergent evolution of several CRF07_BC clusters circulating among men who have sex with men in China.

Loss-of-function variants in FSIP1 identified by targeted sequencing are associated with one particular subtype of mucosal melanoma.

Mucosal melanoma is a tumor caused by the malignant transformation of pigment-producing cells and can arise from any mucosal tissue where melanocytes are present. Due to its rarity, the mucosal melanoma subtype is poorly described, and its genetic characteristics are infrequently studied. The discovery or confirmation of new mucosal melanoma susceptibility genes will provide important insights for the study of its pathogenesis. We performed deep targeted sequencing of 100 previously reported melanoma-related genes in 39 mucosal melanoma samples and a gene-level loss-of-function (LOF) variant enrichment analysis for mucosal melanoma from different incidence sites. We detected 7,589 variants in these samples, and 484 were LOF variants (gain or loss of a stop codon, missense, and splice site). Four different gene-level enrichment analyses revealed that FSIP1 (fibrous sheath interacting protein 1) is a susceptibility gene for oral mucosal melanoma (OR=0.33, PChi=4.05×10-2, Pburden=3.06×10-2, Pskat=3.01×10-2, Pskato=3.01×10-2), whereas the different methods did not detect a significant susceptibility gene for the other subtypes. In our study, a susceptibility gene for oral mucosal melanoma was confirmed in a Chinese Han population, and these findings contribute to a better genetic understanding of mucosal melanoma of different subtypes.

Evolutionary analysis of genes coding for Cysteine-RIch Secretory Proteins (CRISPs) in mammals

BackgroundCysteine-RIch Secretory Proteins (CRISP) are expressed in the reproductive tract of mammalian males and are involved in fertilization and related processes. Due to their important role in sperm performance and sperm-egg interaction, these genes are likely to be exposed to strong selective pressures, including postcopulatory sexual selection and/or male-female coevolution. We here perform a comparative evolutionary analysis of Crisp genes in mammals. Currently, the nomenclature of CRISP genes is confusing, as a consequence of discrepancies between assignments of orthologs, particularly due to numbering of CRISP genes. This may generate problems when performing comparative evolutionary analyses of mammalian clades and species. To avoid such problems, we first carried out a study of possible orthologous relationships and putative origins of the known CRISP gene sequences. Furthermore, and with the aim to facilitate analyses, we here propose a different nomenclature for CRISP genes (EVAC1–4, “EVolutionarily-analyzed CRISP”) to be used in an evolutionary context.ResultsWe found differing selective pressures among Crisp genes. CRISP1/4 (EVAC1) and CRISP2 (EVAC2) orthologs are found across eutherian mammals and seem to be conserved in general, but show signs of positive selection in primate CRISP1/4 (EVAC1). Rodent Crisp1 (Evac3a) seems to evolve under a comparatively more relaxed constraint with positive selection on codon sites. Finally, murine Crisp3 (Evac4), which appears to be specific to the genus Mus, shows signs of possible positive selection. We further provide evidence for sexual selection on the sequence of one of these genes (Crisp1/4) that, unlike others, is thought to be exclusively expressed in male reproductive tissues.ConclusionsWe found differing selective pressures among CRISP genes and sexual selection as a contributing factor in CRISP1/4 gene sequence evolution. Our evolutionary analysis of this unique set of genes contributes to a better understanding of Crisp function in particular and the influence of sexual selection on reproductive mechanisms in general.

Adaptive Molecular Evolution of AKT3 Gene for Positive Diversifying Selection in Mammals.

The V-Akt Murine Thymoma Viral Oncogene Homolog 3 (AKT3) gene is of the serine/threonine-protein kinase family and influences the production of milk fats and cholesterol by acting on the sterol administrative area restricting protein (SREBP). The AKT3 gene is highly preserved in animals, and during lactation in cattle, its expression increases. The AKT3 gene is expressed in the digestive system, mammary gland, and immune cells. A phylogenetic investigation was performed to clarify the evolutionary role of AKT3, by maximum probability. The AKT3 gene sequence data of various mammalian species was evident even with animals undergoing breeding selection. From 39 mammalian species studied, there was a signal of positive diversifying selection with Hominidae at 13Q, 16G, 23R, 24P, 121P, 294K, 327V, 376L, 397K, 445T, and 471F among other codon sites of the AKT3 gene. These sites were codes for amino acids such as arginine, proline, lysine, and leucine indicating major roles for the function of immunological proteins, and in particular, the study highlighted the importance of changes in gene expression of AKT3 on immunity.

Analysis of sequence diversity and selection pressure in HIV-1 clade C gp41 from India.

Evaluation of viral diversity is critical for the rational design of treatment modalities against Human immunodeficiency virus (HIV). Predominated by HIV-1 clade C (HIV-1C), the epidemic in India represents the third largest population infected with HIV-1 globally. Glycoprotein 41 (gp41) is critical for viral replication and is a target for the design of therapeutic strategies. However, documentation of viral diversity of gp41 gene in infected individuals from India remains limited. Present study employed high throughput sequencing to examine variation in gp41 amplicons generated from blood derived viruses in 24 HIV-1C infected individuals from Mumbai, India. Sequence diversity profiles were documented in different functional domains of gp41. Furthermore, through a meta-analysis approach, all reported gp41 sequences from India (N = 70) were compared with those from South Africa (N = 126), country with the largest HIV epidemic globally, also predominated by HIV-1C. A total of 44 positions displayed statistically significant differential (p < 0.05) Shannon entropy in the two regions. This comparison also identified 11 codon sites undergoing distinct selection, 8 of which remained differentially selected in an extended comparison of data from Asia (N = 137) and Africa(N = 383). Assessment of correlated mutation networks associated with differentially selected residues revealed common as well as distinct interaction networks. Furthermore, codon usage analysis revealed 17 differentially selected codons (Mann-Whitney test, p < 0.001) in Asia and Africa. Dissimilar trends in GC content across codon positions were also observed. In depth understanding of these divergent evolutionary signatures through extended analysis with larger data-sets would assist development of effective interventions being considered for HIV-1C.

Genetic analysis of HA1 domain of influenza A/H3N2 viruses isolated in Kenya during the 2007–2013 seasons reveal significant divergence from WHO-recommended vaccine strains

BackgroundInfluenza viruses evolve rapidly and cause regular seasonal epidemics in humans challenging effective vaccination. The virus surface HA glycoprotein is the primary target for the host immune response. Here, we investigated the vaccine efficacy and evolution patterns of human influenza A/H3N2 viruses that circulated in Kenyan in the period before and after the 2009 A/H1N1 pandemic, targeting the HA1 domain. Materials and methodsA hundred and fifteen HA sequences of Kenyan virus viruses were analyzed relative to the corresponding WHO vaccine reference strains using bioinformatics approaches. ResultsOur analyses revealed varied amino acid substitutions at all the five antigenic sites (A–E) of the HA1 domain, with a majority the changes occurring at sites A and B. The Kenyan A/H3N2 viruses isolated during 2007/2008 seasons belonged to A/Brisbane/10/2007-like viruses lineage, while those circulating in 2009–2012 belonged to the lineage of A/Victoria/361/2011-like viruses. The 2013 viruses clustered in clade 3C.3 of the A/Samara/73/2013-like viruses. The mean evolutionary rate of the A/H3N2 viruses analyzed in the study was at 4.17×10−3 (95% HPD=3.09×10−3–5.31×10−3) nucleotide substitutions per site per year, whereas the TMRCA was estimated at 11.18 (95% HPD=9.00–14.12) years ago from 2013. The prediction of vaccine efficacy revealed modest vaccine efficaciousness during 2008, and 2010 influenza seasons, whilst sub-optimal effectiveness was registered in 2007, 2009, 2012 and 2013. Further, the overall selective pressure acting on the HA1 domain was estimated at 0.56 (ω<1), suggesting that a majority of codon sites in the HA1 epitopes were evolving under purifying selection. ConclusionsGenerally, our results highlight the genetic plasticity of A/H3N2 viruses and reveal considerable disparity in vaccine efficaciousness against the A/H3N2 viruses that circulated in Kenya, specifically during 2007, 2009, 2012, and 2013 influenza seasons. Our findings underscore the importance and need for consistent surveillance and molecular characterization of influenza viruses, to inform decision making and enhance early of detection of strains with epidemic/pandemic potential as well as benefit in guiding decisions regarding the appropriate annual influenza vaccine formulations.

Structural and functional insights into a novel two-component endolysin encoded by a single gene in Enterococcus faecalis phage.

Using bacteriophage-derived endolysins as an alternative strategy for fighting drug-resistant bacteria has recently been garnering renewed interest. However, their application is still hindered by their narrow spectra of activity. In our previous work, we demonstrated that the endolysin LysIME-EF1 possesses efficient bactericidal activity against multiple strains of Enterococcus faecalis (E. faecalis). Herein, we observed an 8 kDa fragment and hypothesized that it contributes to LysIME-EF1 lytic activity. To examine our hypothesis, we determined the structure of LysIME-EF1 at 1.75 Å resolution. LysIME-EF1 exhibits a unique architecture in which one full-length LysIME-EF1 forms a tetramer with three additional C-terminal cell-wall binding domains (CBDs) that correspond to the abovementioned 8 kDa fragment. Furthermore, we identified an internal ribosomal binding site (RBS) and alternative start codon within LysIME-EF1 gene, which are demonstrated to be responsible for the translation of the truncated CBD. To elucidate the molecular mechanism for the lytic activity of LysIME-EF1, we combined mutagenesis, lytic activity assays and in vivo animal infection experiments. The results confirmed that the additional LysIME-EF1 CBDs are important for LysIME-EF1 architecture and its lytic activity. To our knowledge, this is the first determined structure of multimeric endolysin encoded by a single gene in E. faecalis phages. As such, it may provide valuable insights into designing potent endolysins against the opportunistic pathogen E. faecalis.

Efficiency of RAPD, ISSR, iPBS, SCoT and phytochemical markers in the genetic relationship study of five native and economical important bamboos of North-East India

10 primers each of random amplified polymorphic DNA (RAPD), inter-simple sequence repeats (ISSR), inter primer binding site (iPBS) and start codon targeted (SCoT) were used to analyze genetic polymorphism and relationship between 50 genotypes of 5 economical important native bamboos (Bambusa cacharensis, B. mizorameana, Dendrocalamus manipureanus, D. hamiltonii and D. sikkimensis) of North-East India. The 40 different primers generated 111, 115, 116 and 138 polymorphic bands for RAPD, ISSR, iPBS and SCoT markers respectively. The comparative analysis of 4 marker systems based on polymorphic information content (PIC), effective multiplex ratio (EMR) and marker index (MI) values showed SCoT to be more informative with higher discriminating power than the other three markers. The correlation value (r) as determined by the Mantel test ranged from 0.60 (SCoT and RAPD) to 0.83 (iPBS and ISSR) indicating a high positive correlation between the markers. The close correspondence between the genetic matrices of RAPD, ISSR, iPBS and SCoT markers revealed the effectiveness of each marker system in determining the genetic relationship between bamboos. UPGMA (Unweighted Pair Group Arithmetic Mean Method) dendrograms generated from DNA marker analysis demonstrated species-specific clustering of different bamboo genotypes. Except for RAPD, the dendrograms of ISSR, iPBS and SCoT markers also showed a close association of bamboo genotypes based on geographical origin. Principal coordinate analysis (PCoA) revealed the distribution of different bamboo genotypes in accordance with the cluster analysis. The cluster grouping based on phytochemical study not only discriminated the different bamboo species but also illustrated a location-specific grouping of the genotypes. The bamboo clustering pattern derived from phytochemical analysis matched closely with the dendrograms generated by the DNA markers. The present investigation established the possibility of using a combined molecular and phytochemical marker approach to determine the genetic relationship between 5 native bamboos of North-East India with high precision.

Single-nucleotide polymorphism characterization of gametocyte development 1 gene in Plasmodium falciparum isolates from Baringo, Uasin Gishu, and Nandi Counties, Kenya

IntroductionPlasmodium falciparum relies on gametocytogenesis to transmit from humans to mosquitoes. Gametocyte development 1 (Pfgdv1) is an upstream activator and epigenetic controller of gametocytogenesis. The emergence of drug resistance is a major public health concern and this requires the development of new strategies that target the transmission of malaria. As a putative drug target, Pfgdv1 has not been characterized to identify its polymorphisms and alleles under selection and how such polymorphisms influence protein structure. MethodsThis study characterized single-nucleotide polymorphisms (SNPs) in primary sequences (n = 30) of Pfgdv1 gene generated from thirty blood samples collected from patients infected with P. falciparum and secondary sequences (n = 216) retrieved from PlasmoDB. ChromasPro, MUSCLE, Tajima's D statistic, SLAC, and STRUM were used in editing raw sequences, performing multiple sequence alignment (MSA), identifying signatures of selection, detecting codon sites under selection pressure, and determining the effect of SNPs, respectively. ResultsMSA of primary and secondary sequences established the existence of five SNPs, consisting of four non-synonymous substitutions (nsSNPs) (p.P217H, p.R398Q, p.H417N, and p.D497E), and a synonymous substitution (p.S514S). The analysis of amino acid changes reveals that p.P217H, p.R398Q, and p.H417N comprise non-conservative changes. Tajima's D statistic showed that these SNPs were under balancing selection, while SLAC analysis identified p.P217H to be under the strongest positive selection. . Further analysis based on thermodynamics indicated that p.P217H has a destabilizing effect, while p.R398Q and p.D497E have stabilizing effects on the protein structure. ConclusionsThe existence of four nsSNPs implies that Pfgdv1 has a minimal diversity in the encoded protein. Selection analysis demonstrates that these nsSNPs are under balancing selection in both local and global populations. However, p.P217H exhibits positive directional selection consistent with previous reports where it showed differentiatial selection of P. falciparum in low and high transmission regions. Therefore, in-silico prediction and experimental determination of protein structure are necessary to evaluate Pfgdv1 as a target candidate for drug design and development.

PCADD: SNV prioritisation in Sus scrofa

BackgroundIn animal breeding, identification of causative genetic variants is of major importance and high economical value. Usually, the number of candidate variants exceeds the number of variants that can be validated. One way of prioritizing probable candidates is by evaluating their potential to have a deleterious effect, e.g. by predicting their consequence. Due to experimental difficulties to evaluate variants that do not cause an amino-acid substitution, other prioritization methods are needed. For human genomes, the prediction of deleterious genomic variants has taken a step forward with the introduction of the combined annotation dependent depletion (CADD) method. In theory, this approach can be applied to any species. Here, we present pCADD (p for pig), a model to score single nucleotide variants (SNVs) in pig genomes.ResultsTo evaluate whether pCADD captures sites with biological meaning, we used transcripts from miRNAs and introns, sequences from genes that are specific for a particular tissue, and the different sites of codons, to test how well pCADD scores differentiate between functional and non-functional elements. Furthermore, we conducted an assessment of examples of non-coding and coding SNVs, which are causal for changes in phenotypes. Our results show that pCADD scores discriminate between functional and non-functional sequences and prioritize functional SNVs, and that pCADD is able to score the different positions in a codon relative to their redundancy. Taken together, these results indicate that based on pCADD scores, regions with biological relevance can be identified and distinguished according to their rate of adaptation.ConclusionsWe present the ability of pCADD to prioritize SNVs in the pig genome with respect to their putative deleteriousness, in accordance to the biological significance of the region in which they are located. We created scores for all possible SNVs, coding and non-coding, for all autosomes and the X chromosome of the pig reference sequence Sscrofa11.1, proposing a toolbox to prioritize variants and evaluate sequences to highlight new sites of interest to explain biological functions that are relevant to animal breeding.

Pervasive Strong Selection at the Level of Codon Usage Bias in Drosophila melanogaster.

Codon usage bias (CUB), where certain codons are used more frequently than expected by chance, is a ubiquitous phenomenon and occurs across the tree of life. The dominant paradigm is that the proportion of preferred codons is set by weak selection. While experimental changes in codon usage have at times shown large phenotypic effects in contrast to this paradigm, genome-wide population genetic estimates have supported the weak selection model. Here we use deep genomic population sequencing of two Drosophila melanogaster populations to measure selection on synonymous sites in a way that allowed us to estimate the prevalence of both weak and strong purifying selection. We find that selection in favor of preferred codons ranges from weak (|Nes| ∼ 1) to strong (|Nes| > 10), with strong selection acting on 10–20% of synonymous sites in preferred codons. While previous studies indicated that selection at synonymous sites could be strong, this is the first study to detect and quantify strong selection specifically at the level of CUB. Further, we find that CUB-associated polymorphism accounts for the majority of strong selection on synonymous sites, with secondary contributions of splicing (selection on alternatively spliced genes, splice junctions, and spliceosome-bound sites) and transcription factor binding. Our findings support a new model of CUB and indicate that the functional importance of CUB, as well as synonymous sites in general, have been underestimated.

The functional role of the C-terminal tail of the human ribosomal protein uS19

The eukaryotic ribosomal protein uS19 has a C-terminal tail that is absent in its bacterial homologue. This tail has been shown to be involved in the formation of the decoding site of human ribosomes. We studied here the previously unexplored functional significance of the 15 C-terminal amino acid residues of human uS19 for the assembly of ribosomes and translation using HEK293-based cell cultures capable of producing FLAG-labeled uS19 (uS19FLAG) or its mutant form deprived of the mentioned amino acid ones. The examination of polysome profiles of cytoplasmic extracts from the respective cells revealed that the deletion of the above uS19 amino acid residues barely affected the assembly and maturation of 40S subunits and the initiation of translation, but completely prevented the formation of polysomes. This implied the crucial importance of the uS19 tail in the elongation process. Analysis of tRNAs associated with 40S subunits and 80S ribosomes containing wild type uS19FLAG or its truncated form showed that the deletion of the C-terminal pentadecapeptide fragment of uS19 did not interfere with the binding of aminoacyl-tRNA (aa-tRNA) at the ribosomal A site. The results led to the conclusion that the transpeptidation, which occurs on the large ribosomal subunit after decoding the A site codon by the incoming aa-tRNA, is the most likely elongation stage, where this uS19 fragment can play a critical role. Our findings suggest that the uS19 tail is a keystone player in the accommodation of aa-tRNA at the A site, which is a pre-requisite for the peptide transfer.

In silicoEvaluation of Substrate Binding Site and Rare Codons in the Structure of CYP152A1

Background:The Cytochromes P450 (CYPs) have an essential role in the oxidation of endogenous and exogenous molecules. The CYPs are identified in all domains of life, but the CYP152A1 from Bacillus subtilis is specially considered for clinical and industrial applications. The molecular cloning of a new type of CYP from Bacillus subtilis was reported, previously. Here, we describe the hidden layer of biological information of the CYP152A1 enzyme, which can help researchers for better understanding of enzyme application. In this study, four rare codons of enzyme, including Arg63, Arg187, Arg276, and Arg338 were identified and evaluated using the bioinformatics web servers.Methods:Through in silico modeling of CYP152A1 via the I-TASSER server, the above-mentioned rare codons were studied in the structure of enzyme that may have an important role in the proper folding of CYP152A1. In the following, the substrate binding site of CYP152A1 was studied by AutoDock Vina, and the heme and palmitic acid were considered as the substrates.Results:The results of docking study elucidated the Arg242 in the active site is closely related to the substrate binding site of CYP152A1, which help us to further clarify the mechanism of the enzyme reaction.Conclusion:Studies of these hidden information’s can enhance our understanding of CYP152A1 folding and protein expression challenges. Moreover, identification of rare codons can help in the rational design of new and effective drugs.

Major Revisions in Arthropod Phylogeny Through Improved Supermatrix, With Support for Two Possible Waves of Land Invasion by Chelicerates.

Deep phylogeny involving arthropod lineages is difficult to recover because the erosion of phylogenetic signals over time leads to unreliable multiple sequence alignment (MSA) and subsequent phylogenetic reconstruction. One way to alleviate the problem is to assemble a large number of gene sequences to compensate for the weakness in each individual gene. Such an approach has led to many robustly supported but contradictory phylogenies. A close examination shows that the supermatrix approach often suffers from two shortcomings. The first is that MSA is rarely checked for reliability and, as will be illustrated, can be poor. The second is that, to alleviate the problem of homoplasy at the third codon position of protein-coding genes due to convergent evolution of nucleotide frequencies, phylogeneticists may remove or degenerate the third codon position but may do it improperly and introduce new biases. We performed extensive reanalysis of one of such “big data” sets to highlight these two problems, and demonstrated the power and benefits of correcting or alleviating these problems. Our results support a new group with Xiphosura and Arachnopulmonata (Tetrapulmonata + Scorpiones) as sister taxa. This favors a new hypothesis in which the ancestor of Xiphosura and the extinct Eurypterida (sea scorpions, of which many later forms lived in brackish or freshwater) returned to the sea after the initial chelicerate invasion of land. Our phylogeny is supported even with the original data but processed with a new “principled” codon degeneration. We also show that removing the 1673 codon sites with both AGN and UCN codons (encoding serine) in our alignment can partially reconcile discrepancies between nucleotide-based and AA-based tree, partly because two sequences, one with AGN and the other with UCN, would be identical at the amino acid level but quite different at the nucleotide level.

Bat IFITM3 restriction depends on S-palmitoylation and a polymorphic site within the CD225 domain.

Host interferon-induced transmembrane proteins (IFITMs) are broad-spectrum antiviral restriction factors. Of these, IFITM3 potently inhibits viruses that enter cells through acidic endosomes, many of which are zoonotic and emerging viruses with bats (order Chiroptera) as their natural hosts. We previously demonstrated that microbat IFITM3 is antiviral. Here, we show that bat IFITMs are characterized by strong adaptive evolution and identify a highly variable and functionally important site-codon 70-within the conserved CD225 domain of IFITMs. Mutation of this residue in microbat IFITM3 impairs restriction of representatives of four different virus families that enter cells via endosomes. This mutant shows altered subcellular localization and reduced S-palmitoylation, a phenotype copied by mutation of conserved cysteine residues in microbat IFITM3. Furthermore, we show that microbat IFITM3 is S-palmitoylated on cysteine residues C71, C72, and C105, mutation of each cysteine individually impairs virus restriction, and a triple C71A-C72A-C105A mutant loses all restriction activity, concomitant with subcellular re-localization of microbat IFITM3 to Golgi-associated sites. Thus, we propose that S-palmitoylation is critical for Chiropteran IFITM3 function and identify a key molecular determinant of IFITM3 S-palmitoylation.

MRNA regions where 80S ribosomes pause during translation elongation in vivo interact with protein uS19, a component of the decoding site.

In eukaryotic ribosomes, the conserved protein uS19, formerly known as S15, extends with its C-terminal tail to the decoding site. The cross-linking of uS19 to the A site codon has been detected using synthetic mRNAs bearing 4-thiouridine (s4U) residues. Here, we showed that the A-site tRNA prevents this cross-linking and that the P site codon does not contact uS19. Next, we focused on determining uS19-mRNA interactions in vivo by applying the photoactivatable-ribonucleoside enhancing cross-linking and immunoprecipitation method to a stable HEK293 cell line producing FLAG-tagged uS19 and grown in a medium containing s4U. We found that when translation was stopped by cycloheximide, uS19 was efficiently cross-linked to mRNA regions with a high frequency of Glu, Lys and, more rarely, Arg codons. The results indicate that the complexes, in which the A site codon is not involved in the formation of the mRNA-tRNA duplex, are present among the cycloheximide-arrested 80S complexes, which implies pausing of elongating ribosomes at the above mRNA regions. Thus, our findings demonstrate that the human ribosomal protein uS19 interacts with mRNAs during translation elongation and highlight the regions of mRNAs where ribosome pausing occurs, bringing new structural and functional insights into eukaryotic translation in vivo.

Episodic positive diversifying selection on key immune system genes in major avian lineages.

The major histocompatibility complex (MHC) of the adaptive immune system and the toll-like receptor (TLR) family of the innate immune system are involved in the detection of foreign invaders, and thus are subject to parasite-driven molecular evolution. Herein, we tested for macroevolutionary signatures of selection in these gene families within and among all three major clades of birds (Paleognathae, Galloanserae, and Neoaves). We characterized evolutionary relationships of representative immune genes (Mhc1 and Tlr2b) and a control gene (ubiquitin, Ubb), using a relatively large and phylogenetically diverse set of species with complete coding sequences (34 orthologous loci for Mhc1, 29 for Tlr2b, and 37 for Ubb). Episodic positive diversifying selection was found in the gene-wide phylogenies of the two immune genes, as well as at specific sites within each gene (8.5% of codon sites in Mhc1 and 2.7% in Tlr2b), but not in the control gene (Ubb). We found 20% of lineages under episodic diversifying selection in Mhc1 versus 9.1% in Tlr2b. For Mhc1, selection was relaxed in the Galloanserae and intensified in the Neoaves relative to the other clades, but no differences were detected among clades in the Tlr2b gene. In summary, we provide evidence of episodic positive diversifying selection in key immune genes and demonstrate differential strengths of selection within Class Aves, with the adaptive gene showing an increased divergence and evolutionary rate over the innate gene, contributing to the growing understanding of vertebrate immune gene evolution.

Human cells adapt to translational errors by modulating protein synthesis rate and protein turnover

ABSTRACTDeregulation of tRNAs, aminoacyl-tRNA synthetases (aaRS) or tRNA modifying enzymes, increase the level of protein synthesis errors (PSE) and are associated with several diseases, but the cause-effect mechanisms of these pathologies remain elusive. To clarify the role of PSE in human biology, we have engineered a HEK293 cell line to overexpress a wild type (Wt) tRNASer and two tRNASer mutants that misincorporate serine at non-cognate codon sites. Then, we followed long-term adaptation to PSE of such recombinant cells by analysing cell viability, protein synthesis rate and activation of protein quality control mechanisms (PQC). Engineered cells showed higher level of misfolded and aggregated proteins; activated the ubiquitin-proteasome system (UPS) and the unfolded protein response (UPR), indicative of proteotoxic stress. Adaptation to PSE involved increased protein turnover, UPR up-regulation and altered protein synthesis rate. Gene expression analysis showed that engineered cells presented recurrent alterations in the endoplasmic reticulum, cell adhesion and calcium homeostasis. Herein, we unveil new phenotypic consequences of protein synthesis errors in human cells and identify the protein quality control processes that are necessary for long-term adaptation to PSE and proteotoxic stress. Our data provide important insight on how chronic proteotoxic stress may cause disease and highlight potential biological pathways that support the association of PSE with disease.