Search

Background: Mitochondrial dysfunction is a hallmark of cancer. However, it is unclear whether it is a cause of cancer. This two-sample Mendelian randomization (MR) analyses, uses genetic instruments to proxy the exposure of mitochondrial dysfunction and cancer summary statistics as outcomes, allowing for causal inferences., Methods: Summary statistics from 18 common cancers (2107-491,974 participants), gene expression, DNA methylation and protein expression quantitative trait loci (eQTL, mQTL and pQTL, respectively, 1000-31,684 participants) on individuals of European ancestry, were included. Genetic variants located within or close to the 1136 mitochondrial-related genes (in cis) and robustly associated with the mitochondrial molecular alterations were used as instrumental variables, and their causal associations with cancers were examined using summary-data-based MR (SMR) analyses. An additional five MR methods were used as sensitivity analyses to confirm the casual associations. A Bayesian test for colocalization between mitochondrial molecular QTLs and cancer risk loci was performed to provide insights into the potential regulatory mechanisms of risk variants on cancers., Findings: We identified potential causal relationships between mitochondrial-related genes and breast, prostate, gastric, lung cancer and melanoma by primary SMR analyses. The sensitivity and the colocalization analyses further refined four genes that have causal effects on three types of cancer. We found strong evidence of positive association of FDPS expression level with breast cancer risk (OR per SD, 0.66; 95% CI, 0.49-0.83; P = 9.77 × 10 -7 ), NSUN4 expression level with both breast cancer risk (OR per SD, 1.05; 95% CI, 1.03-1.07; P = 5.24 × 10 -6 ) and prostate cancer risk (OR per SD, 1.06; 95% CI, 1.03-1.09; P = 1.01 × 10 -5 ), NSUN4 methylation level with both breast and prostate cancer risk, and VARS2 methylation level with lung cancer risk., Interpretations: This data-driven MR study demonstrated the causal role of mitochondrial dysfunction in multiple cancers. Furthermore, this study identified candidate genes that can be the targets of potential pharmacological agents for cancer prevention., Funding: This work was supported by Styrelsen för Allmänna Sjukhusets i Malmö Stiftelse för bekämpande av cancer (20211025)., Competing Interests: Declaration of interests The authors declare that there is no conflict of interest associated with this manuscript., (Copyright © 2022 The Author(s). Published by Elsevier B.V. All rights reserved.)

Environmental factors contribute to risk of bipolar disorder (BD), but how environmental factors impact the development of psychopathology within the context of elevated genetic risk is unknown. We herein sought to identify epigenetic signatures operating in the context of polygenic risk for BD in young people at high familial risk (HR) of BD. Peripheral blood-derived DNA was assayed using Illumina PsychArray, and Methylation-450K or -EPIC BeadChips. Polygenic risk scores (PRS) were calculated using summary statistics from recent genome-wide association studies for BD, major depressive disorder (MDD) and cross-disorder (meta-analysis of eight psychiatric disorders). Unrelated HR participants of European ancestry (n = 103) were stratified based on their BD-PRS score within the HR-population distribution, and the top two quintiles (High-BD-PRS; n = 41) compared against the bottom two quintiles (Low-BD-PRS; n = 41). The High-BD-PRS stratum also had higher mean cross-disorder-PRS and MDD-PRS (ANCOVA p = 0.035 and p = 0.024, respectively). We evaluated DNA methylation differences between High-BD-PRS and Low-BD-PRS strata using linear models. One differentially methylated probe (DMP) (cg00933603; p = 3.54 × 10 -7 ) in VARS2, a mitochondrial aminoacyl-tRNA synthetase, remained significantly hypomethylated after multiple-testing correction. Overall, BD-PRS appeared to broadly impact epigenetic processes, with 1,183 genes mapped to nominal DMPs (p < 0.05); these displayed convergence with genes previously associated with BD, schizophrenia, chronotype, and risk taking. We tested poly-methylomic epigenetic profiles derived from nominal DMPs in two independent samples (n = 54 and n = 82, respectively), and conducted an exploratory evaluation of the effects of family environment, indexing cohesion and flexibility. This study highlights an important interplay between heritable risk and epigenetic factors, which warrant further exploration., (© 2022. The Author(s).)

Mutations in mitochondrial aminoacyl-tRNA synthetases (mtARSs) have been reported in patients with mitochondriopathies: most commonly encephalopathy, but also cardiomyopathy. Through a GWAS, we showed possible associations between mitochondrial valyl-tRNA synthetase (VARS2) dysregulations and non-ischemic cardiomyopathy. We aimed to investigate the possible consequences of VARS2 depletion in zebrafish and cultured HEK293A cells. Transient VARS2 loss-of-function was induced in zebrafish embryos using Morpholinos. The enzymatic activity of VARS2 was measured in VARS2-depleted cells via northern blot. Heterozygous VARS2 knockout was established in HEK293A cells using CRISPR/Cas9 technology. BN-PAGE and SDS-PAGE were used to investigate electron transport chain (ETC) complexes, and the oxygen consumption rate and extracellular acidification rate were measured using a Seahorse XFe96 Analyzer. The activation of the integrated stress response (ISR) and possible disruptions in mitochondrial fatty acid oxidation (FAO) were explored using RT-qPCR and western blot. Zebrafish embryos with transient VARS2 loss-of-function showed features of heart failure as well as indications of CNS and skeletal muscle involvements. The enzymatic activity of VARS2 was significantly reduced in VARS2-depleted cells. Heterozygous VARS2-knockout cells showed a rearrangement of ETC complexes in favor of complexes III 2 , III 2 + IV, and supercomplexes without significant respiratory chain deficiencies. These cells also showed the enhanced activation of the ISR, as indicated by increased eIF-2α phosphorylation and a significant increase in the transcript levels of ATF4, ATF5, and DDIT3 (CHOP), as well as disruptions in FAO. The activation of the ISR and disruptions in mitochondrial FAO may underlie the adaptive changes in VARS2-depleted cells.

Isoleucyl-tRNA synthetase (IleRS), leucyl-tRNA synthetase (LeuRS), and valyl-tRNA synthetase (ValRS) are enzymes that have potential for the determination of l-isoleucine, l-leucine, and l-valine in food products and plasma. However, the disadvantages of these enzymes are their specificity and sensitivity. Here, we examined the substrate specificity of IleRS, LeuRS, and ValRS under various conditions of pyrophosphate amplification to improve their specificity and sensitivity. The amount of pyrophosphate produced in IleRS, LeuRS, and ValRS reactions was amplified after the addition of excess adenosine-5'-triphosphate and magnesium ions, and was approximately 9-, 8-, and 7-fold higher, respectively, for each of the initial l-amino acid substrates (50 μM). However, in addition to their target amino acids, IleRS, LeuRS, and ValRS also reacted with l-valine, l-lysine, and l-threonine, respectively. This substrate misrecognition was overcome by making the reaction pH more acidic and by increasing the magnesium ion concentration. The pyrophosphate amplification in IleRS, LeuRS, and ValRS reactions resulted in the production of p 1 , p 4 -di (adenosine) 5'-tetraphosphate. We also observed a strong positive correlation (R = 0.99) between the amount of pyrophosphate produced and the initial concentration of l-amino acid with 5 and 50 μM l-isoleucine, l-leucine, and l-valine. Our results suggest that amino acid assays using IleRS, LeuRS, and ValRS are promising methods to accurately measure l-valine, l-isoleucine, and l-leucine in food products and plasma., (Copyright © 2022 The Society for Biotechnology, Japan. Published by Elsevier B.V. All rights reserved.)

