Retapamulin-Assisted Ribosome Profiling Reveals the Alternative Bacterial Proteome.

The use of alternative translation initiation sites enables production of more than one protein from a single gene, thereby expanding the cellular proteome. Although several such examples have been serendipitously found in bacteria, genome-wide mapping of alternative translation start sites has been unattainable. We found that the antibiotic retapamulin specifically arrests initiating ribosomes at start codons of the genes. Retapamulin-enhanced Ribo-seq analysis (Ribo-RET) not only allowed mapping of conventional initiation sites at the beginning of the genes, but strikingly, it also revealed putative internal start sites in a number of Escherichia coli genes. Experiments demonstrated that the internal start codons can be recognized by the ribosomes and direct translation initiation in vitro and in vivo. Proteins, whose synthesis is initiated at internal in-frame and out-of-frame start sites, can be functionally important and contribute to the "alternative" bacterial proteome. The internal start sites may also play regulatory roles in gene expression.

Programmed Ribosomal Frameshifting Generates a Copper Transporter and a Copper Chaperone from the Same Gene.

Metal efflux pumps maintain ion homeostasis in the cell. The functions of the transporters are often supported by chaperone proteins, which scavenge the metal ions from the cytoplasm. Although the copper ion transporter CopA has been known in Escherichia coli, no gene for its chaperone had been identified. We show that the CopA chaperone is expressed in E. coli from the same gene that encodes the transporter. Some ribosomes translating copA undergo programmed frameshifting, terminate translation in the -1 frame, and generate the 70 aa-long polypeptide CopA(Z), which helps cells survive toxic copper concentrations. The high efficiency of frameshifting is achieved by the combined stimulatory action of a "slippery" sequence, an mRNA pseudoknot, and the CopA nascent chain. Similar mRNA elements are not only found in the copA genes of other bacteria but are also present in ATP7B, the human homolog of copA, and direct ribosomal frameshifting in vivo.

A role for the Saccharomyces cerevisiae ABCF protein New1 in translation termination/recycling.

Translation is controlled by numerous accessory proteins and translation factors. In the yeast Saccharomyces cerevisiae, translation elongation requires an essential elongation factor, the ABCF ATPase eEF3. A closely related protein, New1, is encoded by a non-essential gene with cold sensitivity and ribosome assembly defect knock-out phenotypes. Since the exact molecular function of New1 is unknown, it is unclear if the ribosome assembly defect is direct, i.e. New1 is a bona fide assembly factor, or indirect, for instance due to a defect in protein synthesis. To investigate this, we employed yeast genetics, cryo-electron microscopy (cryo-EM) and ribosome profiling (Ribo-Seq) to interrogate the molecular function of New1. Overexpression of New1 rescues the inviability of a yeast strain lacking the otherwise strictly essential translation factor eEF3. The structure of the ATPase-deficient (EQ2) New1 mutant locked on the 80S ribosome reveals that New1 binds analogously to the ribosome as eEF3. Finally, Ribo-Seq analysis revealed that loss of New1 leads to ribosome queuing upstream of 3'-terminal lysine and arginine codons, including those genes encoding proteins of the cytoplasmic translational machinery. Our results suggest that New1 is a translation factor that fine-tunes the efficiency of translation termination or ribosome recycling.

Prevalence of drug resistance mutations in HAART patients infected with HIV-1 CRF06_cpx in Estonia.

