Analysis of Common SNPs across Continents Reveals Major Genomic Differences between Human Populations

Common alleles tend to be more ancient than rare alleles. These common SNPs appeared thousands of years ago and reflect intricate human evolution including various adaptations, admixtures, and migration events. Eighty-four thousand abundant region-specific alleles (ARSAs) that are common in one continent but absent in the rest of the world have been characterized by processing 3100 genomes from 230 populations. Also computed were 17,446 polymorphic sites with regional absence of common alleles (RACAs), which are widespread globally but absent in one region. A majority of these region-specific SNPs were found in Africa. America has the second greatest number of ARSAs (3348) and is even ahead of Europe (1911). Surprisingly, East Asia has the highest number of RACAs (10,524) and the lowest number of ARSAs (362). ARSAs and RACAs have distinct compositions of ancestral versus derived alleles in different geographical regions, reflecting their unique evolution. Genes associated with ARSA and RACA SNPs were identified and their functions were analyzed. The core 100 genes shared by multiple populations and associated with region-specific natural selection were examined. The largest part of them (42%) are related to the nervous system. ARSA and RACA SNPs are important for both association and human evolution studies.

CAPN3 c.1746-20C>G variant is hypomorphic for LGMD R1 calpain 3-related

The investigated intronic CAPN3 variant NM_000070.3:c.1746-20C>G occurs in the Central and Eastern Europe with a frequency of >1% and there are conflicting interpretations on its pathogenicity. We collected data on 14 patients carrying the CAPN3 c.1746-20C>G variant in trans position with another CAPN3 pathogenic/likely pathogenic variant. The patients compound heterozygous for the CAPN3 c.1746-20C>G variant presented a phenotype consistent with calpainopathy of mild/medium severity. This variant is most frequent in the North/West regions of Russia and may originate from that area. Molecular studies revealed that different splicing isoforms are produced in the muscle. We hypothesize that c.1746-20C>G is a hypomorphic variant with a reduction of RNA and protein expression and only individuals having a higher ratio of abnormal isoforms are affected. Reclassification of the CAPN3 variant c.1746-20C>G from variant with a conflicting interpretation of pathogenicity to hypomorphic variant explains many unidentified cases of limb girdle muscular dystrophy R1 calpain 3-related in Eastern and Central Europe.

Mapping genomic loci implicates genes and synaptic biology in schizophrenia

Schizophrenia has a heritability of 60-80%1, much of which is attributable to common risk alleles. Here, in a two-stage genome-wide association study of up to 76,755 individuals with schizophrenia and 243,649 control individuals, we report common variant associations at 287 distinct genomic loci. Associations were concentrated in genes that are expressed in excitatory and inhibitory neurons of the central nervous system, but not in other tissues or cell types. Using fine-mapping and functional genomic data, we identify 120 genes (106 protein-coding) that are likely to underpin associations at some of these loci, including 16 genes with credible causal non-synonymous or untranslated region variation. We also implicate fundamental processes related to neuronal function, including synaptic organization, differentiation and transmission. Fine-mapped candidates were enriched for genes associated with rare disruptive coding variants in people with schizophrenia, including the glutamate receptor subunit GRIN2A and transcription factor SP4, and were also enriched for genes implicated by such variants in neurodevelopmental disorders. We identify biological processes relevant to schizophrenia pathophysiology; show convergence of common and rare variant associations in schizophrenia and neurodevelopmental disorders; and provide a resource of prioritized genes and variants to advance mechanistic studies.

A Workflow for Selection of Single Nucleotide Polymorphic Markers for Studying of Genetics of Ischemic Stroke Outcomes

In this paper we propose a workflow for studying the genetic architecture of ischemic stroke outcomes. It develops further the candidate gene approach. The workflow is based on the animal model of brain ischemia, comparative genomics, human genomic variations, and algorithms of selection of tagging single nucleotide polymorphisms (tagSNPs) in genes which expression was changed after ischemic stroke. The workflow starts from a set of rat genes that changed their expression in response to brain ischemia and results in a set of tagSNPs, which represent other SNPs in the human genes analyzed and influenced on their expression as well.

Analysis of shared heritability in common disorders of the brain