Single mutations can confer resistance to antibiotics. Identifying such mutations can help to develop and improve drugs. Here, we systematically screen for candidate quinolone resistance-conferring mutations. We sequenced highly diverse wastewater E. coli and performed a genome-wide association study (GWAS) to determine associations between over 200,000 mutations and quinolone resistance phenotypes. We uncovered 13 statistically significant mutations including 1 located at the active site of the biofilm dispersal gene bdcA and 6 silent mutations in the aminoacyl-tRNA synthetase valS . The study also recovered the known mutations in the topoisomerases gyrase ( gyrA ) and topoisomerase IV ( parC ). In summary, we demonstrate that GWAS effectively and comprehensively identifies resistance mutations without a priori knowledge of targets and mode of action. The results suggest that mutations in the bdcA and valS genes, which are involved in biofilm dispersal and translation, may lead to novel resistance mechanisms.

Isoleucyl-tRNA synthetase (IleRS) is a paradigm for understanding how specificity against smaller hydrophobic substrates evolved in both the synthetic and editing reactions. IleRS misactivates nonproteinogenic norvaline (Nva) and proteinogenic valine (Val), with a 200-fold lower efficiency than the cognate isoleucine (Ile). Translational errors are, however, prevented by IleRS hydrolytic editing. Nva and Val are both smaller than Ile by a single methylene group. How does the removal of one additional methylene group affects IleRS specificity? We found that the nonproteinogenic α-aminobutyrate (Abu) is activated 30-fold less efficiently than Nva and Val, indicating that the removal of the second methylene group comes with a lower penalty. As with Nva and Val, discrimination against Abu predominantly originated from a higher K M . To examine whether increased hydrophobicity could compensate for the loss of van der Waals interactions, we tested fluorinated Abu analogues. We found that fluorination further hampered activation by IleRS, and even more so by the evolutionary-related ValRS. This suggests that hydrophobicity is not a main driving force of substrate binding in these enzymes. Finally, a discrimination factor of 7100 suggests that IleRS is not expected to edit Abu. However, we found that the IleRS editing domain hydrolyzes Abu-tRNA Ile with a rate of 40 s -1 and the introduction of fluorine did not slow down the hydrolysis. This raises interesting questions regarding the mechanism of specificity of the editing domain and its evolution. Understanding what shapes IleRS specificity is also of importance for reengineering translation to accommodate artificial substrates including fluorinated amino acids. ENZYMES: Isoleucyl-tRNA synthetase (EC 6.1.1.5), leucyl-tRNA synthetase (EC 6.1.1.4), valyl-tRNA synthetase (EC 6.1.1.9)., (© 2019 Federation of European Biochemical Societies.)

Consanguineous populations of the Arabian Peninsula, which has seen an uncontrolled rise in type 2 diabetes incidence, are underrepresented in global studies on diabetes genetics. We performed a genome-wide association study on the quantitative trait of fasting plasma glucose (FPG) in unrelated Arab individuals from Kuwait (discovery-cohort:n = 1,353; replication-cohort:n = 1,196). Genome-wide genotyping in discovery phase was performed for 632,375 markers from Illumina HumanOmniExpress Beadchip; and top-associating markers were replicated using candidate genotyping. Genetic models based on additive and recessive transmission modes were used in statistical tests for associations in discovery phase, replication phase, and meta-analysis that combines data from both the phases. A genome-wide significant association with high FPG was found at rs1002487 (RPS6KA1) (p-discovery = 1.64E-08, p-replication = 3.71E-04, p-combined = 5.72E-11; β-discovery = 8.315; β-replication = 3.442; β-combined = 6.551). Further, three suggestive associations (p-values < 8.2E-06) with high FPG were observed at rs487321 (CADPS), rs707927 (VARS and 2Kb upstream of VWA7), and rs12600570 (DHX58); the first two markers reached genome-wide significance in the combined analysis (p-combined = 1.83E-12 and 3.07E-09, respectively). Significant interactions of diabetes traits (serum triglycerides, FPG, and glycated hemoglobin) with homeostatic model assessment of insulin resistance were identified for genotypes heterozygous or homozygous for the risk allele. Literature reports support the involvement of these gene loci in type 2 diabetes etiology.

Pneumonia is a common and severe infectious lung disease. Host genetics, together with underlying medical and lifestyle conditions, determine pneumonia susceptibility. We performed a secondary analysis of the results of two genome-wide studies for pneumonia in 23andMe participants (40 600 cases/90 039 controls) (Tian et al., 2017) and UK Biobank (BB) participants (12 614 cases/324 585 controls) (via the Global Biobank Engine) and used the GTEx database to correlate the results with expression quantitative trait loci (eQTLs) data in lung and whole blood. In the 23andMe pneumonia single nucleotide polymorphism (SNP) set, 177 genotyped SNPs in the human leukocyte antigen (HLA) region satisfied the genome-wide significance level, P ≤ 5·0E-08. Several target genes (e.g. C4A, VARS2, SFTA2, HLA-C, HLA-DQA2) were unidirectionally regulated by many HLA eSNPs associated with a higher risk of pneumonia. In lung, C4A transcript was up-regulated by 291 pneumonia risk alleles spanning the half the HLA region. Among SNPs correlated with the expression levels of SFTA2 and VARS2, approximately 75% overlapped: all risk alleles were associated with VARS2 up-regulation and SFTA2 down-regulation. To find shared gene loci between pneumonia and pulmonary function (PF), we used data from the Global Biobank Engine and literature on genome-wide association studies (GWAS) of PF in general populations. Numerous gene loci overlapped between pneumonia and PF: 28·8% in the BB data set and 49·2% in the 23andMe data set. Enrichment analysis within the database of Genotypes and Phenotypes (dbGaP) and National Human Genome Research Institute-European Bioinformatics Institute (NHGRI-EBI) Catalog of pneumonia and pneumonia/PF gene sets identified significant overlap between these gene sets and genes related to inflammatory, developmental, neuropsychiatric and cardiovascular and obesity-related traits., (© 2019 British Society for Immunology.)

The VARS2 gene encodes a mitochondrial valyl-transfer RNA synthetase which is used in mitochondrial translation. To date, several patients with VARS2 pathogenic variants have been described in the literature. These patients have features of lactic acidosis with encephalomyopathy. We present a case of an infant with lactic acidosis, failure to thrive, and severe primary pulmonary hypertension who was found to be a compound heterozygote for two novel VARS2 variants (c.1940C>T, p.(Thr647Met) and c.2318G>A, p.(Arg773Gln)). The patient was treated with vitamin supplements and a carbohydrate-restricted diet. The lactic acidosis and failure to thrive resolved, and he showed good growth and development. Functional studies and molecular analysis employed a yeast model system and the VAS1 gene (yeast homolog of VARS2). VAS1 genes harboring either one of two mutations corresponding to the two novel variants in the VARS2 gene, exhibited partially reduced function in haploid yeast strains. A combination of both VAS1 variant alleles in a diploid yeast cell exhibited a more significant decrease in oxidative metabolism-dependent growth and in the oxygen consumption rate (reminiscent of the patient who carries two mutant VARS2 alleles). Our results demonstrate the pathogenicity of the biallellic novel VARS2 variants. KEY MESSAGES: • A case of an infant who is a compound heterozygote for two novel VARS2 variants. • This infant displayed lactic acidosis, failure to thrive, and pulmonary hypertension. • Treatment of the patient with a carbohydrate-restricted diet resulted in good growth and development. • Studies with the homologous yeast VAS1 gene showed reduced function of corresponding single mutant in haploid yeast strains. • A combination of both VAS1 variant alleles in diploid yeast exhibited a more significant decrease in function, thereby confirming the pathogenicity of the biallellic novel VARS2 variants.