HIV-1 drug resistance mutations (DRMs) and substitutions were assessed after the failure of the first line non-nucleoside reverse transcriptase inhibitors (NNRTIs) + 2 nucleoside/nucleotide reverse transcriptase inhibitors (NRTIs) treatment regimens (efavirenz [EFV] + lamivudine[3TC] + zidovudine [ZDV] vs. EFV + 3TC + ddI) among the HIV-1 CRF06_cpx infected subjects in Estonia. HIV-1 genomic RNA was sequenced; DRMs and amino acid substitutions were compared in 44 treatment naïve and 45 first-line NNRTI + 2 NRTI treatment failed patients consisting of EFV + 3TC + ZDV (n = 17) and EFV + 3TC + didanosine[ddI] (n = 21) therapy failed sub-populations. At least one DRM was found in 78% of treatment experienced patients. The most common NRTI mutations were M184V (80%), L74V (31%), L74I (17%), K219E (9%), and M184I (9%), NNRTI mutations were K103N (83%), P225H (14%), L100I (11%), and Y188L (11%), reflecting generally the similar pattern of DRMs to that seen in treatment failed subtype B viruses. Sub-population analysis revealed that EFV + 3TC + ddI failed patients had more DRMs compared to EFV + 3TC + ZDV failed patients, especially the ddI DRM L74IV and several additional NNRTI DRMs. Additionally, CRF06_cpx specific mutation E179V and substitutions R32K, K122E, and V200AE were also detected in treatment experienced population. After the failure of the first-line EFV + 3TC + ddI therapy HIV-1 CRF06_cpx viruses develop additional NRTI and NNRTI mutations compared to EFV + 3TC + ZDV regimen. Therefore the usage of EFV + 3TC + ddI in this subtype decreases the options for next regimens containing abacavir, and NNRTI class agents.

The ABCF ATPase New1 resolves translation termination defects associated with specific tRNAArg and tRNALys isoacceptors in the P site.

The efficiency of translation termination is determined by the nature of the stop codon as well as its context. In eukaryotes, recognition of the A-site stop codon and release of the polypeptide are mediated by release factors eRF1 and eRF3, respectively. Translation termination is modulated by other factors which either directly interact with release factors or bind to the E-site and modulate the activity of the peptidyl transferase center. Previous studies suggested that the Saccharomyces cerevisiae ABCF ATPase New1 is involved in translation termination and/or ribosome recycling, however, the exact function remained unclear. Here, we have applied 5PSeq, single-particle cryo-EM and readthrough reporter assays to provide insight into the biological function of New1. We show that the lack of New1 results in ribosomal stalling at stop codons preceded by a lysine or arginine codon and that the stalling is not defined by the nature of the C-terminal amino acid but rather by the identity of the tRNA isoacceptor in the P-site. Collectively, our results suggest that translation termination is inefficient when ribosomes have specific tRNA isoacceptors in the P-site and that the recruitment of New1 rescues ribosomes at these problematic termination contexts.

The Escherichia coli mqsR and ygiT genes encode a new toxin-antitoxin pair.

Toxin-antitoxin (TA) systems are plasmid- or chromosome-encoded protein complexes composed of a stable toxin and a short-lived inhibitor of the toxin. In cultures of Escherichia coli, transcription of toxin-antitoxin genes was induced in a nondividing subpopulation of bacteria that was tolerant to bactericidal antibiotics. Along with transcription of known toxin-antitoxin operons, transcription of mqsR and ygiT, two adjacent genes with multiple TA-like features, was induced in this cell population. Here we show that mqsR and ygiT encode a toxin-antitoxin system belonging to a completely new family which is represented in several groups of bacteria. The mqsR gene encodes a toxin, and ectopic expression of this gene inhibits growth and induces rapid shutdown of protein synthesis in vivo. ygiT encodes an antitoxin, which protects cells from the effects of MqsR. These two genes constitute a single operon which is transcriptionally repressed by the product of ygiT. We confirmed that transcription of this operon is induced in the ampicillin-tolerant fraction of a growing population of E. coli and in response to activation of the HipA toxin. Expression of the MqsR toxin does not kill bacteria but causes reversible growth inhibition and elongation of cells.

Ribosome profiling analysis of eEF3-depleted Saccharomyces cerevisiae.

In addition to the standard set of translation factors common in eukaryotic organisms, protein synthesis in the yeast Saccharomyces cerevisiae requires an ABCF ATPase factor eEF3, eukaryotic Elongation Factor 3. eEF3 is an E-site binder that was originally identified as an essential factor involved in the elongation stage of protein synthesis. Recent biochemical experiments suggest an additional function of eEF3 in ribosome recycling. We have characterised the global effects of eEF3 depletion on translation using ribosome profiling. Depletion of eEF3 results in decreased ribosome density at the stop codon, indicating that ribosome recycling does not become rate limiting when eEF3 levels are low. Consistent with a defect in translation elongation, eEF3 depletion causes a moderate redistribution of ribosomes towards the 5' part of the open reading frames. We observed no E-site codon- or amino acid-specific ribosome stalling upon eEF3 depletion, supporting its role as a general elongation factor. Surprisingly, depletion of eEF3 leads to a relative decrease in P-site proline stalling, which we hypothesise is a secondary effect of generally decreased translation and/or decreased competition for the E-site with eIF5A.

