How Does the Ribosome Fold the Proteome?

Folding of polypeptides begins during their synthesis on ribosomes. This process has evolved as a means for the cell to maintain proteostasis, by mitigating the risk of protein misfolding and aggregation. The capacity to now depict this cellular feat at increasingly higher resolution is providing insight into the mechanistic determinants that promote successful folding. Emerging from these studies is the intimate interplay between protein translation and folding, and within this the ribosome particle is the key player. Its unique structural properties provide a specialized scaffold against which nascent polypeptides can begin to form structure in a highly coordinated, co-translational manner. Here, we examine how, as a macromolecular machine, the ribosome modulates the intrinsic dynamic properties of emerging nascent polypeptide chains and guides them toward their biologically active structures.

The ribosome lowers the entropic penalty of protein folding.

Most proteins fold during biosynthesis on the ribosome1, and co-translational folding energetics, pathways and outcomes of many proteins have been found to differ considerably from those in refolding studies2-10. The origin of this folding modulation by the ribosome has remained unknown. Here we have determined atomistic structures of the unfolded state of a model protein on and off the ribosome, which reveal that the ribosome structurally expands the unfolded nascent chain and increases its solvation, resulting in its entropic destabilization relative to the peptide chain in isolation. Quantitative 19F NMR experiments confirm that this destabilization reduces the entropic penalty of folding by up to 30 kcal mol-1 and promotes formation of partially folded intermediates on the ribosome, an observation that extends to other protein domains and is obligate for some proteins to acquire their active conformation. The thermodynamic effects also contribute to the ribosome protecting the nascent chain from mutation-induced unfolding, which suggests a crucial role of the ribosome in supporting protein evolution. By correlating nascent chain structure and dynamics to their folding energetics and post-translational outcomes, our findings establish the physical basis of the distinct thermodynamics of co-translational protein folding.

Multi-omics identifies large mitoribosomal subunit instability caused by pathogenic MRPL39 variants as a cause of pediatric onset mitochondrial disease.

MRPL39 encodes one of 52 proteins comprising the large subunit of the mitochondrial ribosome (mitoribosome). In conjunction with 30 proteins in the small subunit, the mitoribosome synthesizes the 13 subunits of the mitochondrial oxidative phosphorylation (OXPHOS) system encoded by mitochondrial Deoxyribonucleic acid (DNA). We used multi-omics and gene matching to identify three unrelated individuals with biallelic variants in MRPL39 presenting with multisystem diseases with severity ranging from lethal, infantile-onset (Leigh syndrome spectrum) to milder with survival into adulthood. Clinical exome sequencing of known disease genes failed to diagnose these patients; however quantitative proteomics identified a specific decrease in the abundance of large but not small mitoribosomal subunits in fibroblasts from the two patients with severe phenotype. Re-analysis of exome sequencing led to the identification of candidate single heterozygous variants in mitoribosomal genes MRPL39 (both patients) and MRPL15. Genome sequencing identified a shared deep intronic MRPL39 variant predicted to generate a cryptic exon, with transcriptomics and targeted studies providing further functional evidence for causation. The patient with the milder disease was homozygous for a missense variant identified through trio exome sequencing. Our study highlights the utility of quantitative proteomics in detecting protein signatures and in characterizing gene-disease associations in exome-unsolved patients. We describe Relative Complex Abundance analysis of proteomics data, a sensitive method that can identify defects in OXPHOS disorders to a similar or greater sensitivity to the traditional enzymology. Relative Complex Abundance has potential utility for functional validation or prioritization in many hundreds of inherited rare diseases where protein complex assembly is disrupted.

Shariant platform: Enabling evidence sharing across Australian clinical genetic-testing laboratories to support variant interpretation.