Background: Mitochondrial disorders are heterogeneous clinical syndromes caused by defective activity in the mitochondrial respiratory chain, resulting in a faulty oxidative phosphorylation system. These inherited disorders are individually rare, and furthermore they are phenotypic variables. The genetically characterized mitochondrial disorders are rarely associated with epileptic encephalopathies. Case presentation: We present the clinical phenotype, biochemical analysis, and electrographic and neuro-radiological features of a 5-month-old girl with epileptic encephalopathy, microcephaly, severe psychomotor delay, hypertrophic cardiomyopathy, and abnormal MRI scan. Using whole-genome sequencing technique, compound heterozygous mutations of the VARS2 gene were revealed, with one previously unreported frameshift mutation. Conclusion: Our report extends the phenotypic spectrum of VARS2-related disorders with an initial presentation of epileptic encephalopathy and early death due to malignant arrhythmia.

Background: Mitochondrial respiratory chain consists of five complexes encoded by nuclear and mitochondrial genomes. Mitochondrial aminoacyl-tRNA synthetases are key enzymes in the synthesis of such complexes. Bi-allelic variants of VARS2, a nuclear gene encoding for valyl-tRNA (Val-tRNA) synthetase, are associated to several forms of mitochondrial encephalopathies or cardiomyoencephalopathies. Among these, the rare homozygous c.1100C > T (p.Thr367Ile) mutation variably presents with progressive developmental delay, axial hypotonia, limbs spasticity, drug-resistant epilepsy leading, in some cases, to premature death. Yet only six cases, of which three are siblings, harbouring this homozygous mutation have been described worldwide., Case Presentation: Hereby, we report two additional cases of two non-related young girls from Sardinia, born from non-consanguineous and healthy parents, carrying the aforesaid homozygous VARS2 variant. At onset both the patients presented with worsening psychomotor delay, muscle hypotonia and brisk tendon reflexes. Standard genetic tests were normal, as well as metabolic investigations. Brain MRI showed unspecific progressive abnormalities, such as corpus callosum hypoplasia (patient A) and cerebellar atrophy (patient A and B). Diagnosis was reached by adopting massive parallel next generation sequencing. Notably clinical phenotype of the first patient appears to be milder compared to previous known cases. The second patient eventually developed refractory epilepsy and currently presents with severe global impairment. Because no specific treatment is available as yet, both patients are treated with supporting antioxidant compounds along with symptomatic therapies., Conclusions: Given the paucity of clinical data about this very rare mitochondrial encephalopathy, our report might contribute to broaden the phenotypic spectrum of the disorder. Moreover, noteworthy, three out of five pedigrees so far described belong to the Northern Sardinia ethnicity.

In this study, we report the metabolic consequences of the m.1630 A > G variant in fibroblasts from the symptomatic proband affected with the mitochondrial encephalomyopathy lactic acidosis and stroke-like episode Syndrome and her asymptomatic mother. By long-range PCR followed by massively parallel sequencing of the mitochondrial genome, we accurately measured heteroplasmy in fibroblasts from the proband (89.6%) and her mother (94.8%). Using complementary experimental approaches, we show a functional correlation between manifestation of clinical symptoms and bioenergetic potential. Our mitochondrial morphometric analysis reveals a link between defects of mitochondrial cristae ultrastructure and symptomatic status. Despite near-homoplasmic level of the m.1630A > G variant, the mother's fibroblasts have a normal OXPHOS metabolism, which stands in contrast to the severely impaired OXPHOS response of the proband's fibroblasts. The proband's fibroblasts also exhibit glycolysis at near constitutive levels resulting in a stunted compensatory glycolytic response to offset the severe OXPHOS defect. Whole exome sequencing reveals the presence of a heterozygous nonsense VARS2 variant (p.R334X) exclusively in the proband, which removes two thirds of the VARS2 protein containing key domains interacting with the mt-tRNA val and may play a role in modulating the penetrance of the m.1630A > G variant despite similar near homoplasmic levels. Our transmission electron microscopy study also shows unexpected ultrastructural changes of chromatin suggestive of differential epigenomic regulation between the proband and her mother that may explain the differential OXPHOS response between the proband and her mother. Future study will decipher by which molecular mechanisms the nuclear background influences the penetrance of the m.1630 A > G variant causing MELAS., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Aminoacyl-tRNA synthetases (ARSs) function to transfer amino acids to cognate tRNA molecules, which are required for protein translation. To date, biallelic mutations in 31 ARS genes are known to cause recessive, early-onset severe multi-organ diseases. VARS encodes the only known valine cytoplasmic-localized aminoacyl-tRNA synthetase. Here, we report seven patients from five unrelated families with five different biallelic missense variants in VARS. Subjects present with a range of global developmental delay, epileptic encephalopathy and primary or progressive microcephaly. Longitudinal assessment demonstrates progressive cortical atrophy and white matter volume loss. Variants map to the VARS tRNA binding domain and adjacent to the anticodon domain, and disrupt highly conserved residues. Patient primary cells show intact VARS protein but reduced enzymatic activity, suggesting partial loss of function. The implication of VARS in pediatric neurodegeneration broadens the spectrum of human diseases due to mutations in tRNA synthetase genes.

Aminoacyl tRNA synthetases (ARSs) link specific amino acids with their cognate transfer RNAs in a critical early step of protein translation. Mutations in ARSs have emerged as a cause of recessive, often complex neurological disease traits. Here we report an allelic series consisting of seven novel and two previously reported biallelic variants in valyl-tRNA synthetase (VARS) in ten patients with a developmental encephalopathy with microcephaly, often associated with early-onset epilepsy. In silico, in vitro, and yeast complementation assays demonstrate that the underlying pathomechanism of these mutations is most likely a loss of protein function. Zebrafish modeling accurately recapitulated some of the key neurological disease traits. These results provide both genetic and biological insights into neurodevelopmental disease and pave the way for further in-depth research on ARS related recessive disorders and precision therapies.

The multi-subunit eEF1 complex plays a crucial role in de novo protein synthesis. The central functional component of the complex is eEF1A, which occurs as two independently encoded variants with reciprocal expression patterns: whilst eEF1A1 is widely expressed, eEF1A2 is found only in neurons and muscle. Heterozygous mutations in the gene encoding eEF1A2, EEF1A2, have recently been shown to cause epilepsy, autism, and intellectual disability. The remaining subunits of the eEF1 complex, eEF1Bα, eEF1Bδ, eEF1Bγ, and valyl-tRNA synthetase (VARS), together form the GTP exchange factor for eEF1A and are ubiquitously expressed, in keeping with their housekeeping role. However, mutations in the genes encoding these subunits EEF1B2 (eEF1Bα), EEF1D (eEF1Bδ), and VARS (valyl-tRNA synthetase) have also now been identified as causes of neurodevelopmental disorders. In this review, we describe the mutations identified so far in comparison with the degree of normal variation in each gene, and the predicted consequences of the mutations on the functions of the proteins and their isoforms. We discuss the likely effects of the mutations in the context of the role of protein synthesis in neuronal development., (© 2018 Crown copyright. Human Mutation © 2018 Wiley Periodicals, Inc.)

