Polygenic risk score adds to a clinical risk score in the prediction of cardiovascular disease in a clinical setting.

BACKGROUND AND AIMS: A cardiovascular disease polygenic risk score (CVD-PRS) can stratify individuals into different categories of cardiovascular risk, but whether the addition of a CVD-PRS to clinical risk scores improves the identification of individuals at increased risk in a real-world clinical setting is unknown. METHODS: The Genetics and the Vascular Health Check Study (GENVASC) was embedded within the UK National Health Service Health Check (NHSHC) programme which invites individuals between 40-74 years of age without known CVD to attend an assessment in a UK general practice where CVD risk factors are measured and a CVD risk score (QRISK2) is calculated. Between 2012-2020, 44,141 individuals (55.7% females, 15.8% non-white) who attended an NHSHC in 147 participating practices across two counties in England were recruited and followed. When 195 individuals (cases) had suffered a major CVD event (CVD death, myocardial infarction or acute coronary syndrome, coronary revascularisation, stroke), 396 propensity-matched controls with a similar risk profile were identified, and a nested case-control genetic study undertaken to see if the addition of a CVD-PRS to QRISK2 in the form of an integrated risk tool (IRT) combined with QRISK2 would have identified more individuals at the time of their NHSHC as at high risk (QRISK2 10-year CVD risk of ≥10%), compared with QRISK2 alone. RESULTS: The distribution of the standardised CVD-PRS was significantly different in cases compared with controls (cases mean score .32; controls, -.18, P = 8.28×10-9). QRISK2 identified 61.5% (95% confidence interval [CI]: 54.3%-68.4%) of individuals who subsequently developed a major CVD event as being at high risk at their NHSHC, while the combination of QRISK2 and IRT identified 68.7% (95% CI: 61.7%-75.2%), a relative increase of 11.7% (P = 1×10-4). The odds ratio (OR) of being up-classified was 2.41 (95% CI: 1.03-5.64, P = .031) for cases compared with controls. In individuals aged 40-54 years, QRISK2 identified 26.0% (95% CI: 16.5%-37.6%) of those who developed a major CVD event, while the combination of QRISK2 and IRT identified 38.4% (95% CI: 27.2%-50.5%), indicating a stronger relative increase of 47.7% in the younger age group (P = .001). The combination of QRISK2 and IRT increased the proportion of additional cases identified similarly in women as in men, and in non-white ethnicities compared with white ethnicity. The findings were similar when the CVD-PRS was added to the atherosclerotic cardiovascular disease pooled cohort equations (ASCVD-PCE) or SCORE2 clinical scores. CONCLUSIONS: In a clinical setting, the addition of genetic information to clinical risk assessment significantly improved the identification of individuals who went on to have a major CVD event as being at high risk, especially among younger individuals. The findings provide important real-world evidence of the potential value of implementing a CVD-PRS into health systems.

Dystonia genes functionally converge in specific neurons and share neurobiology with psychiatric disorders.