Chorionic gonadotropin beta-gene variants are associated with recurrent miscarriage in two European populations.

CONTEXT: The incidence of recurrent miscarriage (RM) (>or=3 consecutive pregnancy losses) is estimated as 1-2% in fertile couples. Familial clustering of RM has suggested the contribution of a genetic component. OBJECTIVE: A low level of human chorionic gonadotropin (HCG) in maternal serum during the first trimester of the pregnancy is a clinically accepted risk factor for miscarriage. We sought to study whether variation in chorionic gonadotropin beta-subunit genes (CGBs) expressed in placenta may contribute to the risk of RM. DESIGN: Resequencing of CGB5 and CGB8, the two most actively transcribed loci of the four HCG beta-duplicate genes, was performed. SETTING: A case-control study involving two sample sets, from Estonia (n = 194) and Finland (n = 185), was performed. PATIENTS: RM patients (n = 184) and fertile controls (n = 195) participated in the study. RESULTS: From 71 identified variants in CGB5 and CGB8, 48 polymorphisms were novel. Significant protective effect was associated with two single nucleotide polymorphisms located at identical positions in intron 2 in both CGB5 [P = 0.007; odds ratio (OR) = 0.53] and CGB8 (P = 0.042; OR = 0.15), and with four CGB5 promoter variants (P < 0.03; OR = 0.54-0.58). The carriers of minor alleles had a reduced risk of RM. The haplotype structure of the CGB8 promoter was consistent with balancing selection; a rare mutation in CGB8 initiator element was detected only among patients (n = 3). In addition, three rare nonsynonymous substitutions were identified among RM cases as possible variants increasing the risk of recurrent pregnancy loss. CONCLUSION: The findings encourage studying the functional effect of the identified variants on CGB expression and HCG hormone activity to elucidate further the role of CGB variation in RM.

Phylogenetic distribution of translational GTPases in bacteria.

BACKGROUND: Translational GTPases are a family of proteins in which GTPase activity is stimulated by the large ribosomal subunit. Conserved sequence features allow members of this family to be identified. RESULTS: To achieve accurate protein identification and grouping we have developed a method combining searches with Hidden Markov Model profiles and tree based grouping. We found all the genes for translational GTPases in 191 fully sequenced bacterial genomes. The protein sequences were grouped into nine subfamilies. Analysis of the results shows that three translational GTPases, the translation factors EF-Tu, EF-G and IF2, are present in all organisms examined. In addition, several copies of the genes encoding EF-Tu and EF-G are present in some genomes. In the case of multiple genes for EF-Tu, the gene copies are nearly identical; in the case of multiple EF-G genes, the gene copies have been considerably diverged. The fourth translational GTPase, LepA, the function of which is currently unknown, is also nearly universally conserved in bacteria, being absent from only one organism out of the 191 analyzed. The translation regulator, TypA, is also present in most of the organisms examined, being absent only from bacteria with small genomes.Surprisingly, some of the well studied translational GTPases are present only in a very small number of bacteria. The translation termination factor RF3 is absent from many groups of bacteria with both small and large genomes. The specialized translation factor for selenocysteine incorporation--SelB--was found in only 39 organisms. Similarly, the tetracycline resistance proteins (Tet) are present only in a small number of species. Proteins of the CysN/NodQ subfamily have acquired functions in sulfur metabolism and production of signaling molecules. The genes coding for CysN/NodQ proteins were found in 74 genomes. This protein subfamily is not confined to Proteobacteria, as suggested previously but present also in many other groups of bacteria. CONCLUSION: Four of the translational GTPase subfamilies (IF2, EF-Tu, EF-G and LepA) are represented by at least one member in each bacterium studied, with one exception in LepA. This defines the set of translational GTPases essential for basic cell functions.