Disorders of the brain can exhibit considerable epidemiological comorbidity and often share symptoms, provoking debate about their etiologic overlap. We quantified the genetic sharing of 25 brain disorders from genome-wide association studies of 265,218 patients and 784,643 control participants and assessed their relationship to 17 phenotypes from 1,191,588 individuals. Psychiatric disorders share common variant risk, whereas neurological disorders appear more distinct from one another and from the psychiatric disorders. We also identified significant sharing between disorders and a number of brain phenotypes, including cognitive measures. Further, we conducted simulations to explore how statistical power, diagnostic misclassification, and phenotypic heterogeneity affect genetic correlations. These results highlight the importance of common genetic variation as a risk factor for brain disorders and the value of heritability-based methods in understanding their etiology.

Estimation of Genetic Correlation via Linkage Disequilibrium Score Regression and Genomic Restricted Maximum Likelihood

Genetic correlation is a key population parameter that describes the shared genetic architecture of complex traits and diseases. It can be estimated by current state-of-art methods, i.e., linkage disequilibrium score regression (LDSC) and genomic restricted maximum likelihood (GREML). The massively reduced computing burden of LDSC compared to GREML makes it an attractive tool, although the accuracy (i.e., magnitude of standard errors) of LDSC estimates has not been thoroughly studied. In simulation, we show that the accuracy of GREML is generally higher than that of LDSC. When there is genetic heterogeneity between the actual sample and reference data from which LD scores are estimated, the accuracy of LDSC decreases further. In real data analyses estimating the genetic correlation between schizophrenia (SCZ) and body mass index, we show that GREML estimates based on ∼150,000 individuals give a higher accuracy than LDSC estimates based on ∼400,000 individuals (from combined meta-data). A GREML genomic partitioning analysis reveals that the genetic correlation between SCZ and height is significantly negative for regulatory regions, which whole genome or LDSC approach has less power to detect. We conclude that LDSC estimates should be carefully interpreted as there can be uncertainty about homogeneity among combined meta-datasets. We suggest that any interesting findings from massive LDSC analysis for a large number of complex traits should be followed up, where possible, with more detailed analyses with GREML methods, even if sample sizes are lesser.

A New Baltic Population-Specific Human Genetic Marker in the PMCA4 Gene

OBJECTIVES

The PMCA gene family consists of 4 genes and at least 21 splice variants; among these, the Ca2+ ATPase 4 (PMCA4) gene encodes a plasma membrane protein abundantly expressed in several tissues, including the kidney, heart, and sperm. Knockout of PMCA4 causes infertility due to immotile sperm in mouse models. We therefore investigated variants in this gene for potential association with infertility in groups of Estonian (n = 191) and Latvian (n = 92) men with reduced sperm motility.

METHODS

All exons, exon-intron boundaries, 5' and 3' untranslated regions, and the promoter region of the PMCA4 gene were analysed by direct sequencing for a group of Estonian infertile men. Genotyping of guanine and adenine alleles of rs147729934 was performed, using a custom-designed TaqMan® probe for a group of Latvian infertile men as well as additional groups from Latvia and several groups of people with proven ethnicity from the Baltic region.

RESULTS

Although we did not identify any significant associations between variants in the gene and infertility, our results indicated that in all studied Latvian and Estonian groups the adenine allele of the variant rs147729934 was present at a higher frequency than expected. Analysis of additional samples indicated that the adenine allele of rs147729934 likely originated once in the modern-day Baltic or western Russia area, as the frequency of the minor adenine allele observed in this region is remarkably higher than that in the general European population.

CONCLUSIONS

Our results revealed no significant difference in frequencies of genetic variants in PMCA4 gene between men with normal and those with reduced sperm motility. The adenine allele of the variant rs147729934 is potentially an informative tool for future population studies concerning ancient Baltic and Finno-Ugric history.

Reconstructing genetic history of Siberian and Northeastern European populations