Dystonia is a neurological disorder characterized by sustained or intermittent muscle contractions causing abnormal movements and postures, often occurring in absence of any structural brain abnormality. Psychiatric comorbidities, including anxiety, depression, obsessive-compulsive disorder and schizophrenia, are frequent in patients with dystonia. While mutations in a fast-growing number of genes have been linked to Mendelian forms of dystonia, the cellular, anatomical, and molecular basis remains unknown for most genetic forms of dystonia, as does its genetic and biological relationship to neuropsychiatric disorders. Here we applied an unbiased systems-biology approach to explore the cellular specificity of all currently known dystonia-associated genes, predict their functional relationships, and test whether dystonia and neuropsychiatric disorders share a genetic relationship. To determine the cellular specificity of dystonia-associated genes in the brain, single-nuclear transcriptomic data derived from mouse brain was used together with expression-weighted cell-type enrichment. To identify functional relationships among dystonia-associated genes, we determined the enrichment of these genes in co-expression networks constructed from 10 human brain regions. Stratified linkage-disequilibrium score regression was used to test whether co-expression modules enriched for dystonia-associated genes significantly contribute to the heritability of anxiety, major depressive disorder, obsessive-compulsive disorder, schizophrenia, and Parkinson's disease. Dystonia-associated genes were significantly enriched in adult nigral dopaminergic neurons and striatal medium spiny neurons. Furthermore, 4 of 220 gene co-expression modules tested were significantly enriched for the dystonia-associated genes. The identified modules were derived from the substantia nigra, putamen, frontal cortex, and white matter, and were all significantly enriched for genes associated with synaptic function. Finally, we demonstrate significant enrichments of the heritability of major depressive disorder, obsessive-compulsive disorder and schizophrenia within the putamen and white matter modules, and a significant enrichment of the heritability of Parkinson's disease within the substantia nigra module. In conclusion, multiple dystonia-associated genes interact and contribute to pathogenesis likely through dysregulation of synaptic signalling in striatal medium spiny neurons, adult nigral dopaminergic neurons and frontal cortical neurons. Furthermore, the enrichment of the heritability of psychiatric disorders in the co-expression modules enriched for dystonia-associated genes indicates that psychiatric symptoms associated with dystonia are likely to be intrinsic to its pathophysiology.

Recursive splicing in long vertebrate genes.

It is generally believed that splicing removes introns as single units from precursor messenger RNA transcripts. However, some long Drosophila melanogaster introns contain a cryptic site, known as a recursive splice site (RS-site), that enables a multi-step process of intron removal termed recursive splicing. The extent to which recursive splicing occurs in other species and its mechanistic basis have not been examined. Here we identify highly conserved RS-sites in genes expressed in the mammalian brain that encode proteins functioning in neuronal development. Moreover, the RS-sites are found in some of the longest introns across vertebrates. We find that vertebrate recursive splicing requires initial definition of an 'RS-exon' that follows the RS-site. The RS-exon is then excluded from the dominant mRNA isoform owing to competition with a reconstituted 5' splice site formed at the RS-site after the first splicing step. Conversely, the RS-exon is included when preceded by cryptic promoters or exons that fail to reconstitute an efficient 5' splice site. Most RS-exons contain a premature stop codon such that their inclusion can decrease mRNA stability. Thus, by establishing a binary splicing switch, RS-sites demarcate different mRNA isoforms emerging from long genes by coupling cryptic elements with inclusion of RS-exons.

The Genetic Legacy of Zoroastrianism in Iran and India: Insights into Population Structure, Gene Flow, and Selection.

Zoroastrianism is one of the oldest extant religions in the world, originating in Persia (present-day Iran) during the second millennium BCE. Historical records indicate that migrants from Persia brought Zoroastrianism to India, but there is debate over the timing of these migrations. Here we present genome-wide autosomal, Y chromosome, and mitochondrial DNA data from Iranian and Indian Zoroastrians and neighboring modern-day Indian and Iranian populations and conduct a comprehensive genome-wide genetic analysis in these groups. Using powerful haplotype-based techniques, we find that Zoroastrians in Iran and India have increased genetic homogeneity relative to other sampled groups in their respective countries, consistent with their current practices of endogamy. Despite this, we infer that Indian Zoroastrians (Parsis) intermixed with local groups sometime after their arrival in India, dating this mixture to 690-1390 CE and providing strong evidence that Iranian Zoroastrian ancestry was maintained primarily through the male line. By making use of the rich information in DNA from ancient human remains, we also highlight admixture in the ancestors of Iranian Zoroastrians dated to 570 BCE-746 CE, older than admixture seen in any other sampled Iranian group, consistent with a long-standing isolation of Zoroastrians from outside groups. Finally, we report results, and challenges, from a genome-wide scan to identify genomic regions showing signatures of positive selection in present-day Zoroastrians that might correlate to the prevalence of particular diseases among these communities.

