SNVstory: inferring genetic ancestry from genome sequencing data.

BACKGROUND: Genetic ancestry, inferred from genomic data, is a quantifiable biological parameter. While much of the human genome is identical across populations, it is estimated that as much as 0.4% of the genome can differ due to ancestry. This variation is primarily characterized by single nucleotide variants (SNVs), which are often unique to specific genetic populations. Knowledge of a patient's genetic ancestry can inform clinical decisions, from genetic testing and health screenings to medication dosages, based on ancestral disease predispositions. Nevertheless, the current reliance on self-reported ancestry can introduce subjectivity and exacerbate health disparities. While genomic sequencing data enables objective determination of a patient's genetic ancestry, existing approaches are limited to ancestry inference at the continental level. RESULTS: To address this challenge, and create an objective, measurable metric of genetic ancestry we present SNVstory, a method built upon three independent machine learning models for accurately inferring the sub-continental ancestry of individuals. We also introduce a novel method for simulating individual samples from aggregate allele frequencies from known populations. SNVstory includes a feature-importance scheme, unique among open-source ancestral tools, which allows the user to track the ancestral signal broadcast by a given gene or locus. We successfully evaluated SNVstory using a clinical exome sequencing dataset, comparing self-reported ethnicity and race to our inferred genetic ancestry, and demonstrate the capability of the algorithm to estimate ancestry from 36 different populations with high accuracy. CONCLUSIONS: SNVstory represents a significant advance in methods to assign genetic ancestry, opening the door to ancestry-informed care. SNVstory, an open-source model, is packaged as a Docker container for enhanced reliability and interoperability. It can be accessed from https://github.com/nch-igm/snvstory .

Amino acid substitutions in human growth hormone affect secondary structure and receptor binding.

The interaction between human Growth Hormone (hGH) and hGH Receptor (hGHR) has basic relevance to cancer and growth disorders, and hGH is the scaffold for Pegvisomant, an anti-acromegaly therapeutic. For the latter reason, hGH has been extensively engineered by early workers to improve binding and other properties. We are particularly interested in E174 which belongs to the hGH zinc-binding triad; the substitution E174A is known to significantly increase binding, but to now no explanation has been offered. We generated this and several computationally-selected single-residue substitutions at the hGHR-binding site of hGH. We find that, while many successfully slow down dissociation of the hGH-hGHR complex once bound, they also slow down the association of hGH to hGHR. The E174A substitution induces a change in the Circular Dichroism spectrum that suggests the appearance of coiled-coiling. Here we show that E174A increases affinity of hGH against hGHR because the off-rate is slowed down more than the on-rate. For E174Y (and certain mutations at other sites) the slowdown in on-rate was greater than that of the off-rate, leading to decreased affinity. The results point to a link between structure, zinc binding, and hGHR-binding affinity in hGH.

Multiomics in primary and metastatic breast tumors from the AURORA US network finds microenvironment and epigenetic drivers of metastasis.

The AURORA US Metastasis Project was established with the goal to identify molecular features associated with metastasis. We assayed 55 females with metastatic breast cancer (51 primary cancers and 102 metastases) by RNA sequencing, tumor/germline DNA exome and low-pass whole-genome sequencing and global DNA methylation microarrays. Expression subtype changes were observed in ~30% of samples and were coincident with DNA clonality shifts, especially involving HER2. Downregulation of estrogen receptor (ER)-mediated cell-cell adhesion genes through DNA methylation mechanisms was observed in metastases. Microenvironment differences varied according to tumor subtype; the ER+/luminal subtype had lower fibroblast and endothelial content, while triple-negative breast cancer/basal metastases showed a decrease in B and T cells. In 17% of metastases, DNA hypermethylation and/or focal deletions were identified near HLA-A and were associated with reduced expression and lower immune cell infiltrates, especially in brain and liver metastases. These findings could have implications for treating individuals with metastatic breast cancer with immune- and HER2-targeting therapies.

Association of 17q22 Amplicon Via Cell-Free DNA With Platinum Chemotherapy Response in Metastatic Triple-Negative Breast Cancer.