Association between HIV-1 tropism and CCR5 human haplotype E in a Caucasian population.

BACKGROUND: The influence of the diversity of CCR5 on HIV susceptibility and disease progression has been clearly demonstrated but how the variability of this gene influences the HIV tropism is poorly understood. We investigated whether CCR5 haplotypes are associated with HIV tropism in a Caucasian population. METHODS: We evaluated 161 HIV-positive subjects in a cross-sectional study. CCR5 haplotypes were derived after genotyping 9 CCR2-CCR5 polymorphisms. The HIV subtype was determined by phylogenetic analysis using the maximum likelihood method and viral tropism by the genotypic tropism assay (geno2pheno). Associations between CCR5 haplotypes and viral tropism were determined using logistic regression analyses. Samples from 500 blood donors were used to evaluate the representativeness of HIV-positives in terms of CCR5 haplotype distribution. RESULTS: The distribution of CCR5 haplotypes was similar in HIV-positive subjects and blood donors. The majority of viruses (93.8%) belonged to HIV-1 CRF06_cpx; 7.5% were X4, and the remaining were R5 tropic. X4 tropic viruses were over represented among people with CCR5 human haplotype E (HHE) compared with those without this haplotype (13.0% vs 1.4%; P = 0.006). People possessing CCR5 HHE had 11 times increased odds (odds ratio = 11.00; 95% confidence interval: 1.38 to 87.38) of having X4 tropic viruses than those with non-HHE. After adjusting for antiretroviral (ARV) therapy, neither the presence of HHE nor the use of ARV was associated with X4 tropic viruses. CONCLUSIONS: Our results suggest that CCR5 HHE and ARV treatment might be associated with the presence of HIV-1 X4 tropic viruses.

A modest but significant effect of CGB5 gene promoter polymorphisms in modulating the risk of recurrent miscarriage.

OBJECTIVE: To confirm the effect of single nucleotide polymorphisms (SNPs) in chorionic gonadotropin beta (CGB) genes in modulating the susceptibility to recurrent miscarriage (RM) in Danes and in a meta-analysis across Danes and the discovery samples from Estonia and Finland. DESIGN: Case-control association study, restriction fragment length polymorphism genotyping, resequencing. SETTING: Fertility clinics at the Rigshospitalet, Copenhagen, and Aalborg Hospital, Aalborg, Denmark. PATIENT(S): Four hundred fifty Danish women and men from couples with RM and 119 women with children and no miscarriages in new study. A total of 634 women and men from RM couples and 314 female controls in a combined study of Estonians, Finns, and Danes. INTERVENTION(S): None. MAIN OUTCOME MEASURE(S): Distribution of CGB5 and CGB8 allele and haplotype frequencies in patients and controls. RESULT(S): For the majority of studied SNPs, the allelic and haplotypic distribution differed statistically between the Danish and the previous Estonian-Finnish sample. In Danes, two CGB5 promoter SNPs (c5-155; c5-142) exhibited a nonsignificant trend for higher allele frequency in fertile women compared with RM patients. The meta-analysis of results from three populations confirmed a modest but significant effect on carriage of c5-155C (odds ratio = 0.64; 95% confidence interval [CI] 0.44-0.94) and c5-142A (odds ratio = 0.66; 95% CI, 0.45-0.94) variants in reducing the risk of RM. None of the investigated genetic variants in the CGB8 gene was associated with RM. CONCLUSION(S): Carriage of particular variants in the promoter of the CGB5 gene seems to protect against RM. No common genetic variants in CGB5 and CGB8 were associated with increased RM susceptibility in the studied North European populations.

A computational study of elongation factor G (EFG) duplicated genes: diverged nature underlying the innovation on the same structural template.