Transcriptomic and genetic analyses reveal potential causal drivers for intractable partial epilepsy.

Mesial temporal lobe epilepsy with hippocampal sclerosis represents the most common epilepsy syndrome in adult patients with medically intractable partial epilepsy. Mesial temporal lobe epilepsy is usually regarded as a polygenic and complex disorder, still poorly understood but probably caused and perpetuated by dysregulation of numerous biological networks and cellular functions. The study of gene expression changes by single nucleotide polymorphisms in regulatory elements (expression quantitative trait loci, eQTLs) has been shown to be a powerful complementary approach to the detection and understanding of risk loci by genome-wide association studies. We performed a whole (gene and exon-level) transcriptome analysis on cortical tissue samples (Brodmann areas 20 and 21) from 86 patients with mesial temporal lobe epilepsy with hippocampal sclerosis and 75 neurologically healthy controls. Genome-wide genotyping data from the same individuals (patients and controls) were analysed and paired with the transcriptome data. We report potential epilepsy-risk eQTLs, some of which are specific to tissue from patients with mesial temporal lobe epilepsy with hippocampal sclerosis. We also found large transcriptional and splicing deregulation in mesial temporal lobe epilepsy with hippocampal sclerosis tissue as well as gene networks involving neuronal and glial mechanisms that provide new insights into the cause and maintenance of the seizures. These data (available via the 'Seizubraineac' web-tool resource, www.seizubraineac.org) will facilitate the identification of new therapeutic targets and biomarkers as well as genetic risk variants that could influence epilepsy and pharmacoresistance.

The impact of donor and recipient common clinical and genetic variation on estimated glomerular filtration rate in a European renal transplant population.

Genetic variation across the human leukocyte antigen loci is known to influence renal-transplant outcome. However, the impact of genetic variation beyond the human leukocyte antigen loci is less clear. We tested the association of common genetic variation and clinical characteristics, from both the donor and recipient, with posttransplant eGFR at different time-points, out to 5 years posttransplantation. We conducted GWAS meta-analyses across 10 844 donors and recipients from five European ancestry cohorts. We also analyzed the impact of polygenic risk scores (PRS), calculated using genetic variants associated with nontransplant eGFR, on posttransplant eGFR. PRS calculated using the recipient genotype alone, as well as combined donor and recipient genotypes were significantly associated with eGFR at 1-year posttransplant. Thirty-two percent of the variability in eGFR at 1-year posttransplant was explained by our model containing clinical covariates (including weights for death/graft-failure), principal components and combined donor-recipient PRS, with 0.3% contributed by the PRS. No individual genetic variant was significantly associated with eGFR posttransplant in the GWAS. This is the first study to examine PRS, composed of variants that impact kidney function in the general population, in a posttransplant context. Despite PRS being a significant predictor of eGFR posttransplant, the effect size of common genetic factors is limited compared to clinical variables.

Transethnic differences in GWAS signals: A simulation study.

Genome-wide association studies (GWASs) have allowed researchers to identify thousands of single nucleotide polymorphisms (SNPs) and other variants associated with particular complex traits. Previous studies have reported differences in the strength and even the direction of GWAS signals across different populations. These differences could be due to a combination of (1) lack of power, (2) allele frequency differences, (3) linkage disequilibrium (LD) differences, and (4) true differences in causal variant effect sizes. To determine whether properties (1)-(3) on their own might be sufficient to explain the patterns previously noted in strong GWAS signals, we simulated case-control data of European, Asian and African ancestry, applying realistic allele frequencies and LD from 1000 Genomes data but enforcing equal causal effect sizes across populations. Much of the observed differences in strong GWAS signals could indeed be accounted for by allele frequency and LD differences, enhanced by the Euro-centric SNP bias and lower SNP coverage found in older GWAS panels. While we cannot rule out a role for true transethnic effect size differences, our results suggest that strong causal effects may be largely shared among human populations, motivating the use of transethnic data for fine-mapping.