PURPOSE: To determine whether specific somatic copy-number alterations detectable in circulating tumor DNA (ctDNA) from patients with metastatic triple-negative breast cancer (mTNBC) are associated with sensitivity to platinum chemotherapy. MATERIALS AND METHODS: In this secondary analysis of a large cohort of patients with mTNBC whose ctDNA underwent ultralow-pass whole-genome sequencing, tumor fraction and somatic copy-number alterations were derived with the ichorCNA algorithm. Seventy-two patients were identified who had received a platinum-based chemotherapy regimen in the metastatic setting. Gene-level copy-number analyses were performed with GISTIC2.0. Cytobands were associated with progression-free survival (PFS) to platinum chemotherapy using Cox proportional hazards models. The Cancer Genome Atlas and Molecular Taxonomy of Breast Cancer International Consortium data sets were interrogated for frequency of significant cytobands in primary triple-negative breast cancer (pTNBC) tumors. RESULTS: Among 71 evaluable patients, 17q21 and 17q22 amplifications were most strongly associated with improved PFS with platinum chemotherapy. There were no significant differences in clinicopathologic features or (neo)adjuvant chemotherapy among patients with 17q22 amplification. Patients with 17q22 amplification (n = 17) had longer median PFS with platinum (7.0 v 3.8 months; log-rank P = .015) than patients without 17q22 amplification (n = 54), an effect that remained significant in multivariable analyses (PFS hazard ratio 0.37; 95% CI, 0.16 to 0.84; P = .02). Among 39 patients who received the nonplatinum chemotherapy agent capecitabine, there was no association between 17q22 amplification and capecitabine PFS (log-rank P = .69). In The Cancer Genome Atlas and Molecular Taxonomy of Breast Cancer International Consortium, 17q22 amplification occurred in more than 20% of both pTNBC and mTNBC tumors, whereas 17q21 was more frequently amplified in mTNBC relative to pTNBC (16% v 8.1%, P = .015). CONCLUSION: The 17q22 amplicon, detected by ctDNA, is associated with improved PFS with platinum chemotherapy in patients with mTNBC and warrants further investigation.

A Machine Learning Method to Identify Genetic Variants Potentially Associated With Alzheimer's Disease.

There is hope that genomic information will assist prediction, treatment, and understanding of Alzheimer's disease (AD). Here, using exome data from ∼10,000 individuals, we explore machine learning neural network (NN) methods to estimate the impact of SNPs (i.e., genetic variants) on AD risk. We develop an NN-based method (netSNP) that identifies hundreds of novel potentially protective or at-risk AD-associated SNPs (along with an effect measure); the majority with frequency under 0.01. For case individuals, the number of "protective" (or "at-risk") netSNP-identified SNPs in their genome correlates positively (or inversely) with their age of AD diagnosis and inversely (or positively) with autopsy neuropathology. The effect measure increases correlations. Simulations suggest our results are not due to genetic linkage, overfitting, or bias introduced by netSNP. These findings suggest that netSNP can identify SNPs associated with AD pathophysiology that may assist with the diagnosis and mechanistic understanding of the disease.

EF-P Posttranslational Modification Has Variable Impact on Polyproline Translation in Bacillus subtilis.

Elongation factor P (EF-P) is a ubiquitous translation factor that facilitates translation of polyproline motifs. In order to perform this function, EF-P generally requires posttranslational modification (PTM) on a conserved residue. Although the position of the modification is highly conserved, the structure can vary widely between organisms. In Bacillus subtilis, EF-P is modified at Lys32 with a 5-aminopentanol moiety. Here, we use a forward genetic screen to identify genes involved in 5-aminopentanolylation. Tandem mass spectrometry analysis of the PTM mutant strains indicated that ynbB, gsaB, and ymfI are required for modification and that yaaO, yfkA, and ywlG influence the level of modification. Structural analyses also showed that EF-P can retain unique intermediate modifications, suggesting that 5-aminopentanol is likely directly assembled on EF-P through a novel modification pathway. Phenotypic characterization of these PTM mutants showed that each mutant does not strictly phenocopy the efp mutant, as has previously been observed in other organisms. Rather, each mutant displays phenotypic characteristics consistent with those of either the efp mutant or wild-type B. subtilis depending on the growth condition. In vivo polyproline reporter data indicate that the observed phenotypic differences result from variation in both the severity of polyproline translation defects and altered EF-P context dependence in each mutant. Together, these findings establish a new EF-P PTM pathway and also highlight a unique relationship between EF-P modification and polyproline context dependence.IMPORTANCE Despite the high level of conservation of EF-P, the posttranslational modification pathway that activates EF-P is highly divergent between species. Here, we have identified and characterized in B. subtilis a novel posttranslational modification pathway. This pathway not only broadens the scope of potential EF-P modification strategies, but it also indicates that EF-P modifications can be assembled directly on EF-P. Furthermore, characterization of these PTM mutants has established that an altered modification state can impact both the severity of polyproline translational defects and context dependence.