Sharing genomic variant interpretations across laboratories promotes consistency in variant assertions. A landscape analysis of Australian clinical genetic-testing laboratories in 2017 identified that, despite the national-accreditation-body recommendations encouraging laboratories to submit genotypic data to clinical databases, fewer than 300 variants had been shared to the ClinVar public database. Consultations with Australian laboratories identified resource constraints limiting routine application of manual processes, consent issues, and differences in interpretation systems as barriers to sharing. This information was used to define key needs and solutions required to enable national sharing of variant interpretations. The Shariant platform, using both the GRCh37 and GRCh38 genome builds, was developed to enable ongoing sharing of variant interpretations and associated evidence between Australian clinical genetic-testing laboratories. Where possible, two-way automated sharing was implemented so that disruption to laboratory workflows would be minimized. Terms of use were developed through consultation and currently restrict access to Australian clinical genetic-testing laboratories. Shariant was designed to store and compare structured evidence, to promote and record resolution of inter-laboratory classification discrepancies, and to streamline the submission of variant assertions to ClinVar. As of December 2021, more than 14,000 largely prospectively curated variant records from 11 participating laboratories have been shared. Discrepant classifications have been identified for 11% (28/260) of variants submitted by more than one laboratory. We have demonstrated that co-design with clinical laboratories is vital to developing and implementing a national variant-interpretation sharing effort. This approach has improved inter-laboratory concordance and enabled opportunities to standardize interpretation practices.

Novel CDKL5 targets identified in human iPSC-derived neurons.

CDKL5 Deficiency Disorder (CDD) is a debilitating epileptic encephalopathy disorder affecting young children with no effective treatments. CDD is caused by pathogenic variants in Cyclin-Dependent Kinase-Like 5 (CDKL5), a protein kinase that regulates key phosphorylation events in neurons. For therapeutic intervention, it is essential to understand molecular pathways and phosphorylation targets of CDKL5. Using an unbiased phosphoproteomic approach we identified novel targets of CDKL5, including GTF2I, PPP1R35, GATAD2A and ZNF219 in human iPSC-derived neuronal cells. The phosphoserine residue in the target proteins lies in the CDKL5 consensus motif. We validated direct phosphorylation of GTF2I and PPP1R35 by CDKL5 using complementary approaches. GTF2I controls axon guidance, cell cycle and neurodevelopment by regulating expression of neuronal genes. PPP1R35 is critical for centriole elongation and cilia morphology, processes that are impaired in CDD. PPP1R35 interacts with CEP131, a known CDKL5 phospho-target. GATAD2A and ZNF219 belong to the Nucleosome Remodelling Deacetylase (NuRD) complex, which regulates neuronal activity-dependent genes and synaptic connectivity. In-depth knowledge of molecular pathways regulated by CDKL5 will allow a better understanding of druggable disease pathways to fast-track therapeutic development.

Compound heterozygous variants in SHQ1 are associated with a spectrum of neurological features, including early-onset dystonia.

SHQ1 is essential for biogenesis of H/ACA ribonucleoproteins, a class of molecules important for processing ribosomal RNAs, modifying spliceosomal small nuclear RNAs and stabilizing telomerase. Components of the H/ACA ribonucleoprotein complex have been linked to neurological developmental defects. Here, we report two sibling pairs from unrelated families with compound heterozygous variants in SHQ1. Exome sequencing was used to detect disease causing variants, which were submitted to 'matching' platforms linked to MatchMaker Exchange. Phenotype comparisons supported these matches. The affected individuals present with early-onset dystonia, with individuals from one family displaying additional neurological phenotypes, including neurodegeneration. As a result of cerebrospinal fluid studies suggesting possible abnormal dopamine metabolism, a trial of levodopa replacement therapy was started but no clear response was noted. We show that fibroblasts from affected individuals have dramatic loss of SHQ1 protein. Variants from both families were expressed in Saccharomyces cerevisiae, resulting in a strong reduction in H/ACA snoRNA production and remarkable defects in rRNA processing and ribosome formation. Our study identifies SHQ1 as associated with neurological disease, including early-onset dystonia, and begins to delineate the molecular etiology of this novel condition.