Background: Genetic defects in the mitochondrial aminoacyl-tRNA synthetase are important causes of mitochondrial disorders. VARS2 is one of the genes encoding aminoacyl-tRNA synthetases. Recently, an increasing number of pathogenic variants of VARS2 have been reported., Case Presentation: We report the novel compound heterozygous pathogenic VARS2 mutations c.643 C > T (p. His215Tyr) and c.1354 A > G (p. Met452Val) in a female infant who presented with poor sucking at birth, poor activity, hyporeflexia, hypertonia, persistent pulmonary hypertension of newborn (PPHN), metabolic acidosis, severe lactic acidosis, expansion and hypertrophic cardiomyopathy. These heterozygous mutations were carried individually by the proband's parents and elder sister; the two mutations segregated in the family and were the cause of the disease in the proband.The c.643 C > T (p. His215Tyr) mutation was not described in the ExaC, GNomAD and 1000 Genomes Project databases, and the frequency of c.1354 A > G (p. Met452Val) was < 0.001 in these gene databases.The two mutated amino acids were located in a highly conserved region of the VARS2 protein that is important for its interaction with the cognate tRNA. The two missense mutations were predicted by online tools to be damaging and deleterious., Conclusions: Our report expands the spectrum of known pathogenicVARS2 variants associated with mitochondrial disorders in humans.VARS2 deficiency may cause a severe neonatal presentation with structural cardiac abnormalities.