Elongation Factor P and the Control of Translation Elongation.

Elongation factor P (EF-P) binds to ribosomes requiring assistance with the formation of oligo-prolines. In order for EF-P to associate with paused ribosomes, certain tRNAs with specific d-arm residues must be present in the peptidyl site, e.g., tRNAPro. Once EF-P is accommodated into the ribosome and bound to Pro-tRNAPro, productive synthesis of the peptide bond occurs. The underlying mechanism by which EF-P facilitates this reaction seems to have entropic origins. Maximal activity of EF-P requires a posttranslational modification in Escherichia coli, Pseudomonas aeruginosa, and Bacillus subtilis. Each of these modifications is distinct and ligated onto its respective EF-P through entirely convergent means. Here we review the facets of translation elongation that are controlled by EF-P, with a particular focus on the purpose behind the many different modifications of EF-P.

Carbonyl reduction by YmfI in Bacillus subtilis prevents accumulation of an inhibitory EF-P modification state.

Translation elongation factor P (EF-P) in Bacillus subtilis is required for a form of surface migration called swarming motility. Furthermore, B. subtilis EF-P is post-translationally modified with a 5-aminopentanol group but the pathway necessary for the synthesis and ligation of the modification is unknown. Here we determine that the protein YmfI catalyzes the reduction of EF-P-5 aminopentanone to EF-P-5 aminopentanol. In the absence of YmfI, accumulation of 5-aminopentanonated EF-P is inhibitory to swarming motility. Suppressor mutations that enhanced swarming in the absence of YmfI were found at two positions on EF-P, including one that changed the conserved modification site (Lys 32) and abolished post-translational modification. Thus, while modification of EF-P is thought to be essential for EF-P activity, here we show that in some cases it can be dispensable. YmfI is the first protein identified in the pathway leading to EF-P modification in B. subtilis, and B. subtilis encodes the first EF-P ortholog that retains function in the absence of modification.

Synthesis of rhamnosylated arginine glycopeptides and determination of the glycosidic linkage in bacterial elongation factor P.

A new class of N-linked protein glycosylation - arginine rhamnosylation - has recently been discovered as a critical modification for the function of bacterial elongation factor P (EF-P). Herein, we describe the synthesis of suitably protected α- and ß-rhamnosylated arginine amino acid "cassettes" that can be directly installed into rhamnosylated peptides. Preparation of a proteolytic fragment of Pseudomonas aeruginosa EF-P bearing both α- and ß-rhamnosylated arginine enabled the unequivocal determination of the native glycosidic linkage to be α through 2D NMR and nano-UHPLC-tandem mass spectrometry studies.

Non-canonical roles of tRNAs and tRNA mimics in bacterial cell biology.

Transfer RNAs (tRNAs) are the macromolecules that transfer activated amino acids from aminoacyl-tRNA synthetases to the ribosome, where they are used for the mRNA guided synthesis of proteins. Transfer RNAs are ancient molecules, perhaps even predating the existence of the translation machinery. Albeit old, these molecules are tremendously conserved, a characteristic that is well illustrated by the fact that some bacterial tRNAs are efficient and specific substrates of eukaryotic aminoacyl-tRNA synthetases and ribosomes. Considering their ancient origin and high structural conservation, it is not surprising that tRNAs have been hijacked during evolution for functions outside of translation. These roles beyond translation include synthetic, regulatory and information functions within the cell. Here we provide an overview of the non-canonical roles of tRNAs and their mimics in bacteria, and discuss some of the common themes that arise when comparing these different functions.

Translation Control of Swarming Proficiency in Bacillus subtilis by 5-Amino-pentanolylated Elongation Factor P.

