Genome-wide sequence analyses of ethnic populations across Russia

The Spectrum of Disease-Associated Alleles in Countries with a Predominantly Slavic Population

There are more than 260 million people of Slavic descent worldwide, who reside mainly in Eastern Europe but also represent a noticeable share of the population in the USA and Canada. Slavic populations, particularly Eastern Slavs and some Western Slavs, demonstrate a surprisingly high degree of genetic homogeneity, and, consequently, remarkable contribution of recurrent alleles associated with hereditary diseases. Along with pan-European pathogenic variants with clearly elevated occurrence in Slavic people (e.g., ATP7B c.3207C>A and PAH c.1222C>T), there are at least 52 pan-Slavic germ-line mutations (e.g., NBN c.657_661del and BRCA1 c.5266dupC) as well as several disease-predisposing alleles characteristic of the particular Slavic communities (e.g., Polish SDHD c.33C>A and Russian ARSB c.1562G>A variants). From a clinical standpoint, Slavs have some features of a huge founder population, thus providing a unique opportunity for efficient genetic studies.

Complex trait susceptibilities and population diversity in a sample of 4,145 Russians

The population of Russia consists of more than 150 local ethnicities. The ethnic diversity and geographic origins, which extend from eastern Europe to Asia, make the population uniquely positioned to investigate the shared properties of inherited disease risks between European and Asian ancestries. We present the analysis of genetic and phenotypic data from a cohort of 4,145 individuals collected in three metro areas in western Russia. We show the presence of multiple admixed genetic ancestry clusters spanning from primarily European to Asian and high identity-by-descent sharing with the Finnish population. As a result, there was notable enrichment of Finnish-specific variants in Russia. We illustrate the utility of Russian-descent cohorts for discovery of novel population-specific genetic associations, as well as replication of previously identified associations that were thought to be population-specific in other cohorts. Finally, we provide access to a database of allele frequencies and GWAS results for 464 phenotypes.

Next-generation data filtering in the genomics era

Genomic data are ubiquitous across disciplines, from agriculture to biodiversity, ecology, evolution and human health. However, these datasets often contain noise or errors and are missing information that can affect the accuracy and reliability of subsequent computational analyses and conclusions. A key step in genomic data analysis is filtering - removing sequencing bases, reads, genetic variants and/or individuals from a dataset - to improve data quality for downstream analyses. Researchers are confronted with a multitude of choices when filtering genomic data; they must choose which filters to apply and select appropriate thresholds. To help usher in the next generation of genomic data filtering, we review and suggest best practices to improve the implementation, reproducibility and reporting standards for filter types and thresholds commonly applied to genomic datasets. We focus mainly on filters for minor allele frequency, missing data per individual or per locus, linkage disequilibrium and Hardy-Weinberg deviations. Using simulated and empirical datasets, we illustrate the large effects of different filtering thresholds on common population genetics statistics, such as Tajima's D value, population differentiation (FST), nucleotide diversity (π) and effective population size (Ne).

Association of GNB3, ACE polymorphisms with POAG and NTG

PURPOSE

Primary open-angle glaucoma (POAG) represents the most prevalent form of glaucoma and stands as a foremost contributor to irreversible vision impairment on a global scale. Despite notable strides made in comprehending the genetic underpinnings of POAG, investigations within the context of Russia remain constrained.

METHODS