Long- and short-term outcomes in renal allografts with deceased donors: A large recipient and donor genome-wide association study.

Improvements in immunosuppression have modified short-term survival of deceased-donor allografts, but not their rate of long-term failure. Mismatches between donor and recipient HLA play an important role in the acute and chronic allogeneic immune response against the graft. Perfect matching at clinically relevant HLA loci does not obviate the need for immunosuppression, suggesting that additional genetic variation plays a critical role in both short- and long-term graft outcomes. By combining patient data and samples from supranational cohorts across the United Kingdom and European Union, we performed the first large-scale genome-wide association study analyzing both donor and recipient DNA in 2094 complete renal transplant-pairs with replication in 5866 complete pairs. We studied deceased-donor grafts allocated on the basis of preferential HLA matching, which provided some control for HLA genetic effects. No strong donor or recipient genetic effects contributing to long- or short-term allograft survival were found outside the HLA region. We discuss the implications for future research and clinical application.

Pooled sequencing of 531 genes in inflammatory bowel disease identifies an associated rare variant in BTNL2 and implicates other immune related genes.

The contribution of rare coding sequence variants to genetic susceptibility in complex disorders is an important but unresolved question. Most studies thus far have investigated a limited number of genes from regions which contain common disease associated variants. Here we investigate this in inflammatory bowel disease by sequencing the exons and proximal promoters of 531 genes selected from both genome-wide association studies and pathway analysis in pooled DNA panels from 474 cases of Crohn's disease and 480 controls. 80 variants with evidence of association in the sequencing experiment or with potential functional significance were selected for follow up genotyping in 6,507 IBD cases and 3,064 population controls. The top 5 disease associated variants were genotyped in an extension panel of 3,662 IBD cases and 3,639 controls, and tested for association in a combined analysis of 10,147 IBD cases and 7,008 controls. A rare coding variant p.G454C in the BTNL2 gene within the major histocompatibility complex was significantly associated with increased risk for IBD (p = 9.65x10-10, OR = 2.3[95% CI = 1.75-3.04]), but was independent of the known common associated CD and UC variants at this locus. Rare (<1%) and low frequency (1-5%) variants in 3 additional genes showed suggestive association (p<0.005) with either an increased risk (ARIH2 c.338-6C>T) or decreased risk (IL12B p.V298F, and NICN p.H191R) of IBD. These results provide additional insights into the involvement of the inhibition of T cell activation in the development of both sub-phenotypes of inflammatory bowel disease. We suggest that although rare coding variants may make a modest overall contribution to complex disease susceptibility, they can inform our understanding of the molecular pathways that contribute to pathogenesis.

Network-Informed Gene Ranking Tackles Genetic Heterogeneity in Exome-Sequencing Studies of Monogenic Disease.

Genetic heterogeneity presents a significant challenge for the identification of monogenic disease genes. Whole-exome sequencing generates a large number of candidate disease-causing variants and typical analyses rely on deleterious variants being observed in the same gene across several unrelated affected individuals. This is less likely to occur for genetically heterogeneous diseases, making more advanced analysis methods necessary. To address this need, we present HetRank, a flexible gene-ranking method that incorporates interaction network data. We first show that different genes underlying the same monogenic disease are frequently connected in protein interaction networks. This motivates the central premise of HetRank: those genes carrying potentially pathogenic variants and whose network neighbors do so in other affected individuals are strong candidates for follow-up study. By simulating 1,000 exome sequencing studies (20,000 exomes in total), we model varying degrees of genetic heterogeneity and show that HetRank consistently prioritizes more disease-causing genes than existing analysis methods. We also demonstrate a proof-of-principle application of the method to prioritize genes causing Adams-Oliver syndrome, a genetically heterogeneous rare disease. An implementation of HetRank in R is available via the Website http://sourceforge.net/p/hetrank/.