Two male siblings ages 15 and 10 yr old had similar features of intellectual disability, developmental delay, severe speech impairment, microcephaly, prematurity, and transient elevation of liver enzymes in infancy. Exome sequencing revealed one novel (c.65C>A; p.Ala22Asp) and one ultra-rare (c.3214T>C; p.Phe1072Leu) predicted damaging missense variant in trans in the gene encoding cytoplasmic valyl-tRNA synthetase ( VARS ). Biallelic variants in VARS have previously been associated with a neurodevelopmental disorder characterized by microcephaly, seizures, and cortical atrophy (NDMSCA; MIM #617802). Although our patients have no history of seizures or cortical atrophy, we suggest that the biallelic variants in VARS p.Ala22Asp and p.Phe1072Leu in this family are likely pathogenic and associated with NDMSCA, expanding the clinical phenotype of the condition., (© 2018 Okur et al.; Published by Cold Spring Harbor Laboratory Press.)

Progressive microcephaly and neurodegeneration are genetically heterogenous conditions, largely associated with genes that are essential for the survival of neurons. In this study, we interrogate the genetic etiology of two siblings from a non-consanguineous family with severe early onset of neurological manifestations. Whole exome sequencing identified novel compound heterozygous mutations in VARS that segregated with the proband: a missense (c.3192G>A; p.Met1064Ile) and a splice site mutation (c.1577-2A>G). The VARS gene encodes cytoplasmic valyl-tRNA synthetase (ValRS), an enzyme that is essential during eukaryotic translation. cDNA analysis on patient derived fibroblasts revealed that the splice site acceptor variant allele led to nonsense mediated decay, thus resulting in a null allele. Three-dimensional modeling of ValRS predicts that the missense mutation lies in a highly conserved region and could alter side chain packing, thus affecting tRNA binding or destabilizing the interface between the catalytic and tRNA binding domains. Further quantitation of the expression of VARS showed remarkably reduced levels of mRNA and protein in skin derived fibroblasts. Aminoacylation experiments on patient derived cells showed markedly reduced enzyme activity of ValRS suggesting the mutations to be loss of function. Bi-allelic mutations in cytoplasmic amino acyl tRNA synthetases are well-known for their role in neurodegenerative disorders, yet human disorders associated with VARS mutations have not yet been clinically well characterized. Our study describes the phenotype associated with recessive VARS mutations and further functional delineation of the pathogenicity of novel variants identified, which widens the clinical and genetic spectrum of patients with progressive microcephaly.

In recent years, an increasing number of mitochondrial disorders have been associated with mutations in mitochondrial aminoacyl-tRNA synthetases (mt-aaRSs), which are key enzymes of mitochondrial protein synthesis. Bi-allelic functional variants in VARS2, encoding the mitochondrial valyl tRNA-synthetase, were first reported in a patient with psychomotor delay and epilepsia partialis continua associated with an oxidative phosphorylation (OXPHOS) Complex I defect, before being described in a patient with a neonatal form of encephalocardiomyopathy. Here we provide a detailed genetic, clinical, and biochemical description of 13 patients, from nine unrelated families, harboring VARS2 mutations. All patients except one, who manifested with a less severe disease course, presented at birth exhibiting severe encephalomyopathy and cardiomyopathy. Features included hypotonia, psychomotor delay, seizures, feeding difficulty, abnormal cranial MRI, and elevated lactate. The biochemical phenotype comprised a combined Complex I and Complex IV OXPHOS defect in muscle, with patient fibroblasts displaying normal OXPHOS activity. Homology modeling supported the pathogenicity of VARS2 missense variants. The detailed description of this cohort further delineates our understanding of the clinical presentation associated with pathogenic VARS2 variants and we recommend that this gene should be considered in early-onset mitochondrial encephalomyopathies or encephalocardiomyopathies., (© 2018 The Authors. Human Mutation published by Wiley Periodicals, Inc.)

Background: Most mitochondrial and cytoplasmic aminoacyl-tRNA synthetases (aaRSs) are encoded by nuclear genes. Syndromic disorders resulting from mutation of aaRSs genes display significant phenotypic heterogeneity. We expand aaRSs-related phenotypes through characterization of the clinical and molecular basis of a novel autosomal-recessive syndrome manifesting severe mental retardation, ataxia, speech impairment, epilepsy, short stature, microcephaly, hypogonadism, and growth hormone deficiency., Results: A G>A variant in exon 29 of VARS2 (c.3650G>A) (NM&#95;006295) was identified in the index case. This homozygous variant was confirmed by Sanger sequencing and segregated with disease in the family studied. The c.3650G>A change results in alteration of arginine to histidine at residue 1217 (R1217H) of the mature protein and is predicted to be pathogenic., Conclusions: These findings contribute to a growing list of aaRSs disorders, broadens the spectrum of phenotypes attributable to VARS2 mutations, and provides new insight into genotype-phenotype correlations among the mitochondrial synthetase genes.

Sequence-based association studies are at a critical inflexion point with the increasing availability of exome-sequencing data. A popular test of association is the sequence kernel association test (SKAT). Weights are embedded within SKAT to reflect the hypothesized contribution of the variants to the trait variance. Because the true weights are generally unknown, and so are subject to misspecification, we examined the efficiency of a data-driven weighting scheme. We propose the use of a set of theoretically defensible weighting schemes, of which, we assume, the one that gives the largest test statistic is likely to capture best the allele frequency-functional effect relationship. We show that the use of alternative weights obviates the need to impose arbitrary frequency thresholds. As both the score test and the likelihood ratio test (LRT) may be used in this context, and may differ in power, we characterize the behavior of both tests. The two tests have equal power, if the weights in the set included weights resembling the correct ones. However, if the weights are badly specified, the LRT shows superior power (due to its robustness to misspecification). With this data-driven weighting procedure the LRT detected significant signal in genes located in regions already confirmed as associated with schizophrenia - the PRRC2A (p = 1.020e-06) and the VARS2 (p = 2.383e-06) - in the Swedish schizophrenia case-control cohort of 11,040 individuals with exome-sequencing data. The score test is currently preferred for its computational efficiency and power. Indeed, assuming correct specification, in some circumstances, the score test is the most powerful test. However, LRT has the advantageous properties of being generally more robust and more powerful under weight misspecification. This is an important result given that, arguably, misspecified models are likely to be the rule rather than the exception in weighting-based approaches.

VARS2 encodes a mitochondrial aminoacyl-tRNA-synthetase. Mutations in VARS2 have recently been identified as a cause of mitochondrial encephalomyopathy in three individuals. However, clinical information remained scarce. Exome sequencing lead us to identify compound heterozygous pathogenic VARS2 variants in a boy presenting with severe lactic acidosis, hypertrophic cardiomyopathy, epilepsy, and abnormalities on brain imaging including hypoplasia of corpus callosum and cerebellum as well as a massive lactate peak on MR-spectroscopy. Studies in patient-derived fibroblasts confirmed the functional relevance of the identified VARS2 variants. Our report expands the phenotypic spectrum associated with this rare mitochondrial defect, in that VARS2 deficiency may also cause severe neonatal presentations with cardiac involvement and structural brain abnormalities.

Chloroplasts and mitochondria contain their own genomes and transcriptional and translational systems. Establishing these genetic systems is essential for plant growth and development. Here we characterized a mutant form of a Val-tRNA synthetase (OsValRS2) from Oryza sativa that is targeted to both chloroplasts and mitochondria. A single base change in OsValRS2 caused virescent to albino phenotypes in seedlings and white panicles at heading. We therefore named this mutant white panicle 1 (wp1). Chlorophyll autofluorescence observations and transmission electron microscopy analyses indicated that wp1 mutants are defective in early chloroplast development. RNA-seq analysis revealed that expression of nuclear-encoded photosynthetic genes is significantly repressed, while expression of many chloroplast-encoded genes also changed significantly in wp1 mutants. Western-blot analyses of chloroplast-encoded proteins showed that chloroplast protein levels were reduced in wp1 mutants, although mRNA levels of some genes were higher in wp1 than in wild type. We found that wp1 was impaired in chloroplast ribosome biogenesis. Taken together, our results show that OsValRS2 plays an essential role in chloroplast development and regulating chloroplast ribosome biogenesis., (© 2016 American Society of Plant Biologists. All Rights Reserved.)

Background & Aims: Hepatitis B virus (HBV) infection is the most serious risk factor for chronic hepatitis B (CHB), cirrhosis, and hepatocellular carcinoma. Recently, several genome-wide association studies (GWASs) identified important variants associated with the risk of CHB in Asian populations. Specifically, our previous GWAS identified the VARS2-SFTA2 gene region as one of the genetic risk loci for CHB., Methods: To further characterize this association and to isolate possible causal variants within it, we performed an additional association study by genotyping more SNPs in the vicinity of the VARS2 and SFTA2 genes. In all, 14 SNPs of VARS2-SFTA2 were analysed among a total of 3902 subjects (1046 cases and 2856 controls)., Results: Logistic regression analysis revealed that six SNPs, including the previously reported rs2532932, were significantly associated with the risk of CHB (P = 1.7 × 10(-10) ~0.002). Further linkage disequilibrium and conditional analysis identified two variants (rs9394021 and rs2517459) as new markers of genetic risk factors for CHB rather than the reported SNP from our previous study (rs2532932). To evaluate the cumulative risk for CHB based on all known genetic factors, genetic risk score (GRS) were calculated. As anticipated, the distribution of the number of risk alleles in cases vs. controls clearly differed according to the GRS. Similarly, the odds ratios (ORs) were increased (OR = 0.32-3.97)., Conclusion: Our findings show that common variants in the VARS2-SFTA2 gene region are significantly associated with CHB in a Korean population, which may be useful in further understanding genetic susceptibility to CHB., (© 2014 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.)

Bacterial ribosomes stall on polyproline stretches and require the elongation factor P (EF-P) to relieve the arrest. Yet it remains unclear why evolution has favored the development of EF-P rather than selecting against the occurrence of polyproline stretches in proteins. We have discovered that only a single polyproline stretch is invariant across all domains of life, namely a proline triplet in ValS, the tRNA synthetase, that charges tRNA(Val) with valine. Here, we show that expression of ValS in vivo and in vitro requires EF-P and demonstrate that the proline triplet located in the active site of ValS is important for efficient charging of tRNA(Val) with valine and preventing formation of mischarged Thr-tRNA(Val) as well as efficient growth of E. coli in vivo. We suggest that the critical role of the proline triplet for ValS activity may explain why bacterial cells coevolved the EF-P rescue system.

By way of whole-exome sequencing, we identified a homozygous missense mutation in VARS2 in one subject with microcephaly and epilepsy associated with isolated deficiency of the mitochondrial respiratory chain (MRC) complex I and compound heterozygous mutations in TARS2 in two siblings presenting with axial hypotonia and severe psychomotor delay associated with multiple MRC defects. The nucleotide variants segregated within the families, were absent in Single Nucleotide Polymorphism (SNP) databases and are predicted to be deleterious. The amount of VARS2 and TARS2 proteins and valyl-tRNA and threonyl-tRNA levels were decreased in samples of afflicted patients according to the genetic defect. Expression of the corresponding wild-type transcripts in immortalized mutant fibroblasts rescued the biochemical impairment of mitochondrial respiration and yeast modeling of the VARS2 mutation confirmed its pathogenic role. Taken together, these data demonstrate the role of the identified mutations for these mitochondriopathies. Our study reports the first mutations in the VARS2 and TARS2 genes, which encode two mitochondrial aminoacyl-tRNA synthetases, as causes of clinically distinct, early-onset mitochondrial encephalopathies., (© 2014 The Authors. **Human Mutation published by Wiley Periodicals, Inc.)

Respiratory chain deficiency can result from alterations in mitochondrial and/or cytosolic protein synthesis due to the dual genetic origin of mitochondrial oxidative phosphorylation. In the present paper we report a point mutation (D750G) in the bifunctional VARS (valyl-tRNA synthetase) of the fungus Neurospora crassa, associated with a temperature-sensitive phenotype. Analysis of the mutant strain revealed decreased steady-state levels of VARS and a clear reduction in the rate of mitochondrial protein synthesis. We observed a robust induction of the mitochondrial alternative oxidase with a concomitant decrease in the canonical respiratory pathway, namely in cytochrome b and aa3 content. Furthermore, the mutant strain accumulates the peripheral arm of complex I and depicts decreased levels of complexes III and IV, consistent with severe impairment of the mitochondrial respiratory chain. The phenotypic alterations of the mutant strain are observed at the permissive growth temperature and exacerbated upon increase of the temperature. Surprisingly, glucose-6-phosphate dehydrogenase activities were similar in the wild-type and mutant strains, whereas mitochondrial activities for succinate dehydrogenase and alternative NADH dehydrogenases were increased in the mutant strain, suggesting that the VARSD-G mutation does not affect overall cytosolic protein synthesis. Expression of the wild-type vars gene rescues all of the mutant phenotypes, indicating that the VARSD-G mutation is a loss-of-function mutation that results in a combined respiratory chain deficiency.

The present study analyzed the polymorphisms of DNA repair genes and their impact on the survival of 240 patients with early breast cancer. The genomic DNA was extracted from paraffin-embedded tumor-free tissue or blood, and thirteen polymorphisms in 12 DNA repair genes were determined using the Sequenom Mass array system. Among the target polymorphisms, VARS2 rs2074511 and POLE rs5744857 were found to correlate with survival after curative surgery in a log-rank test. No difference was found in the clinical and tumor characteristics according to the genotypes of these two coding variants, except for a higher incidence of positive ER in patients with the GG genotype of POLE rs5744857 (P=0.025). Meanwhile, a multivariate analysis showed that the GG genotype of VARS2 V552 V (rs2301717) was significantly associated with disease-free survival (DFS; HR=0.298; P=0.044) and marginally with distant DFS (DDFS; HR=0.266; P=0.077). In particular, the VARS2 rs2074511 polymorphism was only associated with survival in patients with triple negative (TN)-type breast cancer (P=0.018 for DFS and 0.042 for DDFS, respectively). In conclusion, VARS2 V552V may be considered as a prognostic factor for survival in patients with early breast cancer.

Previous studies showed that VAS1 of Saccharomyces cerevisiae encodes both cytosolic and mitochondrial forms of valyl-tRNA synthetase (ValRS) through alternative initiation of translation. We show herein that except for Schizosaccharomyces pombe, all yeast species studied contained a single ValRS gene encoding both forms, and all of the mature protein forms deduced from those genes possessed an N-terminal appended domain (Ad) that was absent from their bacterial relatives. In contrast, S. pombe contained two distinct nuclear ValRS genes, one encoding the mitochondrial form and the other its cytosolic counterpart. Although the cytosolic form closely resembles other yeast ValRS sequences (approximately 60% identity), the mitochondrial form exhibits significant divergence from others (approximately 35% identity). Both genes are active and essential for the survival of the yeast. Most conspicuously, the mitochondrial form lacks the characteristic Ad. A phylogenetic analysis further suggested that both forms of S. pombe ValRS are of mitochondrial origin, and the mitochondrial form is ancestral to the cytoplasmic form.

Since evidence indicates that the Last Universal Common Ancestor (LUCA) was phylogenetically closest to Methanopyrus kandleri among living organisms with elucidated genomes, this study has been directed to a search for the most primitive Methanopyrus lineage. For this purpose, the divergence of valyl-tRNA synthetase (ValRS) and isoleucyl-tRNA synthetase (IleRS) was employed as a measure of primitivity. Comparison of Methanopyrus kandleri and the Methanopyrus isolates GC34 and GC37 from the Pacific Ocean and KOL6, TAG1, TAG11, and SNP6 from the Atlantic Ocean established that the Pacific lineages are more primitive than the Atlantic lineages. Both the groups, however, are younger than environmental genomes from the Kairei Field of Central Indian Ridge in the Indian Ocean. These results showed that different Methanopyrus isolates differ significantly with respect to ValRS-IleRS divergence. On this basis, genomes giving rise to the ValRS and IleRS gene fragments from the Central Indian Ridge represent the most primitive Methanopyrus, phylogenetically the oldest living lineage closest to LUCA.

Previous studies showed that cytoplasmic and mitochondrial forms of yeast valyl-tRNA synthetase (ValRS) are specified by the VAS1 gene through alternative initiation of translation. Sequence comparison suggests that the yeast cytoplasmic (or mature mitochondrial) ValRS contains an N-terminal appendage that acts in cis as a nonspecific tRNA-binding domain (TRBD) and is absent from its bacterial relatives. We show here that Escherichia coli ValRS can substitute for the mitochondrial and cytoplasmic functions of VAS1 by fusion of a mitochondrial targeting signal and a TRBD, respectively. In addition, the bacterial ValRS gene can be converted into a dual functional yeast gene encoding both cytoplasmic and mitochondrial activities by fusion of a DNA sequence specifying both the mitochondrial targeting signal and TRBD. In vitro assays suggested that fusion of a nonspecific TRBD to the bacterial enzyme significantly enhanced its yeast tRNA-binding and aminoacylation activities. These results not only underscore the necessity of retaining a TRBD for functioning of a tRNA synthetase in yeast cytoplasm, but also provide insights into the evolution of tRNA synthetase genes.

The main role of aminoacyl-tRNA synthetases (aaRSs) is to transfer the cognate amino acids to the 3'-end of their tRNA by strictly discriminating from non-cognate amino acids. Some aaRSs accomplish this via proofreading and editing mechanisms, among which valyl-tRNA synthetase (ValRS) hydrolyses the non-cognate amino acid, threonine. In ValRS, existence of pre-transfer editing process is still unclear, although crystal structure of editing site with pre-transfer substrate analog (Thr-AMS) was released. In the case of isoleucyl-tRNA synthetase (IleRS), editing mechanism is well studied and mutational analyses revealed the existence of post- and pre-transfer editing mechanisms. Our aim is to investigate the possibility of pre-transfer editing process by performing molecular dynamics (MD) simulation studies. Simulations were carried out for ValRS with pre-transfer substrates (Thr-AMP/Val-AMP) and post-transfer substrates (Thr-A76/Val-A76) to understand their binding pattern. Two important point mutation studies were performed to observe their effect on editing process. This study also intends to compare and contrast the pre-transfer editing with post-transfer editing of ValRS. Interestingly, the MD simulation results revealed that non-cognate substrates (Thr-AMP/Thr-A76) bind more strongly than the cognate substrates (Val-AMP/Val-A76) in both pre- and post-transfer editing respectively. The editing site mutations (Lys270Ala and Asp279Ala) severely affected the binding ability of pre-transfer substrate (Thr-AMP) by different ways. Even though pre- and post-transfer substrates bind to the same site, specific differences were observed which has led us to believe the existence of the pre-transfer editing process in ValRS.

Previous studies showed that valyl-tRNA synthetase of Saccharomyces cerevisiae contains an N-terminal polypeptide extension of 97 residues, which is absent from its bacterial relatives, but is conserved in its mammalian homologues. We showed herein that this appended domain and its human counterpart are both nonspecific tRNA-binding domains (K(d) approximately 0.5 microm). Deletion of the appended domain from the yeast enzyme severely impaired its tRNA binding, aminoacylation, and complementation activities. This N-domain-deleted yeast valyl-tRNA synthetase mutant could be rescued by fusion of the equivalent domain from its human homologue. Moreover, fusion of the N-domain of the yeast enzyme or its human counterpart to Escherichia coli glutaminyl-tRNA synthetase enabled the otherwise "inactive" prokaryotic enzyme to function as a yeast enzyme in vivo. Different from the native yeast enzyme, which showed different affinities toward mixed tRNA populations, the fusion enzyme exhibited similar binding affinities for all yeast tRNAs. These results not only underscore the significance of nonspecific tRNA binding in aminoacylation, but also provide insights into the mechanism of the formation of aminoacyl-tRNAs.

Aminoacyl-tRNA synthetases prevent mistranslation, or genetic code ambiguity, through specialized editing reactions. Mutations that disrupt editing in bacteria adversely affect cell growth and viability, and recent work in the mouse supports the idea that translational errors caused by an editing defect lead to a neurological disease-like phenotype. To further investigate the connection of mistranslation to cell pathology, we introduced an inducible transgene expressing an editing-deficient valyl-tRNA synthetase into mammalian cells. Introducing mistranslation precipitated a disruption of cell morphology and membrane blebbing, accompanied by activation of caspase-3, consistent with an apoptotic response. Addition of a noncanonical amino acid that is misactivated, but not cleared, by the editing-defective enzyme exacerbated these effects. A special ambiguity-detecting sensor provided direct readout of mistranslation in vivo, supporting the possibility that decreased translational fidelity could be associated with disease.

The definition of the node of the last universal common ancestor (LUCA) is justified in a topology of the unrooted universal tree. This definition allows previous analyses based on paralogous proteins to be extended to orthologous ones. In particular, the use of a thermophily index (based on the amino acids' propensity to enter the [hyper] thermophile proteins more frequently) and its correlation with the optimal growth temperature of the various organisms allow inferences to be made on the habitat in which the LUCA lived. The reconstruction of ancestral sequences by means of the maximum likelihood method and their attribution to the set of mesophilic or hyperthermophilic sequences have led to the following conclusions: the LUCA was a hyperthermophile "organism," as were the ancestors of the Archaea and Bacteria domains, while the ancestor of the Eukarya domain was a mesophile. These conclusions are independent of the presence of hyperthermophile bacteria in the sample of sequences used in the analysis and are therefore independent of whether or not these are the first lines of divergence in the Bacteria domain, as observed in the topology of the universal tree of ribosomal RNA. These conclusions are thus more easily understood under the hypothesis that the origin of life took place at a high temperature.

Phylogenetic evidence is presented that primitively amitochondriate eukaryotes containing the nucleus, cytoskeleton, and endomembrane system may have never existed. Instead, the primary host for the mitochondrial progenitor may have been a chimeric prokaryote, created by fusion between an archaebacterium and a eubacterium, in which eubacterial energy metabolism (glycolysis and fermentation) was retained. A Rickettsia-like intracellular symbiont, suggested to be the last common ancestor of the family Rickettsiaceae and mitochondria, may have penetrated such a host (pro-eukaryote), surrounded by a single membrane, due to tightly membrane-associated phospholipase activity, as do present-day rickettsiae. The relatively rapid evolutionary conversion of the invader into an organelle may have occurred in a safe milieu via numerous, often dramatic, changes involving both partners, which resulted in successful coupling of the host glycolysis and the symbiont respiration. Establishment of a potent energy-generating organelle made it possible, through rapid dramatic changes, to develop genuine eukaryotic elements. Such sequential, or converging, global events could fill the gap between prokaryotes and eukaryotes known as major evolutionary discontinuity.

Previous studies showed that yeast VAS1 encodes both the cytoplasmic and mitochondrial forms of valyl-tRNA synthetase (ValRS), using alternative transcription and translation. The ValRS isoforms have identical polypeptide sequences, except for a 46-amino acid leader peptide that functions as a mitochondrial targeting signal. Although the two forms of the enzyme exhibit indistinguishable tRNA specificities in vitro, they cannot substitute for each other in vivo because of their different localizations. Here we show that the 46-residue leader sequence can be divided into two nonoverlapping peptides, each of which retains the ability to target the enzyme into mitochondria. The engineered proteins (with truncated leader sequences) are dual-targeted, rescuing both the cytoplasmic and mitochondrial defects of a vas1 knockout strain. Thus, in addition to alternative splicing and alternative translation initiation as mechanisms by which a single gene can encode cytoplasmic and mitochondrial activities, the inherent characteristics of a single polypeptide may enable it to be distributed simultaneously between two cellular compartments. This mechanism may explain how certain other single genes in Saccharomyces cerevisiae provide dual functions.

Valyl-tRNA synthetase (ValRS) from Escherichia coli undergoes covalent valylation by a donor valyl adenylate synthesized by the enzyme itself. ValRS could also be modified, although to a lesser extent, by the noncognate isosteric substrate L-threonine from a donor threonyl adenylate synthesized by the synthetase itself, or by the nonsubstrate methionine from methionyl adenylate produced by catalytic amounts of methionyl-tRNA synthetase. MALDI mass spectrometry analysis designated lysines 154, 162, 170, 533, 554, 593, 894, 930, and 940 of ValRS as the target residues for the attachment of valine. Following autothreonylation, lysines 162, 170, 178, 277, 291, 554, 580, 593, 861, 894, and 930 were found to be modified. Finally, L-Met-labeled residues were lysines 118, 162, 170, 178, 277, and 938. Alignment of the available ValRS amino acid sequences showed that lysines 277 and 554 are strictly conserved (with the exception concerning replacement of Lys-277 with a methionine or a tyrosine in archaebacteria), suggesting that these residues might be functionally significant. Indeed, lysine 554 of ValRS is the first lysine of the Lys-Met-Ser-Lys-Ser signature of the catalytic site of class I aminoacyl-tRNA synthetases. Lys-277 which is labeled by L-threonine or L-methionine, and not by L-valine, is located at or near the editing site, in the three-dimensional structure of ValRS. The role of lysine 277 was evaluated by site-directed mutagenesis. The Lys277Ala mutant (K277A) exhibited a posttransfer Thr-tRNA(Val) editing rate that was significantly lower than that observed for the wild-type enzyme. In addition, the K277A substitution altered amino acid discrimination in the editing site, resulting in hydrolysis of the correctly charged cognate Val-tRNA(Val). Finally, significant amounts of mischarged Thr-tRNA(Val) were produced by the K277A mutant, and not by wild-type ValRS. Altogether, our results designate Lys-277 as a likely candidate for nucleophilic attack of misacylated tRNA in the editing site of ValRS.

The rules of the genetic code are established in reactions that aminoacylate tRNAs with specific amino acids. Ambiguity in the code is prevented by editing activities whereby incorrect aminoacylations are cleared by specialized hydrolytic reactions of aminoacyl tRNA synthetases. Whereas editing reactions have long been known, their significance for cell viability is still poorly understood. Here we investigated in vitro and in vivo four different mutations in the center for editing that diminish the proofreading activity of valyl-tRNA synthetase (ValRS). The four mutant enzymes were shown to differ quantitatively in the severity of the defect in their ability to clear mischarged tRNA in vitro. Strikingly, in the presence of excess concentrations of alpha-aminobutyrate, one of the amino acids that is misactivated by ValRS, growth of bacterial strains bearing these mutant alleles is arrested. The concentration of misactivated amino acid required for growth arrest correlates inversely in a rank order with the degree of the editing defect seen in vitro. Thus, cell viability depends directly on the suppression of genetic code ambiguity by these specific editing reactions and is finely tuned to any perturbation of these reactions.

Valyl-tRNA synthetase (ValRS) has difficulty differentiating valine from structurally similar non-cognate amino acids, most prominently threonine. To minimize errors in aminoacylation and translation the enzyme catalyzes a proofreading (editing) reaction that is dependent on the presence of cognate tRNA(Val). Editing occurs at a site functionally distinct from the aminoacylation site of ValRS and previous results have shown that the 3'-terminus of tRNA(Val) is recognized differently at the two sites. Here, we extend these studies by comparing the contribution of aminoacylation identity determinants to productive recognition of tRNA(Val) at the aminoacylation and editing sites, and by probing tRNA(Val) for editing determinants that are distinct from those required for aminoacylation. Mutational analysis of Escherichia coli tRNA(Val) and identity switch experiments with non-cognate tRNAs reveal a direct relationship between the ability of a tRNA to be aminoacylated and its ability to stimulate the editing activity of ValRS. This suggests that at least a majority of editing by the enzyme entails prior charging of tRNA and that misacylated tRNA is a transient intermediate in the editing reaction.

Aminoacyl transfer RNA (tRNA) synthetases establish the rules of the genetic code by catalyzing the aminoacylation of tRNAs. For some synthetases, accuracy depends critically on an editing function at a site distinct from the aminoacylation site. Mutants of Escherichia coli that incorrectly charge tRNA(Val) with cysteine were selected after random mutagenesis of the whole chromosome. All mutations obtained were located in the editing site of valyl-tRNA synthetase. More than 20% of the valine in cellular proteins from such an editing mutant organism could be replaced with the noncanonical aminobutyrate, sterically similar to cysteine. Thus, the editing function may have played a central role in restricting the genetic code to 20 amino acids. Disabling this editing function offers a powerful approach for diversifying the chemical composition of proteins and for emulating evolutionary stages of ambiguous translation.

The evolutionary relationships of proteobacteria, which comprise the largest and phenotypically most diverse division among prokaryotes, are examined based on the analyses of available molecular sequence data. Sequence alignments of different proteins have led to the identification of numerous conserved inserts and deletions (referred to as signature sequences), which either are unique characteristics of various proteobacterial species or are shared by only members from certain subdivisions of proteobacteria. These signature sequences provide molecular means to define the proteobacterial phyla and their various subdivisions and to understand their evolutionary relationships to the other groups of eubacteria as well as the eukaryotes. Based on signature sequences that are present in different proteins it is now possible to infer that the various eubacterial phyla evolved from a common ancestor in the following order: low-G+C Gram-positive-->high-G+C Gram-positive-->Deinococcus-Thermus (green nonsulfur bacteria)-->cyanobacteria-->Spirochetes-->Chlamydia-Cytophaga-Aquifex -green sulfur bacteria-->Proteobacteria-1 (epsilon and delta)-->Proteobacteria-2 (alpha)-->Proteobacteria-3 (beta)-->Proteobacteria-4 (gamma). An unexpected but important aspect of the relationship deduced here is that the main eubacterial phyla are related to each other linearly rather than in a tree-like manner, suggesting that the major evolutionary changes within Bacteria have taken place in a directional manner. The identified signatures permit placement of prokaryotes into different groups/divisions and could be used for determinative purposes. These signatures generally support the origin of mitochondria from an alpha-proteobacterium and provide evidence that the nuclear cytosolic homologs of many genes are also derived from proteobacteria.

We characterized two Arabidopsis thaliana cDNAs coding for class I valyl-tRNA synthetase and class II threonyl-tRNA synthetase. The proteins display characteristics of cytosolic enzymes, yet possess an N-terminal extension relative to their prokaryotic homologs. The proximal part of the N-terminal extension is a mitochondrial-targeting signal. Through transient expression of GFP fusions in tobacco cells, we demonstrated that both genes encode the cytosolic and mitochondrial forms of the enzymes by alternative use of two in-frame initiation codons. A long, mitochondrial form of the enzyme is translated from a first initiation codon at reduced levels because of a poor sequence context and a shorter, cytosolic form is translated from a second in-phase AUG, which is in a better context for translation initiation. Primer extension experiments revealed several transcript ends mapping upstream of the first AUG and between the two AUGs. Distal to the mitochondrial transit peptide both valyl-tRNA synthetase and threonyl tRNA synthetase possess an NH2-appended domain compared with their prokaryotic counterparts. This domain's amphiphilic helix is conserved between yeast and A. thaliana valyl-tRNA synthetase, suggesting an important role in translation. Based on the high structural similarities between yeast and A. thaliana valyl-tRNA synthetase, we propose that the acquisition of bifunctionality of valyl-tRNA synthetase predates the divergence of these two organisms.

We have studied the interactions between Escherichia coli tRNAVal and valyl-tRNA synthetase (ValRS) by enzymatic footprinting with nuclease S1 and ribonuclease V1, and by analysis of the aminoacylation kinetics of mutant tRNAVal transcripts. Valyl-tRNA synthetase specifically protects the anticodon loop, the 3' side of the stacked T-stem/acceptor-stem helix, and the 5' side of the anticodon stem of tRNAVal against cleavage by double- and single-strand-specific nucleases. Increased nuclease susceptibility at the ends of the anticodon- and T-stems in the tRNAVal.ValRS complex is indicative of enzyme-induced conformational changes in the tRNA. The most important synthetase recognition determinants are the middle and 3' anticodon nucleotides (A35 and C36, respectively); G20, in the variable pocket, and G45, in the tRNA central core, are minor recognition elements. The discriminator base, position 73, and the anticodon stem also are recognized by ValRS. Replacing wild-type A73 with G73 reduces the aminoacylation efficiency more than 40-fold. However, the C73 and U73 mutants remain good substrates for ValRS, suggesting that guanosine at position 73 acts as a negative determinant. The amino acid acceptor arm of tRNAVal contains no other synthetase recognition nucleotides, but regular A-type RNA helix geometry in the acceptor stem is essential [Liu, M., et al. (1997) Nucleic Acids Res. 25, 4883-4890]. In the anticodon stem, converting the U29:A41 base pair to C29:G41 reduces the aminoacylation efficiency 50-fold. This is apparently due to the rigidity of the anticodon stem caused by the presence of five consecutive C:G base pairs, since the A29:U41 mutant is readily aminoacylated. Identity switch experiments provide additional evidence for a role of the anticodon stem in synthetase recognition. The valine recognition determinants, A35, C36, A73, G20, G45, and a regular A-RNA acceptor helix are insufficient to transform E. coli tRNAPhe into an effective valine acceptor. Replacing the anticodon stem of tRNAPhe with that of tRNAVal, however, converts the tRNA into a good substrate for ValRS. These experiments confirm G45 as a minor ValRS recognition element.

Nuclear-coded valyl-tRNA synthetase (ValRS) of eukaryotes is regarded of mitochondrial origin. Complete ValRS sequences obtained by us from two amitochondriate protists, the diplomonad, Giardia lamblia and the parabasalid, Trichomonas vaginalis were of the eukaryotic type, strongly suggesting an identical history of ValRS in all eukaryotes studied so far. The findings indicate that diplomonads are secondarily amitochondriate and give further evidence for such conclusion reached recently concerning parabasalids. Together with similar findings on other amitochondriate groups (microsporidia and entamoebids), this work provides critical support for the emerging notion that no representatives of the premitochondrial stage of eukaryotic phylogenesis exist among the species living today.

Phylogenetic analyses indicated that a series of paralogous gene pairs, found in two extensive regions on human chromosomal bands 6p21.3 and 9q33-34, were created by at least two independent duplications. The duplicated genes on chromosomal band 6p21.3 include the genes for type 11 collagen alpha2 subunit (COL11A2), NOTCH4 (mouse int-3 homologue), 70 kDa heat shock protein (HSPA1A, HSPA1B, and HSPA1L), valyl-tRNA synthetase 2 (VARS2), complement components (C2 and C4), pre-B cell leukemia transcription factor 2 (PBX2), retinoid X receptor beta (RXRB), NAT/RING3, and four other proteins. Their paralogous genes on chromosomal band 9q33-34 are genes for type 5 collagen alpha1 subunit (COL5A1), NOTCH1, 78 kDa glucose-regulated protein (HSPA5), valyl-tRNA synthetase 1 (VARS1), complement component V (C5), PBX3, retinoid X receptor alpha (RXRA), ORFX/RING3L, and others. Among these, the genes for collagen, complement components, NAT/RING3, PBX, and RXR appear to have been duplicated around the time of vertebrate emergence, supporting the idea that they were duplicated simultaneously at that time. Another group of genes that includes NOTCH and HSP appear to have diverged long before that time. A comparison of the physical maps of these two regions revealed that the genes which duplicated in the same period were arranged in almost the same order in the two regions, with the assumption of a few chromosomal rearrangements. We propose a possible model for the evolution of these regions, taking into account the molecular mechanisms of regional duplication, gene duplication, translocation, and inversion. We also propose that a comparative mapping of paralogous genes within the human genome would be useful for identifying new genes.