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.

Phenylketonuria is an inherited disease caused by mutations in the phenylalanine hydroxylase gene PAH. Different PAH pathogenic variants occur in different ethnic groups with various frequencies and the incidence of the disease itself varies from country to country. In the Caucasus region of Russia, some ethnoses are geographically and culturally isolated from each other. The tradition of monoethnic marriages may cause decreased genetic variability in those populations. In the Karachay-Cherkess Republic (Russia), the highest incidence of phenylketonuria in the world has been detected (1:850 newborns) in the region and 1:332 among the titular nation Karachays. Here, we showed that this phenomenon is due to the widespread prevalence of the p.Arg261* variant. Its allele frequency among Karachay patients with PKU was 68.4% and the carrier frequency in Karachays was 1:16 healthy individuals. PAH haplotype analysis showed a unique common origin. The founder haplotype and mutation “age” were estimated by analyzing the linkage disequilibrium between p.Arg261* and extragenic short tandem repeat loci. The p.Arg261* variant occurred in the Karachays population 10.2 ± 2.7 generations ago (275 ± 73 years) and its spread occurred in parallel with the growth of the population.

Phenylalanine hydroxylase (PAH) deficiency is responsible for most cases of phenylketonuria (PKU). Furthermore, numerous studies on BH4-sensitive PAH deficiency have been conducted. To date, BH4, a cofactor of PAH, has not been used to treat PKU in Russia.Genotype data of patients with PKU can be used to predict their sensitivity to BH4 therapy. A cohort of 2579 patients with PKU from Russia was analyzed for 25 common PAH gene mutations using custom allele-specific multiplex ligation-dependent probe amplification-based technology. A mutation detection rate of 84.1% chromosomes was accomplished. Both pathogenic alleles were identified in 73.1% of patients. The most frequent pathogenic variants were p.Arg408Trp (50.9%), p.Arg261Gln (5.3%), p.Pro281Leu (3.5%), IVS12+1G>A (3.1%), IVS10-11G>A (2.6%), and p.Arg158Leu (2.4%). The exact boundaries of a PAH exon 5 deletion were defined as EX5del4154ins268 (c.442-2913_509+1173del4154ins268). Severe phenotypes prevailed in the cohort, and classical PKU was observed in 71.8% cases. Due to the genotype-based prediction, 55.9% of the probands were non-responders to the BH4-treatment, and 20.2% were potential responders. Analysis of genotype data is useful to predict BH4 response in PKU patients. The high rate of non-responders among Russian patients was due to the high allele frequency of severe PAH mutations.

A timely detection of patients with tetrahydrobiopterin (BH4) -deficient types of hyperphenylalaninemia (HPABH4) is important for assignment of correct therapy, allowing to avoid complications. Often HPABH4 patients receive the same therapy as phenylalanine hydroxylase (PAH) -deficiency (phenylketonuria) patients—dietary treatment—and do not receive substitutive BH4 therapy until the diagnosis is confirmed by molecular genetic means. In this study, we present a cohort of 30 Russian patients with HPABH4 with detected variants in genes causing different types of HPA. Family diagnostics and biochemical urinary pterin spectrum analyses were carried out. HPABH4A is shown to be the prevalent type, 83.3% of all HPABH4 cases. The mutation spectrum for the PTS gene was defined, the most common variants in Russia were p.Thr106Met—32%, p.Asn72Lys—20%, p.Arg9His—8%, p.Ser32Gly—6%. We also detected 7 novel PTS variants and 3 novel QDPR variants. HPABH4 prevalence was estimated to be 0.5–0.9% of all HPA cases in Russia, which is significantly lower than in European countries on average, China, and Saudi Arabia. The results of this research show the necessity of introducing differential diagnostics for HPABH4 into neonatal screening practice.

Prevalence and allelic heterogeneity of hereditary diseases (HDs) could vary significantly in different human populations. Current knowledge of HDs distribution in populations is generally limited to either European data or analyses of isolated populations which were performed several decades ago. Thus, an acknowledgement of the HDs prevalence in different modern open populations is important. The study presents the results of a genetic epidemiological study of hereditary diseases (HDs) in the population of the Karachay-Cherkess Republic (KChR). Clinical screening of a population of 410,367 people for the identification of HDs was conducted. The population surveyed is represented by five major ethnic groups—Karachays, Russians, Circassians, Abazins, Nogais. The study of the populations was carried out in accordance with the proprietary protocol of genetic epidemiological examination designed to identify >3500 HDs easily diagnosed during clinical examination by qualified specialists specializing in the HDs. The protocol consists of the population genetic and medical genetic sections and is intended for comprehensive population analysis based on the data on different genetic systems, including the genes of HDs, DNA polymorphisms, demographic data collected during hospital-based survey. 8950 families (with 10,125 patients) with presumably the HDs were initially identified as a result of the survey and data collection through various sources of registration (from 1156 medical workers from 163 medical institutions). A diagnosis of hereditary pathology was established in 1849 patients (from 1295 families). Two hundred and thirty nosological forms were revealed (in 1857 patients from 1295 families). The total prevalence of HDs was 1:221. Differences between populations and ethnic groups were identified: 1:350 in Russians, 1:195 in Karachays, 1:199 in Circassians, 1:218 in Abazins, 1:135 in Nogais. Frequent diseases were determined, the presence of marked genetic heterogeneity was identified during the confirmatory DNA diagnosis. To explain the reasons for the differentiation of populations by load of HD, a correlation analysis was carried out between the FST (random inbreeding) in populations and HDs load values. This analysis showed genetic drift is probably one of the leading factors determining the differentiation of KChR populations by HDs load. For the first time, the size of the load and spectrum of HDs in the populations of the KChR are determined. We have demonstrated genetic drift to be one of the main factors of the population dynamics in studied population. A significant genetic heterogeneity of HDs, both allelic and locus, was revealed in KChR.