Siberia and Northwestern Russia are home to over 40 culturally and linguistically diverse indigenous ethnic groups, yet genetic variation and histories of peoples from this region are largely uncharacterized. We present deep whole-genome sequencing data (∼38×) from 28 individuals belonging to 14 distinct indigenous populations from that region. We combined these data sets with additional 32 modern-day and 46 ancient human genomes to reconstruct genetic histories of several indigenous Northern Eurasian populations. We found that Siberian and East Asian populations shared 38% of their ancestry with a 45,000-yr-old Ust'-Ishim individual who was previously believed to have no modern-day descendants. Western Siberians trace 57% of their ancestry to ancient North Eurasians, represented by the 24,000-yr-old Siberian Mal'ta boy MA-1. Eastern Siberian populations formed a distinct sublineage that separated from other East Asian populations ∼10,000 yr ago. In addition, we uncovered admixtures between Siberians and Eastern European hunter-gatherers from Samara, Karelia, Hungary, and Sweden (from 8000-6600 yr ago); Yamnaya people (5300-4700 yr ago); and modern-day Northeastern Europeans. Our results provide new insights into genetic histories of Siberian and Northeastern European populations and evidence of ancient gene flow from Siberia into Europe.

Evidence for Genetic Overlap Between Schizophrenia and Age at First Birth in Women

IMPORTANCE

A recently published study of national data by McGrath et al in 2014 showed increased risk of schizophrenia (SCZ) in offspring associated with both early and delayed parental age, consistent with a U-shaped relationship. However, it remains unclear if the risk to the child is due to psychosocial factors associated with parental age or if those at higher risk for SCZ tend to have children at an earlier or later age.

OBJECTIVE

To determine if there is a genetic association between SCZ and age at first birth (AFB) using genetically informative but independently ascertained data sets.

DESIGN, SETTING, AND PARTICIPANTS

This investigation used multiple independent genome-wide association study data sets. The SCZ sample comprised 18 957 SCZ cases and 22 673 controls in a genome-wide association study from the second phase of the Psychiatric Genomics Consortium, and the AFB sample comprised 12 247 genotyped women measured for AFB from the following 4 community cohorts: Estonia (Estonian Genome Center Biobank, University of Tartu), the Netherlands (LifeLines Cohort Study), Sweden (Swedish Twin Registry), and the United Kingdom (TwinsUK). Schizophrenia genetic risk for each woman in the AFB community sample was estimated using genetic effects inferred from the SCZ genome-wide association study.

MAIN OUTCOMES AND MEASURES

We tested if SCZ genetic risk was a significant predictor of response variables based on published polynomial functions that described the relationship between maternal age and SCZ risk in offspring in Denmark. We substituted AFB for maternal age in these functions, one of which was corrected for the age of the father, and found that the fit was superior for the model without adjustment for the father's age.

RESULTS

We observed a U-shaped relationship between SCZ risk and AFB in the community cohorts, consistent with the previously reported relationship between SCZ risk in offspring and maternal age when not adjusted for the age of the father. We confirmed that SCZ risk profile scores significantly predicted the response variables (coefficient of determination R2 = 1.1E-03, P = 4.1E-04), reflecting the published relationship between maternal age and SCZ risk in offspring by McGrath et al in 2014.

CONCLUSIONS AND RELEVANCE

This study provides evidence for a significant overlap between genetic factors associated with risk of SCZ and genetic factors associated with AFB. It has been reported that SCZ risk associated with increased maternal age is explained by the age of the father and that de novo mutations that occur more frequently in the germline of older men are the underlying causal mechanism. This explanation may need to be revised if, as suggested herein and if replicated in future studies, there is also increased genetic risk of SCZ in older mothers.

Modeling Linkage Disequilibrium Increases Accuracy of Polygenic Risk Scores

Polygenic risk scores have shown great promise in predicting complex disease risk and will become more accurate as training sample sizes increase. The standard approach for calculating risk scores involves linkage disequilibrium (LD)-based marker pruning and applying a p value threshold to association statistics, but this discards information and can reduce predictive accuracy. We introduce LDpred, a method that infers the posterior mean effect size of each marker by using a prior on effect sizes and LD information from an external reference panel. Theory and simulations show that LDpred outperforms the approach of pruning followed by thresholding, particularly at large sample sizes. Accordingly, predicted R(2) increased from 20.1% to 25.3% in a large schizophrenia dataset and from 9.8% to 12.0% in a large multiple sclerosis dataset. A similar relative improvement in accuracy was observed for three additional large disease datasets and for non-European schizophrenia samples. The advantage of LDpred over existing methods will grow as sample sizes increase.

Partitioning heritability of regulatory and cell-type-specific variants across 11 common diseases