The study cohort comprised a total of 235 individuals, with 135 of them exhibiting various forms of glaucoma encompassing both POAG and (NTG, while the remaining 100 individuals served as control subjects. Each participant underwent a comprehensive ocular examination to ascertain their ocular health status. Genotyping of the relevant single nucleotide polymorphisms (SNPs) was carried out using the Taq Man genotyping assay. Specifically, the two SNPs under scrutiny were GNB3 rs5443 gene and ACE rs4646994. Statistical analysis was performed to evaluate the association of these SNPs with glaucoma risk.

RESULTS

The presence of the T allele of rs5443 was found to be associated with NTG (p = .004). However, no statistically significant correlation was identified between this SNP and POAG (p = .88).

CONCLUSION

This study provides evidence of an association between the T allele of rs5443 and a reduced susceptibility NTG within the Russian population. These observations augment the comprehension of the genetic underpinnings of glaucoma and hold potential implications for the prospective development of targeted therapeutic interventions.

Hypertrophic Cardiomyopathy in Underrepresented Populations: Clinical and Genetic Landscape Based on a Russian Single-Center Cohort Study

Hypertrophic cardiomyopathy (HCM) is a common inherited cardiac disorder characterized by marked clinical and genetic heterogeneity. Ethnic groups underrepresented in studies may have distinctive characteristics. We sought to evaluate the clinical and genetic landscape of Russian HCM patients. A total of 193 patients (52% male; 95% Eastern Slavic origin; median age 56 years) were clinically evaluated, including genetic testing, and prospectively followed to document outcomes. As a result, 48% had obstructive HCM, 25% had HCM in family, 21% were asymptomatic, and 68% had comorbidities. During 2.8 years of follow-up, the all-cause mortality rate was 2.86%/year. A total of 5.7% received an implantable cardioverter-defibrillator (ICD), and 21% had septal reduction therapy. A sequencing analysis of 176 probands identified 64 causative variants in 66 patients (38%); recurrent variants were MYBPC3 p.Q1233* (8), MYBPC3 p.R346H (2), MYH7 p.A729P (2), TPM1 p.Q210R (3), and FLNC p.H1834Y (2); 10 were multiple variant carriers (5.7%); 5 had non-sarcomeric HCM, ALPK3, TRIM63, and FLNC. Thin filament variant carriers had a worse prognosis for heart failure (HR = 7.9, p = 0.007). In conclusion, in the Russian HCM population, the low use of ICD and relatively high mortality should be noted by clinicians; some distinct recurrent variants are suspected to have a founder effect; and family studies on some rare variants enriched worldwide knowledge in HCM.

Transferability of the PRS estimates for height and BMI obtained from the European ethnic groups to the Western Russian populations

Genetic data plays an increasingly important role in modern medicine. Decrease in the cost of sequencing with subsequent increase in imputation accuracy, and the accumulation of large amounts of high-quality genetic data enable the creation of polygenic risk scores (PRSs) to perform genotype-phenotype associations. The accuracy of phenotype prediction primarily depends on the overall trait heritability, Genome-wide association studies cohort size, and the similarity of genetic background between the base and the target cohort. Here we utilized 8,664 high coverage genomic samples collected across Russia by "Evogen", a Russian biomedical company, to evaluate the predictive power of PRSs based on summary statistics established on cohorts of European ancestry for basic phenotypic traits, namely height and BMI. We have demonstrated that the PRSs calculated for selected traits in three distinct Russian populations, recapitulate the predictive power from the original studies. This is evidence that GWAS summary statistics calculated on cohorts of European ancestry are transferable onto at least some ethnic groups in Russia.

The Pioneer Advantage: Filling the blank spots on the map of genome diversity in Europe

Documenting genome diversity is important for the local biomedical communities and instrumental in developing precision and personalized medicine. Currently, tens of thousands of whole-genome sequences from Europe are publicly available, but most of these represent populations of developed countries of Europe. The uneven distribution of the available data is further impaired by the lack of data sharing. Recent whole-genome studies in Eastern Europe, one in Ukraine and one in Russia, demonstrated that local genome diversity and population structure from Eastern Europe historically had not been fully represented. An unexpected wealth of genomic variation uncovered in these studies was not so much a consequence of high variation within their population, but rather due to the "pioneer advantage." We discovered more variants because we were the first to prospect in the Eastern European genome pool. This simple comparison underscores the importance of removing the remaining geographic genome deserts from the rest of the world map of the human genome diversity.

Peculiarity of Pomors of Onega Peninsula and Winter Coast in the genetic context of Northern Europe

The peculiarity of the Russian North gene pool has long become scientific fact, but has yet to receive informative explanation. Genetic drift cannot be the only contributing factor in the striking genetic differences between not only northern Russian populations and the southern ones, but among individual northern populations as well. Studying Russian North gene pools previously underrepresented in scientific literature may help understand this phenomenon. The work aimed to perform a subtotal study of the gene pool of the Arkhangelsk Oblast Pomors (Onega Coast, Summer Coast, the western fragment of the Winter Coast; n = 130) using a panel of 60 Y-chromosome SNP markers through multidimensional scaling and mapping of genetic distances. The frequencies of 14 identified haplogroups differ drastically in Pomor populations: haplogroups I1, R1a, and N3 each comprise a quarter of the total Pomor gene pool, I2-P37.2, and R1b each comprise about 8%, and the rest of the haplogroups are rare. The Onega Coast Pomors showed genetic similarity to a wide range of North-Eastern Europe Finnic-speaking populations, as well as to Russian populations with a strong pre-Slavic substratum. The Summer Coast Pomors are close to the Scandinavian gene pools, and the Winter Coast Pomors are similar only to specific Finn and Swede populations. None of the Pomor populations demonstrate genetic similarity with the Novgorod Oblast Russian populations, with which the origin of the Pomors is traditionally associated. The genetic distances between Pomor populations are so great, they are comparable to the general range of variability between the Eastern Slavic, Baltic, and Finno-Ugric peoples of the region. The reasons for such pronounced originality of Pomor populations presumably include, along with genetic drift, the gene pool of each population being underlied by a different pre-Slavic substrate, with later gene flows as an additional factor.

Identification of Pathogenic Variant Burden and Selection of Optimal Diagnostic Method Is a Way to Improve Carrier Screening for Autosomal Recessive Diseases

Cystic fibrosis, phenylketonuria, alpha-1 antitrypsin deficiency, and sensorineural hearing loss are among the most common autosomal recessive diseases, which require carrier screening. The evaluation of population allele frequencies (AF) of pathogenic variants in genes associated with these conditions and the choice of the best genotyping method are the necessary steps toward development and practical implementation of carrier-screening programs. We performed custom panel genotyping of 3821 unrelated participants from two Russian population representative samples and three patient groups using real-time polymerase chain reaction (PCR) and next generation sequencing (NGS). The custom panel included 115 known pathogenic variants in the CFTR, PAH, SERPINA1, and GJB2 genes. Overall, 38 variants were detected. The comparison of genotyping platforms revealed the following advantages of real-time PCR: relatively low cost, simple genotyping data analysis, and easier detection of large indels, while NGS showed better accuracy of variants identification and capability for detection of additional pathogenic variants in adjacent regions. A total of 23 variants had significant differences in estimated AF comparing with non-Finnish Europeans from gnomAD. This study provides new AF data for variants associated with the studied disorders and the comparison of genotyping methods for carrier screening.

Genotype imputation and polygenic score estimation in northwestern Russian population

Numerous studies demonstrated the lack of transferability of polygenic score (PGS) models across populations and the problem arising from unequal presentation of ancestries across genetic studies. However, even within European ancestry there are ethnic groups that are rarely presented in genetic studies. For instance, Russians, being one of the largest, diverse, and yet understudied group in Europe. In this study, we evaluated the reliability of genotype imputation for the Russian cohort by testing several commonly used imputation reference panels (e.g. HRC, 1000G, HGDP). HRC, in comparison with two other panels, showed the most accurate results based on both imputation accuracy and allele frequency concordance between masked and imputed genotypes. We built polygenic score models based on GWAS results from the UK biobank, measured the explained phenotypic variance in the Russian cohort attributed to polygenic scores for 11 phenotypes, collected in the clinic for each participant, and finally explored the role of allele frequency discordance between the UK biobank and the study cohort in the resulting PGS performance.

Urine TERT promoter mutations-based tumor DNA detection in patients with bladder cancer: A pilot study

Telomerase reverse transcriptase (TERT) promoter mutations are the most frequent genetic events in bladder cancer (BC). The aim of the present pilot study was to evaluate the diagnostic potential of urine TERT promoter mutations-based liquid biopsy in patients with an ongoing oncological process, as well as in post-resection patients at risk of BC recurrence. A total of 60 patients were enrolled, of whom 27 patients had histologically proven BC; 23 had no signs of BC (control group); and 10 patients underwent transurethral malignancy resection 3-6 months prior to urine donation ('second look' group). Urine TERT promoter mutations were detected using Droplet Digital PCR. Receiver operating characteristic curve analysis revealed significant diagnostic power of the present approach (area under the curve: -0.768). At the cut-off value of tumor DNA fraction 0.34%, the sensitivity and specificity were 55.56 and 100%, respectively. In the positive samples, tumor DNA fraction varied significantly from 0.59 to 48.77%. In the 'second look' group, tumor DNA was detected in 4/10 patients, highlighting the possibility of BC recurrence with its fraction ranging only from 0.90 to 6.61%. Therefore, urine TERT promoter mutations-based liquid biopsy appears to be a promising tool for BC diagnosis and surveillance. The main study will include recruitment of additional patients, extension of the mutation panel, prolonged follow-up of the post-resection patients, as well as screening of industrial workers exposed to specific carcinogens.

Targeted Sequencing of 242 Clinically Important Genes in the Russian Population From the Ivanovo Region

We performed a targeted sequencing of 242 clinically important genes mostly associated with cardiovascular diseases in a representative population sample of 1,658 individuals from the Ivanovo region northeast of Moscow. Approximately 11% of 11,876 detected variants were not found in the Single Nucleotide Polymorphism Database (dbSNP) or reported earlier in the Russian population. Most novel variants were singletons and doubletons in our sample, and virtually no novel alleles presumably specific for the Russian population were able to reach the frequencies above 0.1-0.2%. The overwhelming majority (99.3%) of variants detected in this study in three or more copies were shared with other populations. We found two dominant and seven recessive known pathogenic variants with allele frequencies significantly increased compared to those in the gnomAD non-Finnish Europeans. Of the 242 targeted genes, 28 were in the list of 59 genes for which the American College of Medical Genetics and Genomics (ACMG) recommended the reporting of incidental findings. Based on the number of variants detected in the sequenced subset of ACMG59 genes, we approximated the prevalence of known pathogenic and novel or rare protein-truncating variants in the complete set of ACMG59 genes in the Ivanovo population at 1.4 and 2.8%, respectively. We analyzed the available clinical data and observed the incomplete penetrance of known pathogenic variants in the 28 ACMG59 genes: only 1 individual out of 12 with such variants had the phenotype most likely related to the variant. When known pathogenic and novel or rare protein-truncating variants were considered together, the overall rate of confirmed phenotypes was about 19%, with maximum in the subset of novel protein-truncating variants. We report three novel protein truncating variants in APOB and one in MYH7 observed in individuals with hypobetalipoproteinemia and hypertrophic cardiomyopathy, respectively. Our results provide a valuable reference for the clinical interpretation of gene sequencing in Russian and other populations.

SLC6A3 (DAT1) as a Novel Candidate Biomarker Gene for Suicidal Behavior

It has been previously shown that the serotonin and dopamine neurotransmitter systems might influence the predisposition to suicidal behavior. This study aims to estimate the contribution of 11 polymorphisms in the genes SLC6A4 (5HTT), HTR1A, HTR2A, HTR1B, SLC6A3 (DAT1), DRD4, DRD2, COMT, and BDNF to suicidal behavior and severity of symptoms of depression and anxiety in the Russian population. The study was performed on 100 patients with repeated suicide attempts and 154 controls. We first found an association between SLC6A3 (DAT1) 40 bp VNTR locus and suicidal behavior. This association was significant; when using the codominant (p = 0.006), dominant (p = 0.001), overdominant (p = 0.004), and log-additive (p = 0.004) models, LL genotype played a protective role (OR = 0.48, 0.29–0.82, p = 0.005). Difference in the distribution of COMT rs4680 genotypes was significant in the codominant (p = 0.04), dominant (p = 0.013), and log-additive (p = 0.02) models, and AA genotype might protect against suicide (OR = 0.49, 0.26–0.91, p = 0.025). SLC6A4 5-HTTLPR + rs25531 locus was significant in the recessive model (p = 0.024), and also affected the severity of symptoms of depression (p = 0.044) and personal anxiety (p = 0.029). Our results suggest that allelic variants of SLC6A3, COMT, and SLC6A4 genes might be considered as risk factors for suicidal attempts.

Genome diversity in Ukraine

Background

The main goal of this collaborative effort is to provide genome-wide data for the previously underrepresented population in Eastern Europe, and to provide cross-validation of the data from genome sequences and genotypes of the same individuals acquired by different technologies. We collected 97 genome-grade DNA samples from consented individuals representing major regions of Ukraine that were consented for public data release. BGISEQ-500 sequence data and genotypes by an Illumina GWAS chip were cross-validated on multiple samples and additionally referenced to 1 sample that has been resequenced by Illumina NovaSeq6000 S4 at high coverage.

Results

The genome data have been searched for genomic variation represented in this population, and a number of variants have been reported: large structural variants, indels, copy number variations, single-nucletide polymorphisms, and microsatellites. To our knowledge, this study provides the largest to-date survey of genetic variation in Ukraine, creating a public reference resource aiming to provide data for medical research in a large understudied population.

Conclusions

Our results indicate that the genetic diversity of the Ukrainian population is uniquely shaped by evolutionary and demographic forces and cannot be ignored in future genetic and biomedical studies. These data will contribute a wealth of new information bringing forth a wealth of novel, endemic and medically related alleles.

Whole-exome sequencing provides insights into monogenic disease prevalence in Northwest Russia

BACKGROUND

Allele frequency data from large exome and genome aggregation projects such as the Genome Aggregation Database (gnomAD) are of ultimate importance to the interpretation of medical resequencing data. However, allele frequencies might significantly differ in poorly studied populations that are underrepresented in large-scale projects, such as the Russian population.

METHODS

In this work, we leveraged our access to a large dataset of 694 exome samples to analyze genetic variation in the Northwest Russia. We compared the spectrum of genetic variants to the dbSNP build 151, and made estimates of ClinVar-based autosomal recessive (AR) disease allele prevalence as compared to gnomAD r. 2.1.

RESULTS

An estimated 9.3% of discovered variants were not present in dbSNP. We report statistically significant overrepresentation of pathogenic variants for several Mendelian disorders, including phenylketonuria (PAH, rs5030858), Wilson's disease (ATP7B, rs76151636), factor VII deficiency (F7, rs36209567), kyphoscoliosis type of Ehlers-Danlos syndrome (FKBP14, rs542489955), and several other recessive pathologies. We also make primary estimates of monogenic disease incidence in the population, with retinal dystrophy, cystic fibrosis, and phenylketonuria being the most frequent AR pathologies.

CONCLUSION

Our observations demonstrate the utility of population-specific allele frequency data to the diagnosis of monogenic disorders using high-throughput technologies.