Assessing models for genetic prediction of complex traits: a comparison of visualization and quantitative methods.

BACKGROUND: In silico models have recently been created in order to predict which genetic variants are more likely to contribute to the risk of a complex trait given their functional characteristics. However, there has been no comprehensive review as to which type of predictive accuracy measures and data visualization techniques are most useful for assessing these models. METHODS: We assessed the performance of the models for predicting risk using various methodologies, some of which include: receiver operating characteristic (ROC) curves, histograms of classification probability, and the novel use of the quantile-quantile plot. These measures have variable interpretability depending on factors such as whether the dataset is balanced in terms of numbers of genetic variants classified as risk variants versus those that are not. RESULTS: We conclude that the area under the curve (AUC) is a suitable starting place, and for models with similar AUCs, violin plots are particularly useful for examining the distribution of the risk scores.

Resolving the polymorphism-in-probe problem is critical for correct interpretation of expression QTL studies.

Polymorphisms in the target mRNA sequence can greatly affect the binding affinity of microarray probe sequences, leading to false-positive and false-negative expression quantitative trait locus (QTL) signals with any other polymorphisms in linkage disequilibrium. We provide the most complete solution to this problem, by using the latest genome and exome sequence reference data to identify almost all common polymorphisms (frequency >1% in Europeans) in probe sequences for two commonly used microarray panels (the gene-based Illumina Human HT12 array, which uses 50-mer probes, and exon-based Affymetrix Human Exon 1.0 ST array, which uses 25-mer probes). We demonstrate the impact of this problem using cerebellum and frontal cortex tissues from 438 neuropathologically normal individuals. We find that although only a small proportion of the probes contain polymorphisms, they account for a large proportion of apparent expression QTL signals, and therefore result in many false signals being declared as real. We find that the polymorphism-in-probe problem is insufficiently controlled by previous protocols, and illustrate this using some notable false-positive and false-negative examples in MAPT and PRICKLE1 that can be found in many eQTL databases. We recommend that both new and existing eQTL data sets should be carefully checked in order to adequately address this issue.

Genetic signatures reveal high-altitude adaptation in a set of ethiopian populations.

The Tibetan and Andean Plateaus and Ethiopian highlands are the largest regions to have long-term high-altitude residents. Such populations are exposed to lower barometric pressures and hence atmospheric partial pressures of oxygen. Such "hypobaric hypoxia" may limit physical functional capacity, reproductive health, and even survival. As such, selection of genetic variants advantageous to hypoxic adaptation is likely to have occurred. Identifying signatures of such selection is likely to help understanding of hypoxic adaptive processes. Here, we seek evidence of such positive selection using five Ethiopian populations, three of which are from high-altitude areas in Ethiopia. As these populations may have been recipients of Eurasian gene flow, we correct for this admixture. Using single-nucleotide polymorphism genotype data from multiple populations, we find the strongest signal of selection in BHLHE41 (also known as DEC2 or SHARP1). Remarkably, a major role of this gene is regulation of the same hypoxia response pathway on which selection has most strikingly been observed in both Tibetan and Andean populations. Because it is also an important player in the circadian rhythm pathway, BHLHE41 might also provide insights into the mechanisms underlying the recognized impacts of hypoxia on the circadian clock. These results support the view that Ethiopian, Andean, and Tibetan populations living at high altitude have adapted to hypoxia differently, with convergent evolution affecting different genes from the same pathway.

Conditional analysis identifies three novel major histocompatibility complex loci associated with psoriasis.