Regulatory and coding variants are known to be enriched with associations identified by genome-wide association studies (GWASs) of complex disease, but their contributions to trait heritability are currently unknown. We applied variance-component methods to imputed genotype data for 11 common diseases to partition the heritability explained by genotyped SNPs (hg(2)) across functional categories (while accounting for shared variance due to linkage disequilibrium). Extensive simulations showed that in contrast to current estimates from GWAS summary statistics, the variance-component approach partitions heritability accurately under a wide range of complex-disease architectures. Across the 11 diseases DNaseI hypersensitivity sites (DHSs) from 217 cell types spanned 16% of imputed SNPs (and 24% of genotyped SNPs) but explained an average of 79% (SE = 8%) of hg(2) from imputed SNPs (5.1× enrichment; p = 3.7 × 10(-17)) and 38% (SE = 4%) of hg(2) from genotyped SNPs (1.6× enrichment, p = 1.0 × 10(-4)). Further enrichment was observed at enhancer DHSs and cell-type-specific DHSs. In contrast, coding variants, which span 1% of the genome, explained <10% of hg(2) despite having the highest enrichment. We replicated these findings but found no significant contribution from rare coding variants in independent schizophrenia cohorts genotyped on GWAS and exome chips. Our results highlight the value of analyzing components of heritability to unravel the functional architecture of common disease.

Genetic characterization of northeastern Italian population isolates in the context of broader European genetic diversity

Population genetic studies on European populations have highlighted Italy as one of genetically most diverse regions. This is possibly due to the country's complex demographic history and large variability in terrain throughout the territory. This is the reason why Italy is enriched for population isolates, Sardinia being the best-known example. As the population isolates have a great potential in disease-causing genetic variants identification, we aimed to genetically characterize a region from northeastern Italy, which is known for isolated communities. Total of 1310 samples, collected from six geographically isolated villages, were genotyped at >145 000 single-nucleotide polymorphism positions. Newly genotyped data were analyzed jointly with the available genome-wide data sets of individuals of European descent, including several population isolates. Despite the linguistic differences and geographical isolation the village populations still show the greatest genetic similarity to other Italian samples. The genetic isolation and small effective population size of the village populations is manifested by higher levels of genomic homozygosity and elevated linkage disequilibrium. These estimates become even more striking when the detected substructure is taken into account. The observed level of genetic isolation in Friuli-Venezia Giulia region is more extreme according to several measures of isolation compared with Sardinians, French Basques and northern Finns, thus proving the status of an isolate.

Pharmacogenomics of cisplatin-based chemotherapy in ovarian cancer patients of different ethnic origins

AIM

There is a substantial difference between Asians and Caucasians in their reaction to platinum drugs. To determine whether population-related genomics contribute to differences in patient outcomes, pharmacogenomic relevance of polymorphisms in some genes, the protein activities of which may affect aspects of cisplatin metabolism, were evaluated.

PATIENTS & METHODS

Nineteen polymorphisms in ten genes were tested for correlations with the efficacy and toxicity of a cisplatin-cyclophosphamide regimen in Yakut and Russian patients with ovarian cancer.

RESULTS

The CYP2E1 7632T>A polymorphism was associated with progression-free survival (p = 0.015) in Yakuts. In Russians, progression-free survival was correlated with the GSTP1 Ile105Val polymorphism (p = 0.004). Yakut patients with the GSTT1-null genotype had a higher risk for nephrotoxicity. By contrast, in the Russian group, nephrotoxicity was more frequent among patients with heterozygous ERCC1 genotypes. Severe emesis in Yakuts was independently associated with two polymorphisms in the CYP2E1 gene but in Russians, it was more common in patients with the GSTT1-null genotype. Differences in genotypic correlations with anemia were also observed.

CONCLUSION

Significant differences in genotype distribution between Russian and Yakut women were observed for ten of the 19 polymorphisms, but none of them seemed to be a clear casual candidate and further studies involving more markers are required.

Genetic structure of Europeans: a view from the North-East