BACKGROUND: Elongation factor G (EFG) is a core translational protein that catalyzes the elongation and recycling phases of translation. A more complex picture of EFG's evolution and function than previously accepted is emerging from analyzes of heterogeneous EFG family members. Whereas the gene duplication is postulated to be a prominent factor creating functional novelty, the striking divergence between EFG paralogs can be interpreted in terms of innovation in gene function. METHODOLOGY/PRINCIPAL FINDINGS: We present a computational study of the EFG protein family to cover the role of gene duplication in the evolution of protein function. Using phylogenetic methods, genome context conservation and insertion/deletion (indel) analysis we demonstrate that the EFG gene copies form four subfamilies: EFG I, spdEFG1, spdEFG2, and EFG II. These ancient gene families differ by their indispensability, degree of divergence and number of indels. We show the distribution of EFG subfamilies and describe evidences for lateral gene transfer and recent duplications. Extended studies of the EFG II subfamily concern its diverged nature. Remarkably, EFG II appears to be a widely distributed and a much-diversified subfamily whose subdivisions correlate with phylum or class borders. The EFG II subfamily specific characteristics are low conservation of the GTPase domain, domains II and III; absence of the trGTPase specific G2 consensus motif "RGITI"; and twelve conserved positions common to the whole subfamily. The EFG II specific functional changes could be related to changes in the properties of nucleotide binding and hydrolysis and strengthened ionic interactions between EFG II and the ribosome, particularly between parts of the decoding site and loop I of domain IV. CONCLUSIONS/SIGNIFICANCE: Our work, for the first time, comprehensively identifies and describes EFG subfamilies and improves our understanding of the function and evolution of EFG duplicated genes.

Explaining the imperfection of the molecular clock of hominid mitochondria.

The molecular clock of mitochondrial DNA has been extensively used to date various genetic events. However, its substitution rate among humans appears to be higher than rates inferred from human-chimpanzee comparisons, limiting the potential of interspecies clock calibrations for intraspecific dating. It is not well understood how and why the substitution rate accelerates. We have analyzed a phylogenetic tree of 3057 publicly available human mitochondrial DNA coding region sequences for changes in the ratios of mutations belonging to different functional classes. The proportion of non-synonymous and RNA genes substitutions has reduced over hundreds of thousands of years. The highest mutation ratios corresponding to fast acceleration in the apparent substitution rate of the coding sequence have occurred after the end of the Last Ice Age. We recalibrate the molecular clock of human mtDNA as 7990 years per synonymous mutation over the mitochondrial genome. However, the distribution of substitutions at synonymous sites in human data significantly departs from a model assuming a single rate parameter and implies at least 3 different subclasses of sites. Neutral model with 3 synonymous substitution rates can explain most, if not all, of the apparent molecular clock difference between the intra- and interspecies levels. Our findings imply the sluggishness of purifying selection in removing the slightly deleterious mutations from the human as well as the Neandertal and chimpanzee populations. However, for humans, the weakness of purifying selection has been further exacerbated by the population expansions associated with the out-of Africa migration and the end of the Last Ice Age.

Segmental duplications and gene conversion: Human luteinizing hormone/chorionic gonadotropin beta gene cluster.

Segmental duplicons (>1 kb) of high sequence similarity (>90%) covering >5% of the human genome are characterized by complex sequence variation. Apart from a few well-characterized regions (MHC, beta-globin), the diversity and linkage disequilibrium (LD) patterns of duplicons and the role of gene conversion in shaping them have been poorly studied. To shed light on these issues, we have re-sequenced the human Luteinizing Hormone/Chorionic Gonadotropin beta (LHB/CGB) cluster (19q13.32) of three population samples (Estonians, Mandenka, and Han). The LHB/CGB cluster consists of seven duplicated genes critical in human reproduction. In the LHB/CGB region, high sequence diversity, concentration of gene-conversion acceptor sites, and strong LD colocalize with peripheral genes, whereas central loci are characterized by lower variation, gene-conversion donor activity, and breakdown of LD between close markers. The data highlight an important role of gene conversion in spreading polymorphisms among duplicon copies and generating LD around them. The directionality of gene-conversion events seems to be determined by the localization of a predicted recombination "hotspot" and "warm spot" in the vicinity of the most active acceptor genes at the periphery of the cluster. The data suggest that enriched crossover activity in direct and inverted segmental repeats is in accordance with the formation of palindromic secondary structures promoting double-strand breaks rather than fixed DNA sequence motifs. Also, this first detailed coverage of sequence diversity and structure of the LHB/CGB gene cluster will pave the way for studying the identified polymorphisms as well as potential genomic rearrangements in association with an individual's reproductive success.