Pathogenic variants in nucleoporin TPR (translocated promoter region, nuclear basket protein) cause severe intellectual disability in humans.

The nuclear pore complex (NPC) is a multi-protein complex that regulates the trafficking of macromolecules between the nucleus and cytoplasm. Genetic variants in components of the NPC have been shown to cause a range of neurological disorders, including intellectual disability and microcephaly. Translocated promoter region, nuclear basket protein (TPR) is a critical scaffolding element of the nuclear facing interior of the NPC. Here, we present two siblings with biallelic variants in TPR who present with a phenotype of microcephaly, ataxia and severe intellectual disability. The variants result in a premature truncation variant, and a splice variant leading to a 12-amino acid deletion respectively. Functional analyses in patient fibroblasts demonstrate significantly reduced TPR levels, and decreased TPR-containing NPC density. A compensatory increase in total NPC levels was observed, and decreased global RNA intensity in the nucleus. The discovery of variants that partly disable TPR function provide valuable insight into this essential protein in human disease, and our findings suggest that TPR variants are the cause of the siblings' neurological disorder.

A DNA repair disorder caused by de novo monoallelic DDB1 variants is associated with a neurodevelopmental syndrome.

The DNA damage-binding protein 1 (DDB1) is part of the CUL4-DDB1 ubiquitin E3 ligase complex (CRL4), which is essential for DNA repair, chromatin remodeling, DNA replication, and signal transduction. Loss-of-function variants in genes encoding the complex components CUL4 and PHIP have been reported to cause syndromic intellectual disability with hypotonia and obesity, but no phenotype has been reported in association with DDB1 variants. Here, we report eight unrelated individuals, identified through Matchmaker Exchange, with de novo monoallelic variants in DDB1, including one recurrent variant in four individuals. The affected individuals have a consistent phenotype of hypotonia, mild to moderate intellectual disability, and similar facies, including horizontal or slightly bowed eyebrows, deep-set eyes, full cheeks, a short nose, and large, fleshy and forward-facing earlobes, demonstrated in the composite face generated from the cohort. Digital anomalies, including brachydactyly and syndactyly, were common. Three older individuals have obesity. We show that cells derived from affected individuals have altered DDB1 function resulting in abnormal DNA damage signatures and histone methylation following UV-induced DNA damage. Overall, our study adds to the growing family of neurodevelopmental phenotypes mediated by disruption of the CRL4 ubiquitin ligase pathway and begins to delineate the phenotypic and molecular effects of DDB1 misregulation.

Clustered mutations in the GRIK2 kainate receptor subunit gene underlie diverse neurodevelopmental disorders.

Kainate receptors (KARs) are glutamate-gated cation channels with diverse roles in the central nervous system. Bi-allelic loss of function of the KAR-encoding gene GRIK2 causes a nonsyndromic neurodevelopmental disorder (NDD) with intellectual disability and developmental delay as core features. The extent to which mono-allelic variants in GRIK2 also underlie NDDs is less understood because only a single individual has been reported previously. Here, we describe an additional eleven individuals with heterozygous de novo variants in GRIK2 causative for neurodevelopmental deficits that include intellectual disability. Five children harbored recurrent de novo variants (three encoding p.Thr660Lys and two p.Thr660Arg), and four children and one adult were homozygous for a previously reported variant (c.1969G>A [p.Ala657Thr]). Individuals with shared variants had some overlapping behavioral and neurological dysfunction, suggesting that the GRIK2 variants are likely pathogenic. Analogous mutations introduced into recombinant GluK2 KAR subunits at sites within the M3 transmembrane domain (encoding p.Ala657Thr, p.Thr660Lys, and p.Thr660Arg) and the M3-S2 linker domain (encoding p.Ile668Thr) had complex effects on functional properties and membrane localization of homomeric and heteromeric KARs. Both p.Thr660Lys and p.Thr660Arg mutant KARs exhibited markedly slowed gating kinetics, similar to p.Ala657Thr-containing receptors. Moreover, we observed emerging genotype-phenotype correlations, including the presence of severe epilepsy in individuals with the p.Thr660Lys variant and hypomyelination in individuals with either the p.Thr660Lys or p.Thr660Arg variant. Collectively, these results demonstrate that human GRIK2 variants predicted to alter channel function are causative for early childhood development disorders and further emphasize the importance of clarifying the role of KARs in early nervous system development.