Psoriasis is a common, chronic, inflammatory skin disorder. A number of genetic loci have been shown to confer risk for psoriasis. Collectively, these offer an integrated model for the inherited basis for susceptibility to psoriasis that combines altered skin barrier function together with the dysregulation of innate immune pathogen sensing and adap-tive immunity. The major histocompatibility complex (MHC) harbours the psoriasis susceptibility region which exhibits the largest effect size, driven in part by variation contained on the HLA-Cw*0602 allele. However, the resolution of the number and genomic location of potential independent risk loci are hampered by extensive linkage disequilibrium across the region. We leveraged the power of large psoriasis case and control data sets and the statistical approach of conditional analysis to identify potential further association signals distributed across the MHC. In addition to the major loci at HLA-C (P = 2.20 × 10(-236)), we observed and replicated four additional independent signals for disease association, three of which are novel. We detected evidence for association at SNPs rs2507971 (P = 6.73 × 10(-14)), rs9260313 (P = 7.93 × 10(-09)), rs66609536 (P = 3.54 × 10(-07)) and rs380924 (P = 6.24 × 10(-06)), located within the class I region of the MHC, with each observation replicated in an independent sample (P ≤ 0.01). The previously identified locus is close to MICA, the other three lie near MICB, HLA-A and HCG9 (a non-coding RNA gene). The identification of disease associations with both MICA and MICB is particularly intriguing, since each encodes an MHC class I-related protein with potent immunological function.

MAPT expression and splicing is differentially regulated by brain region: relation to genotype and implication for tauopathies.

The MAPT (microtubule-associated protein tau) locus is one of the most remarkable in neurogenetics due not only to its involvement in multiple neurodegenerative disorders, including progressive supranuclear palsy, corticobasal degeneration, Parksinson's disease and possibly Alzheimer's disease, but also due its genetic evolution and complex alternative splicing features which are, to some extent, linked and so all the more intriguing. Therefore, obtaining robust information regarding the expression, splicing and genetic regulation of this gene within the human brain is of immense importance. In this study, we used 2011 brain samples originating from 439 individuals to provide the most reliable and coherent information on the regional expression, splicing and regulation of MAPT available to date. We found significant regional variation in mRNA expression and splicing of MAPT within the human brain. Furthermore, at the gene level, the regional distribution of mRNA expression and total tau protein expression levels were largely in agreement, appearing to be highly correlated. Finally and most importantly, we show that while the reported H1/H2 association with gene level expression is likely to be due to a technical artefact, this polymorphism is associated with the expression of exon 3-containing isoforms in human brain. These findings would suggest that contrary to the prevailing view, genetic risk factors for neurodegenerative diseases at the MAPT locus are likely to operate by changing mRNA splicing in different brain regions, as opposed to the overall expression of the MAPT gene.

A systematic evaluation of the performance and properties of the UK Biobank Polygenic Risk Score (PRS) Release.

We assess the UK Biobank (UKB) Polygenic Risk Score (PRS) Release, a set of PRSs for 28 diseases and 25 quantitative traits that has been made available on the individuals in UKB, using a unified pipeline for PRS evaluation. We also release a benchmarking software tool to enable like-for-like performance evaluation for different PRSs for the same disease or trait. Extensive benchmarking shows the PRSs in the UKB Release to outperform a broad set of 76 published PRSs. For many of the diseases and traits we also validate the PRS algorithms in a separate cohort (100,000 Genomes Project). The availability of PRSs for 53 traits on the same set of individuals also allows a systematic assessment of their properties, and the increased power of these PRSs increases the evidence for their potential clinical benefit.

A polygenic risk score added to a QRISK®2 cardiovascular disease risk calculator demonstrated robust clinical acceptance and clinical utility in the primary care setting.