Elongation factor P (EF-P) accelerates diprolyl synthesis and requires a posttranslational modification to maintain proteostasis. Two phylogenetically distinct EF-P modification pathways have been described and are encoded in the majority of Gram-negative bacteria, but neither is present in Gram-positive bacteria. Prior work suggested that the EF-P-encoding gene (efp) primarily supports Bacillus subtilis swarming differentiation, whereas EF-P in Gram-negative bacteria has a more global housekeeping role, prompting our investigation to determine whether EF-P is modified and how it impacts gene expression in motile cells. We identified a 5-aminopentanol moiety attached to Lys(32) of B. subtilis EF-P that is required for swarming motility. A fluorescent in vivo B. subtilis reporter system identified peptide motifs whose efficient synthesis was most dependent on 5-aminopentanol EF-P. Examination of the B. subtilis genome sequence showed that these EF-P-dependent peptide motifs were represented in flagellar genes. Taken together, these data show that, in B. subtilis, a previously uncharacterized posttranslational modification of EF-P can modulate the synthesis of specific diprolyl motifs present in proteins required for swarming motility.

Cyclic Rhamnosylated Elongation Factor P Establishes Antibiotic Resistance in Pseudomonas aeruginosa.

UNLABELLED: Elongation factor P (EF-P) is a ubiquitous bacterial protein that is required for the synthesis of poly-proline motifs during translation. In Escherichia coli and Salmonella enterica, the posttranslational ß-lysylation of Lys34 by the PoxA protein is critical for EF-P activity. PoxA is absent from many bacterial species such as Pseudomonas aeruginosa, prompting a search for alternative EF-P posttranslation modification pathways. Structural analyses of P. aeruginosa EF-P revealed the attachment of a single cyclic rhamnose moiety to an Arg residue at a position equivalent to that at which ß-Lys is attached to E. coli EF-P. Analysis of the genomes of organisms that both lack poxA and encode an Arg32-containing EF-P revealed a highly conserved glycosyltransferase (EarP) encoded at a position adjacent to efp. EF-P proteins isolated from P. aeruginosa ΔearP, or from a ΔrmlC::acc1 strain deficient in dTDP-L-rhamnose biosynthesis, were unmodified. In vitro assays confirmed the ability of EarP to use dTDP-L-rhamnose as a substrate for the posttranslational glycosylation of EF-P. The role of rhamnosylated EF-P in translational control was investigated in P. aeruginosa using a Pro4-green fluorescent protein (Pro4GFP) in vivo reporter assay, and the fluorescence was significantly reduced in Δefp, ΔearP, and ΔrmlC::acc1 strains. ΔrmlC::acc1, ΔearP, and Δefp strains also displayed significant increases in their sensitivities to a range of antibiotics, including ertapenem, polymyxin B, cefotaxim, and piperacillin. Taken together, our findings indicate that posttranslational rhamnosylation of EF-P plays a key role in P. aeruginosa gene expression and survival. IMPORTANCE: Infections with pathogenic Salmonella, E. coli, and Pseudomonas isolates can all lead to infectious disease with potentially fatal sequelae. EF-P proteins contribute to the pathogenicity of the causative agents of these and other diseases by controlling the translation of proteins critical for modulating antibiotic resistance, motility, and other traits that play key roles in establishing virulence. In Salmonella spp. and E. coli, the attachment of ß-Lys is required for EF-P activity, but the proteins required for this posttranslational modification pathway are absent from many organisms. Instead, bacteria such as P. aeruginosa activate EF-P by posttranslational modification with rhamnose, revealing a new role for protein glycosylation that may also prove useful as a target for the development of novel antibiotics.

Elongation factor-P at the crossroads of the host-endosymbiont interface.

Elongation factor P (EF-P) is an ancient bacterial translational factor that aids the ribosome in polymerizing oligo-prolines. EF-P structurally resembles tRNA and binds in-between the exit and peptidyl sites of the ribosome to accelerate the intrinsically slow reaction of peptidyl-prolyl bond formation. Recent studies have identified in separate organisms, two evolutionarily convergent EF-P post-translational modification systems (EPMS), split predominantly between gammaproteobacteria, and betaproteobacteria. In both cases EF-P receives a post-translational modification, critical for its function, on a highly conserved residue that protrudes into the peptidyl-transfer center of the ribosome. EPMSs are comprised of a gene(s) that synthesizes the precursor molecule used in modifying EF-P, and a gene(s) encoding an enzyme that reacts with the precursor molecule to catalyze covalent attachment to EF-P. However, not all organisms genetically encode a complete EPMS. For instance, some symbiotic bacteria harbor efp and the corresponding gene that enzymatically attaches the modification, but lack the ability to synthesize the substrate used in the modification reaction. Here we highlight the recent discoveries made regarding EPMSs, with a focus on how these incomplete modification pathways shape or have been shaped by the endosymbiont-host relationship.