Gene selection for genomic newborn screening: Moving toward consensus?

PURPOSE: Gene selection for genomic newborn screening (gNBS) underpins the validity, acceptability, and ethical application of this technology. Existing gNBS gene lists are highly variable despite being based on shared principles of gene-disease validity, treatability, and age of onset. This study aimed to curate a gNBS gene list that builds upon existing efforts and provide a core consensus list of gene-disease pairs assessed by multiple expert groups worldwide. METHODS: Our multidisciplinary expert team curated a gene list using an open platform and multiple existing curated resources. We included severe treatable disorders with age of disease onset <5 years with established gene-disease associations and reliable variant detection. We compared the final list with published lists from 5 other gNBS projects to determine consensus genes and to identify areas of discrepancy. RESULTS: We reviewed 1279 genes and 604 met our inclusion criteria. Metabolic conditions comprised the largest group (25%), followed by immunodeficiencies (21%) and endocrine disorders (15%). We identified 55 consensus genes included by all 6 gNBS research projects. Common reasons for discrepancy included variable definitions of treatability and strength of gene-disease association. CONCLUSION: We have identified a consensus gene list for gNBS that can be used as a basis for systematic harmonization efforts internationally.

A national education program for rapid genomics in pediatric acute care: building workforce confidence, competence and capability.

PURPOSE: To develop and evaluate a scalable national program to build confidence, competence and capability in the use of rapid genomic testing in the acute pediatric setting. METHODS: We used theory-informed approaches to design a modular, adaptive program of blended learning aimed at diverse professional groups involved in acute pediatric care. The program comprised four online learning modules and an online workshop and was centered on case-based learning. We evaluated the program using the Kirkpatrick four-level model of training evaluation and report our findings using the RISE2 guidelines for genomics education and evaluation. RESULTS: Two hundred and two participants engaged with at least one component of the program. Participants self-reported increased confidence in using rGT, (p<0.001) and quiz responses objectively demonstrated increased competence (e.g., correct responses to a question on pre-test counseling increased from 30% to 64%; p<0.001). Additionally, their capability in applying genomic principles to simulated clinical cases increased (p<0.001), as did their desire to take on more responsibility for performing rGT. The clinical interpretation of more complex test results (such as negative results or variants of uncertain significance) appeared to be more challenging, indicating a need for targeted education in this area. CONCLUSION: The program format was effective in delivering multidisciplinary and wide-scale genomics education in the acute care context. The modular approach we have developed now lends itself to application in other medical specialties or areas of healthcare.

Further delineation of short-chain enoyl-CoA hydratase deficiency in the Pacific population.