AIMS: The aim of the study was to assess the real-world feasibility, acceptability, and impact of an integrated risk tool for cardiovascular disease (CVD IRT, combining the standard QRISK®2 risk algorithm with a polygenic risk score), implemented within routine primary practice in the UK National Health Service. METHODS AND RESULTS: The Healthcare Evaluation of Absolute Risk Testing Study (NCT05294419) evaluated participants undergoing primary care health checks. Both QRISK2 and CVD IRT scores were returned to the healthcare providers (HCPs), who then communicated the results to participants. The primary outcome of the study was feasibility of CVD IRT implementation. Secondary outcomes included changes in CVD risk (QRISK2 vs. CVD IRT) and impact of the CVD IRT on clinical decision-making. A total of 832 eligible participants (median age 55 years, 62% females, 97.5% White ethnicity) were enrolled across 12 UK primary care practices. Cardiovascular disease IRT scores were obtained on 100% of the blood samples. Healthcare providers stated that the CVD IRT could be incorporated into routine primary care in a straightforward manner in 90.7% of reports. Participants stated they were 'likely' or 'very likely' to recommend the use of this test to their family or friends in 86.9% of reports. Participants stated that the test was personally useful (98.8%) and that the results were easy to understand (94.6%). When CVD IRT exceeded QRISK2, HCPs planned changes in management for 108/388 (27.8%) of participants and 47% (62/132) of participants with absolute risk score changes of >2%. CONCLUSION: Amongst HCPs and participants who agreed to the trial of genetic data for refinement of clinical risk prediction in primary care, we observed that CVD IRT implementation was feasible and well accepted. The CVD IRT results were associated with planned changes in prevention strategies.

When a standard cardiovascular risk tool, as currently used in National Health Service Health Checks, was expanded to include genetic risk information, it was well accepted by both participants and healthcare providers and generated impactful changes in planned clinical decision-making.Most participants found the test useful and easy to understand, and healthcare providers found it straightforward to use in most cases.When risk was increased by the addition of genetic information, this influenced planned management decisions.

Initial assessment of the pathogenic mechanisms of the recently identified Alzheimer risk Loci.

Recent genome wide association studies have identified CLU, CR1, ABCA7 BIN1, PICALM and MS4A6A/MS4A6E in addition to the long established APOE, as loci for Alzheimer's disease. We have systematically examined each of these loci to assess whether common coding variability contributes to the risk of disease. We have also assessed the regional expression of all the genes in the brain and whether there is evidence of an eQTL explaining the risk. In agreement with other studies we find that coding variability may explain the ABCA7 association, but common coding variability does not explain any of the other loci. We were not able to show that any of the loci had eQTLs within the power of this study. Furthermore the regional expression of each of the loci did not match the pattern of brain regional distribution in Alzheimer pathology. Although these results are mainly negative, they allow us to start defining more realistic alternative approaches to determine the role of all the genetic loci involved in Alzheimer's disease.

Natural selection on EPAS1 (HIF2alpha) associated with low hemoglobin concentration in Tibetan highlanders.

By impairing both function and survival, the severe reduction in oxygen availability associated with high-altitude environments is likely to act as an agent of natural selection. We used genomic and candidate gene approaches to search for evidence of such genetic selection. First, a genome-wide allelic differentiation scan (GWADS) comparing indigenous highlanders of the Tibetan Plateau (3,200-3,500 m) with closely related lowland Han revealed a genome-wide significant divergence across eight SNPs located near EPAS1. This gene encodes the transcription factor HIF2alpha, which stimulates production of red blood cells and thus increases the concentration of hemoglobin in blood. Second, in a separate cohort of Tibetans residing at 4,200 m, we identified 31 EPAS1 SNPs in high linkage disequilibrium that correlated significantly with hemoglobin concentration. The sex-adjusted hemoglobin concentration was, on average, 0.8 g/dL lower in the major allele homozygotes compared with the heterozygotes. These findings were replicated in a third cohort of Tibetans residing at 4,300 m. The alleles associating with lower hemoglobin concentrations were correlated with the signal from the GWADS study and were observed at greatly elevated frequencies in the Tibetan cohorts compared with the Han. High hemoglobin concentrations are a cardinal feature of chronic mountain sickness offering one plausible mechanism for selection. Alternatively, as EPAS1 is pleiotropic in its effects, selection may have operated on some other aspect of the phenotype. Whichever of these explanations is correct, the evidence for genetic selection at the EPAS1 locus from the GWADS study is supported by the replicated studies associating function with the allelic variants.