The non-canonical hydroxylase structure of YfcM reveals a metal ion-coordination motif required for EF-P hydroxylation.

EF-P is a bacterial tRNA-mimic protein, which accelerates the ribosome-catalyzed polymerization of poly-prolines. In Escherichia coli, EF-P is post-translationally modified on a conserved lysine residue. The post-translational modification is performed in a two-step reaction involving the addition of a ß-lysine moiety and the subsequent hydroxylation, catalyzed by PoxA and YfcM, respectively. The ß-lysine moiety was previously shown to enhance the rate of poly-proline synthesis, but the role of the hydroxylation is poorly understood. We solved the crystal structure of YfcM and performed functional analyses to determine the hydroxylation mechanism. In addition, YfcM appears to be structurally distinct from any other hydroxylase structures reported so far. The structure of YfcM is similar to that of the ribonuclease YbeY, even though they do not share sequence homology. Furthermore, YfcM has a metal ion-coordinating motif, similar to YbeY. The metal ion-coordinating motif of YfcM resembles a 2-His-1-carboxylate motif, which coordinates an Fe(II) ion and forms the catalytic site of non-heme iron enzymes. Our findings showed that the metal ion-coordinating motif of YfcM plays an essential role in the hydroxylation of the ß-lysylated lysine residue of EF-P. Taken together, our results suggested the potential catalytic mechanism of hydroxylation by YfcM.

Oxidation of cellular amino acid pools leads to cytotoxic mistranslation of the genetic code.

Aminoacyl-tRNA synthetases use a variety of mechanisms to ensure fidelity of the genetic code and ultimately select the correct amino acids to be used in protein synthesis. The physiological necessity of these quality control mechanisms in different environments remains unclear, as the cost vs benefit of accurate protein synthesis is difficult to predict. We show that in Escherichia coli, a non-coded amino acid produced through oxidative damage is a significant threat to the accuracy of protein synthesis and must be cleared by phenylalanine-tRNA synthetase in order to prevent cellular toxicity caused by mis-synthesized proteins. These findings demonstrate how stress can lead to the accumulation of non-canonical amino acids that must be excluded from the proteome in order to maintain cellular viability.

(R)-ß-lysine-modified elongation factor P functions in translation elongation.

Post-translational modification of bacterial elongation factor P (EF-P) with (R)-ß-lysine at a conserved lysine residue activates the protein in vivo and increases puromycin reactivity of the ribosome in vitro. The additional hydroxylation of EF-P at the same lysine residue by the YfcM protein has also recently been described. The roles of modified and unmodified EF-P during different steps in translation, and how this correlates to its physiological role in the cell, have recently been linked to the synthesis of polyproline stretches in proteins. Polysome analysis indicated that EF-P functions in translation elongation, rather than initiation as proposed previously. This was further supported by the inability of EF-P to enhance the rate of formation of fMet-Lys or fMet-Phe, indicating that the role of EF-P is not to specifically stimulate formation of the first peptide bond. Investigation of hydroxyl-(ß)-lysyl-EF-P showed 30% increased puromycin reactivity but no differences in dipeptide synthesis rates when compared with the ß-lysylated form. Unlike disruption of the other genes required for EF-P modification, deletion of yfcM had no phenotypic consequences in Salmonella. Taken together, our findings indicate that EF-P functions in translation elongation, a role critically dependent on post-translational ß-lysylation but not hydroxylation.

Identification and characterization of a putative arginine kinase homolog from Myxococcus xanthus required for fruiting body formation and cell differentiation.

Arginine kinases catalyze the reversible transfer of a high-energy phosphoryl group from ATP to l-arginine to form phosphoarginine, which is used as an energy buffer in insects, crustaceans, and some unicellular organisms. It plays an analogous role to that of phosphocreatine in vertebrates. Recently, putative arginine kinases were identified in several bacterial species, including the social Gram-negative soil bacterium Myxococcus xanthus. It is still unclear what role these proteins play in bacteria and whether they have evolved to acquire novel functions in the species in which they are found. In this study, we biochemically purified and characterized a putative M. xanthus arginine kinase, Ark, and demonstrated that it has retained the ability to catalyze the phosphorylation of arginine by using ATP. We also constructed a null mutation in the ark gene and demonstrated its role in both certain stress responses and development.