Short-chain enoyl-coA hydratase (SCEH) deficiency due to biallelic pathogenic ECHS1 variants was first reported in 2014 in association with Leigh syndrome (LS) and increased S-(2-carboxypropyl)cysteine excretion. It is potentially treatable with a valine-restricted, high-energy diet and emergency regimen. Recently, Simon et al. described four Samoan children harbouring a hypomorphic allele (c.489G > A, p.Pro163=) associated with reduced levels of normally-spliced mRNA. This synonymous variant, missed on standard genomic testing, is prevalent in the Samoan population (allele frequency 0.17). Patients with LS and one ECHS1 variant were identified in NZ and Australian genomic and clinical databases. ECHS1 sequence data were interrogated for the c.489G > A variant and clinical data were reviewed. Thirteen patients from 10 families were identified; all had Pacific ancestry including Samoan, Maori, Cook Island Maori, and Tokelauan. All developed bilateral globus pallidi lesions, excluding one pre-symptomatic infant. Symptom onset was in early childhood, and was triggered by illness or starvation in 9/13. Four of 13 had exercise-induced dyskinesia, 9/13 optic atrophy and 6/13 nystagmus. Urine S-(2-carboxypropyl)cysteine-carnitine and other SCEH-related metabolites were normal or mildly increased. Functional studies demonstrated skipping of exon four and markedly reduced ECHS1 protein. These data provide further support for the pathogenicity of this ECHS1 variant which is also prevalent in Maori, Cook Island Maori, and Tongan populations (allele frequency 0.14-0.24). It highlights the need to search for a second variant in apparent heterozygotes with an appropriate phenotype, and has implications for genetic counselling in family members who are heterozygous for the more severe ECHS1 alleles. SYNOPSIS: Short-chain enoyl-CoA hydratase deficiency is a frequent cause of Leigh-like disease in Maori and wider-Pacific populations, due to the high carrier frequency of a hypomorphic ECHS1 variant c.489G > A, p.[Pro163=, Phe139Valfs*65] that may be overlooked by standard genomic testing.

De novo enhancer deletion of LMX1B produces a mild nail-patella clinical phenotype.

Critical genes involved in embryonic development are often transcription factors, regulating many downstream genes. LMX1B is a homeobox gene that is involved in formation of the limbs, eyes and kidneys, heterozygous loss-of-function sequence variants and deletions cause Nail-Patella syndrome. Most of the reported variants are localised within the gene's coding sequence, however, approximately 5%-10% of affected individuals do not have a pathogenic variant identified within this region. In this study, we present a family with four affected individuals across two generations with a deletion spanning a conserved upstream LMX1B-binding sequence. This deletion is de novo in the mother of three affected children. Furthermore, in this family, the manifestations appear limited to the nails and limbs, and therefore may reflect an attenuated phenotype of the classic Nail-Patella phenotype that includes ophthalmological and renal manifestations.

The collective burden of childhood dementia: a scoping review.

