Translational Selection Is Ubiquitous in Prokaryotes

Codon usage bias in prokaryotic genomes is largely a consequence of background substitution patterns in DNA, but highly expressed genes may show a preference towards codons that enable more efficient and/or accurate translation. We introduce a novel approach based on supervised machine learning that detects effects of translational selection on genes, while controlling for local variation in nucleotide substitution patterns represented as sequence composition of intergenic DNA. A cornerstone of our method is a Random Forest classifier that outperformed previous distance measure-based approaches, such as the codon adaptation index, in the task of discerning the (highly expressed) ribosomal protein genes by their codon frequencies. Unlike previous reports, we show evidence that translational selection in prokaryotes is practically universal: in 460 of 461 examined microbial genomes, we find that a subset of genes shows a higher codon usage similarity to the ribosomal proteins than would be expected from the local sequence composition. These genes constitute a substantial part of the genome—between 5% and 33%, depending on genome size—while also exhibiting higher experimentally measured mRNA abundances and tending toward codons that match tRNA anticodons by canonical base pairing. Certain gene functional categories are generally enriched with, or depleted of codon-optimized genes, the trends of enrichment/depletion being conserved between Archaea and Bacteria. Prominent exceptions from these trends might indicate genes with alternative physiological roles; we speculate on specific examples related to detoxication of oxygen radicals and ammonia and to possible misannotations of asparaginyl–tRNA synthetases. Since the presence of codon optimizations on genes is a valid proxy for expression levels in fully sequenced genomes, we provide an example of an “adaptome” by highlighting gene functions with expression levels elevated specifically in thermophilic Bacteria and Archaea

Evolutionary Rate Covariation Identifies New Members of a Protein Network Required for Drosophila melanogaster Female Post-Mating Responses

Seminal fluid proteins transferred from males to females during copulation are required for full fertility and can exert dramatic effects on female physiology and behavior. In Drosophila melanogaster, the seminal protein sex peptide (SP) affects mated females by increasing egg production and decreasing receptivity to courtship. These behavioral changes persist for several days because SP binds to sperm that are stored in the female. SP is then gradually released, allowing it to interact with its female-expressed receptor. The binding of SP to sperm requires five additional seminal proteins, which act together in a network. Hundreds of uncharacterized male and female proteins have been identified in this species, but individually screening each protein for network function would present a logistical challenge. To prioritize the screening of these proteins for involvement in the SP network, we used a comparative genomic method to identify candidate proteins whose evolutionary rates across the Drosophila phylogeny co-vary with those of the SP network proteins. Subsequent functional testing of 18 co-varying candidates by RNA interference identified three male seminal proteins and three female reproductive tract proteins that are each required for the long-term persistence of SP responses in females. Molecular genetic analysis showed the three new male proteins are required for the transfer of other network proteins to females and for SP to become bound to sperm that are stored in mated females. The three female proteins, in contrast, act downstream of SP binding and sperm storage. These findings expand the number of seminal proteins required for SP's actions in the female and show that multiple female proteins are necessary for the SP response. Furthermore, our functional analyses demonstrate that evolutionary rate covariation is a valuable predictive tool for identifying candidate members of interacting protein networks. © 2014 Findlay et al

Multiple novel prostate cancer susceptibility signals identified by fine-mapping of known risk loci among Europeans

Genome-wide association studies (GWAS) have identified numerous common prostate cancer (PrCa) susceptibility loci. We have fine-mapped 64 GWAS regions known at the conclusion of the iCOGS study using large-scale genotyping and imputation in 25 723 PrCa cases and 26 274 controls of European ancestry. We detected evidence for multiple independent signals at 16 regions, 12 of which contained additional newly identified significant associations. A single signal comprising a spectrum of correlated variation was observed at 39 regions; 35 of which are now described by a novel more significantly associated lead SNP, while the originally reported variant remained as the lead SNP only in 4 regions. We also confirmed two association signals in Europeans that had been previously reported only in East-Asian GWAS. Based on statistical evidence and linkage disequilibrium (LD) structure, we have curated and narrowed down the list of the most likely candidate causal variants for each region. Functional annotation using data from ENCODE filtered for PrCa cell lines and eQTL analysis demonstrated significant enrichment for overlap with bio-features within this set. By incorporating the novel risk variants identified here alongside the refined data for existing association signals, we estimate that these loci now explain ∼38.9% of the familial relative risk of PrCa, an 8.9% improvement over the previously reported GWAS tag SNPs. This suggests that a significant fraction of the heritability of PrCa may have been hidden during the discovery phase of GWAS, in particular due to the presence of multiple independent signals within the same regio

Functional interplay of DnaE polymerase, DnaG primase and DnaC helicase within a ternary complex, and primase to polymerase hand-off during lagging strand DNA replication in Bacillus subtilis

Bacillus subtilis has two replicative DNA polymerases. PolC is a processive high-fidelity replicative polymerase, while the error-prone DnaEBs extends RNA primers before hand-off to PolC at the lagging strand. We show that DnaEBs interacts with the replicative helicase DnaC and primase DnaG in a ternary complex. We characterize their activities and analyse the functional significance of their interactions using primase, helicase and primer extension assays, and a ‘stripped down’ reconstituted coupled assay to investigate the coordinated displacement of the parental duplex DNA at a replication fork, synthesis of RNA primers along the lagging strand and hand-off to DnaEBs. The DnaG– DnaEBs hand-off takes place after de novo polymerization of only two ribonucleotides by DnaG, and does not require other replication proteins. Furthermore, the fidelity of DnaEBs is improved by DnaC and DnaG, likely via allosteric effects induced by direct protein–protein interactions that lower the efficiency of nucleotide mis-incorporations and/or the efficiency of extension of mis-aligned primers in the catalytic site of DnaEBs. We conclude that de novo RNA primer synthesis by DnaG and initial primer extension by DnaEBs are carried out by a lagging strand–specific subcomplex comprising DnaG, DnaEBs and DnaC, which stimulates chromosomal replication with enhanced fidelity

A KCNK16 mutation causing TALK-1 gain-of-function is associated with maturity-onset diabetes of the young

Maturity-onset diabetes of the young (MODY) is a heterogeneous group of monogenic disorders of impaired pancreatic β cell function. The mechanisms underlying MODY include β cell KATP channel dysfunction (e.g., KCNJ11 [MODY13] or ABCC8 [MODY12] mutations); however, no other β cell channelopathies have been associated with MODY to date. Here, we have identified a nonsynonymous coding variant in KCNK16 (NM_001135105: c.341T>C, p.Leu114Pro) segregating with MODY. KCNK16 is the most abundant and β cell–restricted K(+) channel transcript, encoding the two-pore-domain K(+) channel TALK-1. Whole-cell K(+) currents demonstrated a large gain of function with TALK-1 Leu114Pro compared with TALK-1 WT, due to greater single-channel activity. Glucose-stimulated membrane potential depolarization and Ca(2+) influx were inhibited in mouse islets expressing TALK-1 Leu114Pro with less endoplasmic reticulum Ca(2+) storage. TALK-1 Leu114Pro significantly blunted glucose-stimulated insulin secretion compared with TALK-1 WT in mouse and human islets. These data suggest that KCNK16 is a previously unreported gene for MODY