Using principal component (PC) analysis, we studied the genetic constitution of 3,112 individuals from Europe as portrayed by more than 270,000 single nucleotide polymorphisms (SNPs) genotyped with the Illumina Infinium platform. In cohorts where the sample size was >100, one hundred randomly chosen samples were used for analysis to minimize the sample size effect, resulting in a total of 1,564 samples. This analysis revealed that the genetic structure of the European population correlates closely with geography. The first two PCs highlight the genetic diversity corresponding to the northwest to southeast gradient and position the populations according to their approximate geographic origin. The resulting genetic map forms a triangular structure with a) Finland, b) the Baltic region, Poland and Western Russia, and c) Italy as its vertexes, and with d) Central- and Western Europe in its centre. Inter- and intra- population genetic differences were quantified by the inflation factor lambda (λ) (ranging from 1.00 to 4.21), fixation index (F_st) (ranging from 0.000 to 0.023), and by the number of markers exhibiting significant allele frequency differences in pair-wise population comparisons. The estimated lambda was used to assess the real diminishing impact to association statistics when two distinct populations are merged directly in an analysis. When the PC analysis was confined to the 1,019 Estonian individuals (0.1% of the Estonian population), a fine structure emerged that correlated with the geography of individual counties. With at least two cohorts available from several countries, genetic substructures were investigated in Czech, Finnish, German, Estonian and Italian populations. Together with previously published data, our results allow the creation of a comprehensive European genetic map that will greatly facilitate inter-population genetic studies including genome wide association studies (GWAS).

Development of a single tube 640-plex genotyping method for detection of nucleic acid variations on microarrays

Detection of DNA sequence variation is critical to biomedical applications, including disease genetic identification, diagnosis and treatment, drug discovery and forensic analysis. Here, we describe an arrayed primer extension-based genotyping method (APEX-2) that allows multiplex (640-plex) DNA amplification and detection of single nucleotide polymorphisms (SNPs) and mutations on microarrays via four-color single-base primer extension. The founding principle of APEX-2 multiplex PCR requires two oligonucleotides per SNP/mutation to generate amplicons containing the position of interest. The same oligonucleotides are then subsequently used as immobilized single-base extension primers on a microarray. The method described here is ideal for SNP or mutation detection analysis, molecular diagnostics and forensic analysis. This robust genetic test has minimal requirements: two primers, two spots on the microarray and a low cost four-color detection system for the targeted site; and provides an advantageous alternative to high-density platforms and low-density detection systems.

Regional differences in the genetic variability of Finno-Ugric speaking Komi populations

The Komi (Komi-Zyryan) people are one of the most numerous ethnic groups belonging to the Finno-Ugric linguistic community. They occupy an extensive territory in north Russia to the west of the Ural Mountains, in the northeast of the East European Plain. This is an area of long-term interactions between Europeans and North Asians. Genetic variability was evaluated in two geographically distinct populations, the Izhemski and Priluzski Komi. We searched for polymorphisms of the TP53 gene (a 16-bp duplication in intron 3 and three RFLPs: for Bsh1236I at codon 72, for MspI in intron 6, and for BamHI in the 3' flanking region) and for variable number tandem repeat (VNTR) polymorphisms of locus D1S80 and of the 3' untranslated region of the gene for apolipoprotein B (ApoB). Some data from our previous studies of TP53, 3'ApoB, and D1S80 variability were involved in the comparison of Komi with other Eastern European populations. Multidimensional scaling analysis of genetic distances was used for the evaluation of genetic relationships between populations. The results revealed some affinity between Priluzski Komi and Eastern Slavonic populations, and significant segregation of Izhemski Komi from other ethnic groups studied. The unique genetic features of Izhemski Komi may have been determined by their ethnogenesis or the pressure of environmental factors, such as special nutrition and adaptation to extreme climatic conditions.

Mitochondrial DNA Variations in Russian and Belorussian Populations

The sequence of the first hypervariable segment (HVS-I) of mitochondrial DNA (mtDNA) was determined in 251 individuals from three eastern Slavonic populations, two Russian and one Belorussian. Within HVS-I, 78 polymorphic positions were revealed. Within-population diversity of HVS-I varies slightly among three samples; its estimates do not differ strongly from those for European populations. Haplotype diversity for three populations calculated in this study is 0.949; mean pairwise differences estimate is 3.59. To assign mtDNA sequences to major phylogenetic clusters, haplogroup-specific restriction polymorphisms were selectively typed in most samples. The haplogroup distribution in the total Eastern Slavonic sample is similar to that reported for the European sample. However, the separate consideration of three Slavonic samples reveals the complicated structure of the mitochondrial gene pool in the Eastern European area. Data of this study support the proposed model of the origin of modern Eastern Slavs, which implies the admixture of ancient Slavonic tribes with pre-Slavonic populations of Eastern Europe. These data should contribute to general studies of mitochondrial DNA variations in Europe.