Childhood dementia is a devastating and under-recognized group of disorders with a high level of unmet need. Typically monogenic in origin, this collective of individual neurodegenerative conditions are defined by a progressive impairment of neurocognitive function, presenting in childhood and adolescence. This scoping review aims to clarify definitions and conceptual boundaries of childhood dementia and quantify the collective disease burden. A literature review identified conditions that met the case definition. An expert clinical working group reviewed and ratified inclusion. Epidemiological data were extracted from published literature and collective burden modelled. One hundred and seventy genetic childhood dementia disorders were identified. Of these, 25 were analysed separately as treatable conditions. Collectively, currently untreatable childhood dementia was estimated to have an incidence of 34.5 per 100 000 (1 in 2900 births), median life expectancy of 9 years and prevalence of 5.3 per 100 000 persons. The estimated number of premature deaths per year is similar to childhood cancer (0-14 years) and approximately 70% of those deaths will be prior to adulthood. An additional 49.8 per 100 000 births are attributable to treatable conditions that would cause childhood dementia if not diagnosed early and stringently treated. A relational database of the childhood dementia disorders has been created and will be continually updated as new disorders are identified (https://knowledgebase.childhooddementia.org/). We present the first comprehensive overview of monogenic childhood dementia conditions and their collective epidemiology. Unifying these conditions, with consistent language and definitions, reinforces motivation to advance therapeutic development and health service supports for this significantly disadvantaged group of children and their families.

Common sequence motifs of nascent chains engage the ribosome surface and trigger factor.

In the cell, the conformations of nascent polypeptide chains during translation are modulated by both the ribosome and its associated molecular chaperone, trigger factor. The specific interactions that underlie these modulations, however, are still not known in detail. Here, we combine protein engineering, in-cell and in vitro NMR spectroscopy, and molecular dynamics simulations to explore how proteins interact with the ribosome during their biosynthesis before folding occurs. Our observations of α-synuclein nascent chains in living Escherichia coli cells reveal that ribosome surface interactions dictate the dynamics of emerging disordered polypeptides in the crowded cytosol. We show that specific basic and aromatic motifs drive such interactions and directly compete with trigger factor binding while biasing the direction of the nascent chain during its exit out of the tunnel. These results reveal a structural basis for the functional role of the ribosome as a scaffold with holdase characteristics and explain how handover of the nascent chain to specific auxiliary proteins occurs among a host of other factors in the cytosol.

Nature and Regulation of Protein Folding on the Ribosome.

Co-translational protein folding is an essential process by which cells ensure the safe and efficient production and assembly of new proteins in their functional native states following biosynthesis on the ribosome. In this review, we describe recent progress in probing the changes during protein synthesis of the free energy landscapes that underlie co-translational folding and discuss the critical coupling between these landscapes and the rate of translation that ultimately determines the success or otherwise of the folding process. Recent developments have revealed a variety of mechanisms by which both folding and translation can be modulated or regulated, and we discuss how these effects are utilised by the cell to optimise the outcome of protein biosynthesis.

Bi-allelic LoF NRROS Variants Impairing Active TGF-ß1 Delivery Cause a Severe Infantile-Onset Neurodegenerative Condition with Intracranial Calcification.

Negative regulator of reactive oxygen species (NRROS) is a leucine-rich repeat-containing protein that uniquely associates with latent transforming growth factor beta-1 (TGF- ß1) and anchors it on the cell surface; this anchoring is required for activation of TGF-ß1 in macrophages and microglia. We report six individuals from four families with bi-allelic variants in NRROS. All affected individuals had neurodegenerative disease with refractory epilepsy, developmental regression, and reduced white matter volume with delayed myelination. The clinical course in affected individuals began with normal development or mild developmental delay, and the onset of seizures occurred within the first year of life, followed by developmental regression. Intracranial calcification was detected in three individuals. The phenotypic features in affected individuals are consistent with those observed in the Nrros knockout mouse, and they overlap with those seen in the human condition associated with TGF-ß1 deficiency. The disease-causing NRROS variants involve two significant functional NRROS domains. These variants result in aberrant NRROS proteins with impaired ability to anchor latent TGF-ß1 on the cell surface. Using confocal microscopy in HEK293T cells, we demonstrate that wild-type and mutant NRROS proteins co-localize with latent TGF-ß1 intracellularly. However, using flow cytometry, we show that our mutant NRROS proteins fail to anchor latent TGF-ß1 at the cell surface in comparison to wild-type NRROS. Moreover, wild-type NRROS rescues the defect of our disease-associated mutants in presenting latent TGF-ß1 to the cell surface. Taken together, our findings suggest that loss of NRROS function causes a severe childhood-onset neurodegenerative condition with features suggestive of a disordered response to inflammation.

Bi-allelic ADARB1 Variants Associated with Microcephaly, Intellectual Disability, and Seizures.

The RNA editing enzyme ADAR2 is essential for the recoding of brain transcripts. Impaired ADAR2 editing leads to early-onset epilepsy and premature death in a mouse model. Here, we report bi-allelic variants in ADARB1, the gene encoding ADAR2, in four unrelated individuals with microcephaly, intellectual disability, and epilepsy. In one individual, a homozygous variant in one of the double-stranded RNA-binding domains (dsRBDs) was identified. In the others, variants were situated in or around the deaminase domain. To evaluate the effects of these variants on ADAR2 enzymatic activity, we performed in vitro assays with recombinant proteins in HEK293T cells and ex vivo assays with fibroblasts derived from one of the individuals. We demonstrate that these ADAR2 variants lead to reduced editing activity on a known ADAR2 substrate. We also demonstrate that one variant leads to changes in splicing of ADARB1 transcript isoforms. These findings reinforce the importance of RNA editing in brain development and introduce ADARB1 as a genetic etiology in individuals with intellectual disability, microcephaly, and epilepsy.

Bi-allelic Variations of SMO in Humans Cause a Broad Spectrum of Developmental Anomalies Due to Abnormal Hedgehog Signaling.

The evolutionarily conserved hedgehog (Hh) pathway is essential for organogenesis and plays critical roles in postnatal tissue maintenance and renewal. A unique feature of the vertebrate Hh pathway is that signal transduction requires the primary cilium (PC) where major pathway components are dynamically enriched. These factors include smoothened (SMO) and patched, which constitute the core reception system for sonic hedgehog (SHH) as well as GLI transcription factors, the key mediators of the pathway. Here, we report bi-allelic loss-of-function variations in SMO in seven individuals from five independent families; these variations cause a wide phenotypic spectrum of developmental anomalies affecting the brain (hypothalamic hamartoma and microcephaly), heart (atrioventricular septal defect), skeleton (postaxial polydactyly, narrow chest, and shortening of long bones), and enteric nervous system (aganglionosis). Cells derived from affected individuals showed normal ciliogenesis but severely altered Hh-signal transduction as a result of either altered PC trafficking or abnormal activation of the pathway downstream of SMO. In addition, Hh-independent GLI2 accumulation at the PC tip in cells from the affected individuals suggests a potential function of SMO in regulating basal ciliary trafficking of GLI2 when the pathway is off. Thus, loss of SMO function results in abnormal PC dynamics of key components of the Hh signaling pathway and leads to a large continuum of malformations in humans.

The Genetic Landscape and Epidemiology of Phenylketonuria.

Phenylketonuria (PKU), caused by variants in the phenylalanine hydroxylase (PAH) gene, is the most common autosomal-recessive Mendelian phenotype of amino acid metabolism. We estimated that globally 0.45 million individuals have PKU, with global prevalence 1:23,930 live births (range 1:4,500 [Italy]-1:125,000 [Japan]). Comparing genotypes and metabolic phenotypes from 16,092 affected subjects revealed differences in disease severity in 51 countries from 17 world regions, with the global phenotype distribution of 62% classic PKU, 22% mild PKU, and 16% mild hyperphenylalaninemia. A gradient in genotype and phenotype distribution exists across Europe, from classic PKU in the east to mild PKU in the southwest and mild hyperphenylalaninemia in the south. The c.1241A>G (p.Tyr414Cys)-associated genotype can be traced from Northern to Western Europe, from Sweden via Norway, to Denmark, to the Netherlands. The frequency of classic PKU increases from Europe (56%) via Middle East (71%) to Australia (80%). Of 758 PAH variants, c.1222C>T (p.Arg408Trp) (22.2%), c.1066-11G>A (IVS10-11G>A) (6.4%), and c.782G>A (p.Arg261Gln) (5.5%) were most common and responsible for two prevalent genotypes: p.[Arg408Trp];[Arg408Trp] (11.4%) and c.[1066-11G>A];[1066-11G>A] (2.6%). Most genotypes (73%) were compound heterozygous, 27% were homozygous, and 55% of 3,659 different genotypes occurred in only a single individual. PAH variants were scored using an allelic phenotype value and correlated with pre-treatment blood phenylalanine concentrations (n = 6,115) and tetrahydrobiopterin loading test results (n = 4,381), enabling prediction of both a genotype-based phenotype (88%) and tetrahydrobiopterin responsiveness (83%). This study shows that large genotype databases enable accurate phenotype prediction, allowing appropriate targeting of therapies to optimize clinical outcome.