In vitro data suggest a role for PMS2 Kozak sequence mutations in Lynch Syndrome risk

Lynch Syndrome (LS) is the most common hereditary cancer syndrome. Heterozygous loss-of-function variants in PMS2 are linked to LS. While these variants are not directly cancer causing, reduced PMS2 function results in the accumulation of somatic variants and increased cancer risk over time due to DNA mismatch repair dysfunction. It is reasonable that other types of genetic variation that impact the expression of PMS2 may also contribute to cancer risk. The Kozak sequence is a highly conserved translation initiation motif among higher eukaryotes and is defined as the nine basepairs upstream of the translation start codon through the first four bases of the translated sequence (5'-[GTT]GCATCCATGG-3'; human PMS2 NM_000535.7). While Kozak sequence variants in PMS2 have been reported in ClinVar in suspected hereditary cancer patients, all variants upstream of the translation start site are currently classified as variants of undetermined significance (VUSs). We hypothesized that variants significantly disrupting the Kozak sequence of PMS2 would decrease PMS2 protein expression, contributing to increased cancer risk over time. Using a dual luciferase reporter plasmid and site-directed mutagenesis, we generated the wild-type human PMS2 and the ClinVar VUSs within the PMS2 Kozak sequence. Besides the c.1A>C variant, which is already known to be pathogenic, we implicate six additional variants as ACMG/AMP pathogenic supporting (PP) variants and classify ten as benign supporting (BP). In summary, we present a method developed for the classification of human PMS2 Kozak sequence variants that can contribute to the re-classification of VUSs identified in patients.

Large-scale genome-wide association study of 398,238 women unveils seven novel loci associated with high-grade serous epithelial ovarian cancer risk

Background

Nineteen genomic regions have been associated with high-grade serous ovarian cancer (HGSOC). We used data from the Ovarian Cancer Association Consortium (OCAC), Consortium of Investigators of Modifiers of BRCA1/BRCA2 (CIMBA), UK Biobank (UKBB), and FinnGen to identify novel HGSOC susceptibility loci and develop polygenic scores (PGS).

Methods

We analyzed >22 million variants for 398,238 women. Associations were assessed separately by consortium and meta-analysed. OCAC and CIMBA data were used to develop PGS which were trained on FinnGen data and validated in UKBB and BioBank Japan.

Results

Eight novel variants were associated with HGSOC risk. An interesting discovery biologically was finding that TP53 3'-UTR SNP rs78378222 was associated with HGSOC (per T allele relative risk (RR)=1.44, 95%CI:1.28-1.62, P=1.76×10-9). The optimal PGS included 64,518 variants and was associated with an odds ratio of 1.46 (95%CI:1.37-1.54) per standard deviation in the UKBB validation (AUROC curve=0.61, 95%CI:0.59-0.62).

Conclusions

This study represents the largest GWAS for HGSOC to date. The results highlight that improvements in imputation reference panels and increased sample sizes can identify HGSOC associated variants that previously went undetected, resulting in improved PGS. The use of updated PGS in cancer risk prediction algorithms will then improve personalized risk prediction for HGSOC.

CAGI, the Critical Assessment of Genome Interpretation, establishes progress and prospects for computational genetic variant interpretation methods

BACKGROUND

The Critical Assessment of Genome Interpretation (CAGI) aims to advance the state-of-the-art for computational prediction of genetic variant impact, particularly where relevant to disease. The five complete editions of the CAGI community experiment comprised 50 challenges, in which participants made blind predictions of phenotypes from genetic data, and these were evaluated by independent assessors.

RESULTS

Performance was particularly strong for clinical pathogenic variants, including some difficult-to-diagnose cases, and extends to interpretation of cancer-related variants. Missense variant interpretation methods were able to estimate biochemical effects with increasing accuracy. Assessment of methods for regulatory variants and complex trait disease risk was less definitive and indicates performance potentially suitable for auxiliary use in the clinic.

CONCLUSIONS

Results show that while current methods are imperfect, they have major utility for research and clinical applications. Emerging methods and increasingly large, robust datasets for training and assessment promise further progress ahead.

Functional EPAS1/HIF2A missense variant is associated with hematocrit in Andean highlanders

Hypoxia-inducible factor pathway genes are linked to adaptation in both human and nonhuman highland species. EPAS1, a notable target of hypoxia adaptation, is associated with relatively lower hemoglobin concentration in Tibetans. We provide evidence for an association between an adaptive EPAS1 variant (rs570553380) and the same phenotype of relatively low hematocrit in Andean highlanders. This Andean-specific missense variant is present at a modest frequency in Andeans and absent in other human populations and vertebrate species except the coelacanth. CRISPR-base-edited human cells with this variant exhibit shifts in hypoxia-regulated gene expression, while metabolomic analyses reveal both genotype and phenotype associations and validation in a lowland population. Although this genocopy of relatively lower hematocrit in Andean highlanders parallels well-replicated findings in Tibetans, it likely involves distinct pathway responses based on a protein-coding versus noncoding variants, respectively. These findings illuminate how unique variants at EPAS1 contribute to the same phenotype in Tibetans and a subset of Andean highlanders despite distinct evolutionary trajectories.

Characterizing prostate cancer risk through multi-ancestry genome-wide discovery of 187 novel risk variants

The transferability and clinical value of genetic risk scores (GRSs) across populations remain limited due to an imbalance in genetic studies across ancestrally diverse populations. Here we conducted a multi-ancestry genome-wide association study of 156,319 prostate cancer cases and 788,443 controls of European, African, Asian and Hispanic men, reflecting a 57% increase in the number of non-European cases over previous prostate cancer genome-wide association studies. We identified 187 novel risk variants for prostate cancer, increasing the total number of risk variants to 451. An externally replicated multi-ancestry GRS was associated with risk that ranged from 1.8 (per standard deviation) in African ancestry men to 2.2 in European ancestry men. The GRS was associated with a greater risk of aggressive versus non-aggressive disease in men of African ancestry (P = 0.03). Our study presents novel prostate cancer susceptibility loci and a GRS with effective risk stratification across ancestry groups.

Detection of distant relatedness in biobanks for identification of undiagnosed carriers of a Mendelian disease variant: application to Long QT Syndrome

Rare genetic diseases are typically studied in referral populations, resulting in underdiagnosis and biased assessment of penetrance and phenotype. To address this, we developed a generalizable method of genotype inference based on distant relatedness and deployed this to identify undiagnosed Type 5 Long QT Syndrome (LQT5) rare variant carriers in a non-referral population. We identified 9 LQT5 families referred to a single specialty clinic, each carrying p.Asp76Asn, the most common LQT5 variant. We uncovered recent common ancestry and a single shared haplotype among probands. Application to a non-referral population of 69,819 BioVU biobank subjects identified 22 additional subjects sharing this haplotype, subsequently confirmed to carry p.Asp76Asn. Referral and non-referral carriers had prolonged QTc compared to controls, and, among carriers, QTc polygenic score additively associated with QTc prolongation. Thus, our novel analysis of shared chromosomal segments identified undiagnosed cases of genetic disease and refined the understanding of LQT5 penetrance and phenotype.

Evidence of Novel Susceptibility Variants for Prostate Cancer and a Multiancestry Polygenic Risk Score Associated with Aggressive Disease in Men of African Ancestry

BACKGROUND

Genetic factors play an important role in prostate cancer (PCa) susceptibility.

OBJECTIVE

To discover common genetic variants contributing to the risk of PCa in men of African ancestry.

DESIGN, SETTING, AND PARTICIPANTS

We conducted a meta-analysis of ten genome-wide association studies consisting of 19378 cases and 61620 controls of African ancestry.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

Common genotyped and imputed variants were tested for their association with PCa risk. Novel susceptibility loci were identified and incorporated into a multiancestry polygenic risk score (PRS). The PRS was evaluated for associations with PCa risk and disease aggressiveness.

RESULTS AND LIMITATIONS

Nine novel susceptibility loci for PCa were identified, of which seven were only found or substantially more common in men of African ancestry, including an African-specific stop-gain variant in the prostate-specific gene anoctamin 7 (ANO7). A multiancestry PRS of 278 risk variants conferred strong associations with PCa risk in African ancestry studies (odds ratios [ORs] >3 and >5 for men in the top PRS decile and percentile, respectively). More importantly, compared with men in the 40-60% PRS category, men in the top PRS decile had a significantly higher risk of aggressive PCa (OR = 1.23, 95% confidence interval = 1.10-1.38, p = 4.4 × 10-4).

CONCLUSIONS

This study demonstrates the importance of large-scale genetic studies in men of African ancestry for a better understanding of PCa susceptibility in this high-risk population and suggests a potential clinical utility of PRS in differentiating between the risks of developing aggressive and nonaggressive disease in men of African ancestry.

PATIENT SUMMARY

In this large genetic study in men of African ancestry, we discovered nine novel prostate cancer (PCa) risk variants. We also showed that a multiancestry polygenic risk score was effective in stratifying PCa risk, and was able to differentiate risk of aggressive and nonaggressive disease.

Germline genetic variants and pediatric rhabdomyosarcoma outcomes: a report from the Children's Oncology Group

BACKGROUND

Relative to other pediatric cancers, survival for rhabdomyosarcoma (RMS) has not improved in recent decades, suggesting the need to enhance risk stratification. Therefore, we conducted a genome-wide association study for event-free survival (EFS) and overall survival (OS) to identify genetic variants associated with outcomes in individuals with RMS.

METHODS

The study included 920 individuals with newly diagnosed RMS who were enrolled in Children's Oncology Group protocols. To assess the association of each single nucleotide polymorphism (SNP) with EFS and OS, we estimated hazard ratios (HRs) and 95% confidence intervals (CIs) using multivariable Cox proportional hazards models, adjusted for clinical covariates. All statistical tests were two sided. We also performed stratified analyses by histological subtype (alveolar and embryonal RMS) and carried out sensitivity analyses of statistically significant SNPs by PAX3/7-FOXO1 fusion status and genetic ancestry group.

RESULTS

We identified that rs17321084 was associated with worse EFS (HR = 2.01, 95% CI = 1.59 to 2.53, P = 5.39 × 10-9) and rs10094840 was associated with worse OS (HR = 1.84, 95% CI = 1.48 to 2.27, P = 2.13 × 10-8). Using publicly available data, we found that rs17321084 lies in a binding region for transcription factors GATA2 and GATA3, and rs10094840 is associated with SPAG1 and RNF19A expression. We also identified that CTNNA3 rs2135732 (HR = 3.75, 95% CI = 2.34 to 5.99, P = 3.54 × 10-8) and MED31 rs74504320 (HR = 3.21, 95% CI = 2.12 to 4.86, P = 3.60 × 10-8) were associated with worse OS among individuals with alveolar RMS.

CONCLUSIONS

We demonstrated that common germline variants are associated with EFS and OS among individuals with RMS. Additional replication and investigation of these SNP effects may further support their consideration in risk stratification protocols.

Detection of distant familial relatedness in biobanks for identification of undiagnosed carriers of a Mendelian disease variant: application to Long QT syndrome

Importance

The diagnosis and study of rare genetic disease is often limited to referral populations, leading to underdiagnosis and a biased assessment of penetrance and phenotype.

Objective

To develop a generalizable method of genotype inference based on distant relatedness and to deploy this to identify undiagnosed Type 5 Long QT Syndrome (LQT5) rare variant carriers in a non-referral population.

Participants

We identified 9 LQT5 probands and 3 first-degree relatives referred to a single Genetic Arrhythmia clinic, each carrying D76N (p.Asp76Asn), the most common variant implicated in LQT5. The non-referral population consisted of 69,879 ancestry-matched subjects in BioVU, a large biobank that links electronic health records to dense array data. Participants were enrolled from 2007-2022. Data analysis was performed in 2022.

Exposures

We developed and applied a novel approach to genotype inference (Distant Relatedness for Identification and Variant Evaluation, or DRIVE) to identify shared, identical-by-descent (IBD) large chromosomal segments in array data.

Main Outcomes and Measures

We sought to establish genetic relatedness among the probands and to use genomic segments underlying D76N to identify other potential carriers in BioVU. We then further studied the role of D76N in LQT5 pathogenesis.

Results

Genetic reconstruction of pedigrees and distant relatedness detection among clinic probands using DRIVE revealed shared recent common ancestry and identified a single long shared haplotype. Interrogation of the non-referral population in BioVU identified a further 23 subjects sharing this haplotype, and sequencing confirmed D76N carrier status in 22, all previously undiagnosed with LQT5. The QTc was prolonged in D76N carriers compared to BioVU controls, with 40% penetrance of QTc ≥ 480 msec. Among D76N carriers, a QTc polygenic score was additively associated with QTc prolongation.

Conclusions and Relevance

Detection of IBD shared chromosomal segments around D76N enabled identification of distantly related and previously undiagnosed rare-variant carriers, demonstrated the contribution of polygenic risk to monogenic disease penetrance, and further established LQT5 as a primary arrhythmia disorder. Analysis of shared chromosomal regions spanning disease-causing mutations can identify undiagnosed cases of genetic diseases.

Abstract 3508: Association between clonal hematopoiesis and risk of prostate cancer in a large sample of African ancestry men

Clonal hematopoiesis of indeterminate potential (CHIP) has been associated with inflammation, which is a risk factor for cancer, including prostate cancer. We previously reported weak evidence of an association between CHIP and prostate cancer risk in men of European ancestry. However, little is known for African ancestry populations. We investigated the association of age-related CHIP with overall and aggressive prostate cancer risk in a large whole-exome sequencing study of 12,049 African ancestry men, including 7,176 prostate cancer cases (of which 3,283 had aggressive disease and 1,074 had metastatic disease) and 4,873 controls. Somatic variant calling was carried out using GATK Mutect2, and only variants with minor allele frequencies (MAF) <0.1% and a variant allelic fraction (VAF) >5% were included. Variants with a MAF ≥0.1% in gnomAD were excluded. CHIP variants were identified from a list of pre-specified mutations in 74 genes. Associations were tested using regression models adjusting for age, sub-study, and top 10 principal components, with statistical significance tested by the likelihood ratio test and applying a Bonferroni correction to account for multiple testing. In total, 998 variants in 57 CHIP genes were identified. Consistent with previous results, we observed a strong association between CHIP and age at blood draw. CHIP genes in aggregate were not statistically significantly associated with risks of total (OR=1.12, 95% CI=0.97-1.28), aggressive (OR=1.14, 95% CI=0.92-1.43) or metastatic (OR=1.17, 95% CI=0.91-1.49) prostate cancer. We observed that carriers of variants in DNMT3A, which is the gene that harbors the most CHIP driver mutations, had a nominally elevated risk of prostate cancer compared to non-carriers (OR=1.35, 95% CI=1.08-1.68, p=0.007). Additionally, carriers of variants in EZH2, which is implicated in cancer progression, showed a suggestive association with aggressive prostate cancer (OR=7.33, 95% CI=1.01-53.21, p=0.029). After adjusting for age at blood draw, CHIP genes in aggregate were not associated with age at prostate cancer diagnosis. However, we found that EZH2 variants carriers were diagnosed 12.9 years earlier on average than non-carriers (95% CI=6.1-19.7, adjusted p=0.01). A prostate cancer polygenic risk score (PRS) constructed using 269 risk variants was not associated with CHIP carrier status in aggregate (OR=0.99, 95% CI=0.92-1.06, p=0.70) or with any individual gene (all adjusted p>0.05). In summary, overall CHIP is not likely to be a risk factor of prostate cancer or aggressive disease in men of African ancestry. However, our results do confirm the association of CHIP in DNMT3A with prostate cancer risk as reported in previous studies in men of European ancestry. Future work will be needed to evaluate the biological causality of DNMT3A- and EZH2- related CHIP on prostate cancer.

Citation Format: Anqi Wang, Yili Xu, Xin Sheng, Raymond Hughley, Ben Adusei, Mohamed Jalloh, Serigne Magueye Gueye, Andrew A Adjei, James Mensah, Pedro W. Fernandez, Akin Olupelumi Adebiyi, Oseremen Inokhoife Aisuodionoe-Shadrach, Lindsay Petersen, Maureen Joffe, Jeannette T. Bensen, James L. Mohler, Jack A. Taylor, Eboneé N. Butler, Sue A. Ingles, Benjamin A. Rybicki, Janet L. Stanford, Wei Zheng, Sonja I. Berndt, Chad D. Huff, Joseph Lachance, Luc Multigner, Caroline Andrews, Timothy R. Rebbeck, Laurent Brureau, Stephen J. Chanock, Adam de Smith, Fei Chen, Burcu F. Darst, David V. Conti, Christopher A. Haiman. Association between clonal hematopoiesis and risk of prostate cancer in a large sample of African ancestry men [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 3508.

Rewired m6A epitranscriptomic networks link mutant p53 to neoplastic transformation

N6-methyladenosine (m6A), one of the most prevalent mRNA modifications in eukaryotes, plays a critical role in modulating both biological and pathological processes. However, it is unknown whether mutant p53 neomorphic oncogenic functions exploit dysregulation of m6A epitranscriptomic networks. Here, we investigate Li-Fraumeni syndrome (LFS)-associated neoplastic transformation driven by mutant p53 in iPSC-derived astrocytes, the cell-of-origin of gliomas. We find that mutant p53 but not wild-type (WT) p53 physically interacts with SVIL to recruit the H3K4me3 methyltransferase MLL1 to activate the expression of m6A reader YTHDF2, culminating in an oncogenic phenotype. Aberrant YTHDF2 upregulation markedly hampers expression of multiple m6A-marked tumor-suppressing transcripts, including CDKN2B and SPOCK2, and induces oncogenic reprogramming. Mutant p53 neoplastic behaviors are significantly impaired by genetic depletion of YTHDF2 or by pharmacological inhibition using MLL1 complex inhibitors. Our study reveals how mutant p53 hijacks epigenetic and epitranscriptomic machinery to initiate gliomagenesis and suggests potential treatment strategies for LFS gliomas.

Polygenic risk impacts PDGFRA mutation penetrance in non-syndromic cleft lip and palate

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is a common, severe craniofacial malformation that imposes significant medical, psychosocial and financial burdens. NSCL/P is a multifactorial disorder with genetic and environmental factors playing etiologic roles. Currently, only 25% of the genetic variation underlying NSCL/P has been identified by linkage, candidate gene and genome-wide association studies. In this study, whole-genome sequencing and genome-wide genotyping followed by polygenic risk score (PRS) and linkage analyses were used to identify the genetic etiology of NSCL/P in a large three-generation family. We identified a rare missense variant in PDGFRA (c.C2740T; p.R914W) as potentially etiologic in a gene-based association test using pVAAST (P = 1.78 × 10-4) and showed decreased penetrance. PRS analysis suggested that variant penetrance was likely modified by common NSCL/P risk variants, with lower scores found among unaffected carriers. Linkage analysis provided additional support for PRS-modified penetrance, with a 7.4-fold increase in likelihood after conditioning on PRS. Functional characterization experiments showed that the putatively causal variant was null for signaling activity in vitro; further, perturbation of pdgfra in zebrafish embryos resulted in unilateral orofacial clefting. Our findings show that a rare PDGFRA variant, modified by additional common NSCL/P risk variants, have a profound effect on NSCL/P risk. These data provide compelling evidence for multifactorial inheritance long postulated to underlie NSCL/P and may explain some unusual familial patterns.

Prostate cancer risk stratification improvement across multiple ancestries with new polygenic hazard score

BACKGROUND

Prostate cancer risk stratification using single-nucleotide polymorphisms (SNPs) demonstrates considerable promise in men of European, Asian, and African genetic ancestries, but there is still need for increased accuracy. We evaluated whether including additional SNPs in a prostate cancer polygenic hazard score (PHS) would improve associations with clinically significant prostate cancer in multi-ancestry datasets.

METHODS

In total, 299 SNPs previously associated with prostate cancer were evaluated for inclusion in a new PHS, using a LASSO-regularized Cox proportional hazards model in a training dataset of 72,181 men from the PRACTICAL Consortium. The PHS model was evaluated in four testing datasets: African ancestry, Asian ancestry, and two of European Ancestry-the Cohort of Swedish Men (COSM) and the ProtecT study. Hazard ratios (HRs) were estimated to compare men with high versus low PHS for association with clinically significant, with any, and with fatal prostate cancer. The impact of genetic risk stratification on the positive predictive value (PPV) of PSA testing for clinically significant prostate cancer was also measured.

RESULTS

The final model (PHS290) had 290 SNPs with non-zero coefficients. Comparing, for example, the highest and lowest quintiles of PHS290, the hazard ratios (HRs) for clinically significant prostate cancer were 13.73 [95% CI: 12.43-15.16] in ProtecT, 7.07 [6.58-7.60] in African ancestry, 10.31 [9.58-11.11] in Asian ancestry, and 11.18 [10.34-12.09] in COSM. Similar results were seen for association with any and fatal prostate cancer. Without PHS stratification, the PPV of PSA testing for clinically significant prostate cancer in ProtecT was 0.12 (0.11-0.14). For the top 20% and top 5% of PHS290, the PPV of PSA testing was 0.19 (0.15-0.22) and 0.26 (0.19-0.33), respectively.

CONCLUSIONS

We demonstrate better genetic risk stratification for clinically significant prostate cancer than prior versions of PHS in multi-ancestry datasets. This is promising for implementing precision-medicine approaches to prostate cancer screening decisions in diverse populations.

A whole-exome case-control association study to characterize the contribution of rare coding variation to pancreatic cancer risk

Pancreatic cancer is a deadly disease that accounts for approximately 5% of cancer deaths worldwide, with a dismal 5-year survival rate of 10%. Known genetic risk factors explain only a modest proportion of the heritable risk of pancreatic cancer. We conducted a whole-exome case-control sequencing study in 1,591 pancreatic cancer cases and 2,134 cancer-free controls of European ancestry. In our gene-based analysis, ATM ranked first, with a genome-wide significant p value of 1 × 10^-8. The odds ratio for protein-truncating variants in ATM was 24, which is substantially higher than prior estimates, although ours includes a broad 95% confidence interval (4.0-1000). SIK3 was the second highest ranking gene (p = 3.84 × 10^-6, false discovery rate or FDR = 0.032). We observed nominally significant association signals in several genes of a priori interest, including BRCA2 (p = 4.3 × 10^-4), STK11 (p = 0.003), PALB2 (p = 0.019), and TP53 (p = 0.037), and reported risk estimates for known pathogenic variants and variants of uncertain significance (VUS) in these genes. The rare variants in established susceptibility genes explain approximately 24% of log familial relative risk, which is comparable to the contribution from established common susceptibility variants (17%). In conclusion, this study provides new insights into the genetic susceptibility of pancreatic cancer, refining rare variant risk estimates in known pancreatic cancer susceptibility genes and identifying SIK3 as a novel candidate susceptibility gene. This study highlights the prominent importance of ATM truncating variants and the underappreciated role of VUS in pancreatic cancer etiology.

A Rare Germline HOXB13 Variant Contributes to Risk of Prostate Cancer in Men of African Ancestry

A rare African ancestry-specific germline deletion variant in HOXB13 (X285K, rs77179853) was recently reported in Martinican men with early-onset prostate cancer. Given the role of HOXB13 germline variation in prostate cancer, we investigated the association between HOXB13 X285K and prostate cancer risk in a large sample of 22 361 African ancestry men, including 11 688 prostate cancer cases. The risk allele was present only in men of West African ancestry, with an allele frequency in men that ranged from 0.40% in Ghana and 0.31% in Nigeria to 0% in Uganda and South Africa, with a range of frequencies in men with admixed African ancestry from North America and Europe (0-0.26%). HOXB13 X285K was associated with 2.4-fold increased odds of prostate cancer (95% confidence interval [CI] = 1.5-3.9, p = 2 × 10^-4), with greater risk observed for more aggressive and advanced disease (Gleason ≥8: odds ratio [OR] = 4.7, 95% CI = 2.3-9.5, p = 2 × 10^-5; stage T3/T4: OR = 4.5, 95% CI = 2.0-10.0, p = 2 × 10^-4; metastatic disease: OR = 5.1, 95% CI = 1.9-13.7, p = 0.001). We estimated that the allele arose in West Africa 1500-4600 yr ago. Further analysis is needed to understand how the HOXB13 X285K variant impacts the HOXB13 protein and function in the prostate. Understanding who carries this mutation may inform prostate cancer screening in men of West African ancestry. PATIENT SUMMARY: A rare African ancestry-specific germline deletion in HOXB13, found only in men of West African ancestry, was reported to be associated with an increased risk of overall and advanced prostate cancer. Understanding who carries this mutation may help inform screening for prostate cancer in men of West African ancestry.

Pedigree reconstruction and distant pairwise relatedness estimation from genome sequence data: A demonstration in a population of rhesus macaques (Macaca mulatta)

A primary challenge in the analysis of free‐ranging animal populations is the accurate estimation of relatedness among individuals. Many aspects of population analysis rely on knowledge of relatedness patterns, including socioecology, demography, heritability and gene mapping analyses, wildlife conservation and the management of breeding colonies. Methods for determining relatedness using genome‐wide data have improved our ability to determine kinship and reconstruct pedigrees in humans. However, methods for reconstructing complex pedigree structures and estimating distant relatedness (beyond third‐degree) have not been widely applied to other species. We sequenced the genomes of 150 male rhesus macaques from the Tulane National Primate Research Center colony to estimate pairwise relatedness, reconstruct closely related pedigrees, estimate more distant relationships and augment colony records. Methods for determining relatedness developed for human genetic data were applied and evaluated in the analysis of nonhuman primates, including identity‐by‐descent‐based methods for pedigree reconstruction and shared segment‐based inference of more distant relatedness. We compared the genotype‐based pedigrees and estimated relationships to available colony pedigree records and found high concordance (95.5% agreement) between expected and identified relationships for close relatives. In addition, we detected distant relationships not captured in colony records, including some as distant as twelfth‐degree. Furthermore, while deep sequence coverage is preferable, we show that this approach can also provide valuable information when only low‐coverage (5×) sequence data is available. Our findings demonstrate the value of these methods for determination of relatedness in various animal populations, with diverse applications to conservation biology, evolutionary and ecological research and biomedical studies.

Host genetic effects in pneumonia

Given the coronavirus disease 2019 (COVID-19) pandemic, investigations into host susceptibility to infectious diseases and downstream sequelae have never been more relevant. Pneumonia is a lung disease that can cause respiratory failure and hypoxia and is a common complication of infectious diseases, including COVID-19. Few genome-wide association studies (GWASs) of host susceptibility and severity of pneumonia have been conducted. We performed GWASs of pneumonia susceptibility and severity in the Vanderbilt University biobank (BioVU) with linked electronic health records (EHRs), including Illumina Expanded Multi-Ethnic Global Array (MEGAEX)-genotyped European ancestry (EA, n= 69,819) and African ancestry (AA, n = 15,603) individuals. Two regions of large effect were identified: the CFTR locus in EA (rs113827944; OR = 1.84, p value = 1.2 × 10-36) and HBB in AA (rs334 [p.Glu7Val]; OR = 1.63, p value = 3.5 × 10-13). Mutations in these genes cause cystic fibrosis (CF) and sickle cell disease (SCD), respectively. After removing individuals diagnosed with CF and SCD, we assessed heterozygosity effects at our lead variants. Further GWASs after removing individuals with CF uncovered an additional association in R3HCC1L (rs10786398; OR = 1.22, p value = 3.5 × 10-8), which was replicated in two independent datasets: UK Biobank (n = 459,741) and 7,985 non-overlapping BioVU subjects, who are genotyped on arrays other than MEGAEX. This variant was also validated in GWASs of COVID-19 hospitalization and lung function. Our results highlight the importance of the host genome in infectious disease susceptibility and severity and offer crucial insight into genetic effects that could potentially influence severity of COVID-19 sequelae.

Trans-ancestry genome-wide association meta-analysis of prostate cancer identifies new susceptibility loci and informs genetic risk prediction

Prostate cancer is a highly heritable disease with large disparities in incidence rates across ancestry populations. We conducted a multiancestry meta-analysis of prostate cancer genome-wide association studies (107,247 cases and 127,006 controls) and identified 86 new genetic risk variants independently associated with prostate cancer risk, bringing the total to 269 known risk variants. The top genetic risk score (GRS) decile was associated with odds ratios that ranged from 5.06 (95% confidence interval (CI), 4.84-5.29) for men of European ancestry to 3.74 (95% CI, 3.36-4.17) for men of African ancestry. Men of African ancestry were estimated to have a mean GRS that was 2.18-times higher (95% CI, 2.14-2.22), and men of East Asian ancestry 0.73-times lower (95% CI, 0.71-0.76), than men of European ancestry. These findings support the role of germline variation contributing to population differences in prostate cancer risk, with the GRS offering an approach for personalized risk prediction.

A Germline Variant at 8q24 Contributes to Familial Clustering of Prostate Cancer in Men of African Ancestry

Although men of African ancestry have a high risk of prostate cancer (PCa), no genes or mutations have been identified that contribute to familial clustering of PCa in this population. We investigated whether the African ancestry-specific PCa risk variant at 8q24, rs72725854, is enriched in men with a PCa family history in 9052 cases, 143 cases from high-risk families, and 8595 controls of African ancestry. We found the risk allele to be significantly associated with earlier age at diagnosis, more aggressive disease, and enriched in men with a PCa family history (32% of high-risk familial cases carried the variant vs 23% of cases without a family history and 12% of controls). For cases with two or more first-degree relatives with PCa who had at least one family member diagnosed at age <60 yr, the odds ratios for TA heterozygotes and TT homozygotes were 3.92 (95% confidence interval [CI] = 2.13-7.22) and 33.41 (95% CI = 10.86-102.84), respectively. Among men with a PCa family history, the absolute risk by age 60 yr reached 21% (95% CI = 17-25%) for TA heterozygotes and 38% (95% CI = 13-65%) for TT homozygotes. We estimate that in men of African ancestry, rs72725854 accounts for 32% of the total familial risk explained by all known PCa risk variants. PATIENT SUMMARY: We found that rs72725854, an African ancestry-specific risk variant, is more common in men with a family history of prostate cancer and in those diagnosed with prostate cancer at younger ages. Men of African ancestry may benefit from the knowledge of their carrier status for this genetic risk variant to guide decisions about prostate cancer screening.

Abstract 3517: A germline variant at 8q24 contributes to familial clustering of prostate cancer in men of African ancestry

Germline variation at 8q24 is the strongest risk factor for prostate cancer (PCa) across all racial and ethnic populations. While most 8q24 associations have been observed across populations, rs72725854 [T risk allele frequency ~6%] is only found in men of African ancestry and is the strongest known genome-wide association signal for PCa in this population. We investigated whether the T allele of rs72725854 is associated with PCa family history and age at diagnosis, characteristics known to have a strong genetic component. Analyses were performed using a sample of 9,052 cases and 8,595 controls from the African Ancestry Prostate Cancer (AAPC) GWAS Consortium and the ELLIPSE/PRACTICAL OncoArray Consortium. Participants were unselected for PCa family history. Among cases, 23.7% carried at least one copy of the T allele versus 11.6% of controls. The OR was 2.29 (95% CI=2.10–2.49) for TA heterozygotes and 5.04 (95% CI=3.36–7.55) for TT homozygotes. The percentage of cases carrying the T allele was significantly greater for men with a PCa family history (27.4% vs. 22.7% without a family history, p=0.002) and for men diagnosed &lt;60 (28.2% vs. 21.6% if over ≥60, p=0.002). The mean age at diagnosis for men with the TT genotype was 61.1 years (sd=8.7), compared to 62.7 (sd=9.1) for TA heterozygotes and 64.3 (sd=8.9) for AA homozygotes (p=5.7E-14). Carrier frequency was highest among men with both a positive family history and an early diagnosis (30.8%). The T allele was also over-represented in cases with more advanced PCa, with carrier frequencies ranging from 26% for lethal PCa (metastatic disease, PSA&gt;100 ng/ml or death from PCa), 25.4% for high-risk disease (stage T3/T4, Gleason 8-10, or PSA=20-100 ng/ml), 24.6% for intermediate-risk disease (Gleason=7, stage T1/T2, and PSA=10-20 ng/ml), and 21.4% for low-risk disease (Gleason&lt;7, stage T1/T2, and PSA&lt;10 ng/ml) (p=0.026). We also examined whether the risk allele is over-represented in 144 men from PCa families, with multiple first- and/or second-degree relatives with PCa or men diagnosed with PCa ≤55 years old. Among affected probands, 32.7% carried the risk allele, with 3.5% being homozygous carriers. The OR for TA heterozygotes and TT homozygotes was 3.41 (95% CI=2.33–4.98) and 11.06 (95% CI=3.92–31.18), respectively. Among men without a family history, the absolute risk for PCa by age 60 for non-risk allele carriers was 4.3%, compared to 9.0% and 15.6% for TA heterozygotes and TT homozygotes, respectively. Absolute risks by age 60 were higher among men with a family history of prostate cancer, reaching 9.0% for non-risk allele carriers, compared to 20.8% and 37.7% for TA heterozygotes and TT homozygotes, respectively. Given the high PCa risk conveyed by rs72725854 and the greater frequency of the allele in men with more aggressive and lethal disease, carriers of the risk allele would benefit from earlier and more regular PSA screening.

Citation Format: Burcu F. Darst, Jeannette T. Bensen, Sue A. Ingles, Benjamin A. Rybicki, Barbara Nemesure, Esther M. John, Jay H. Fowke, Victoria L. Stevens, Sonja I. Berndt, Chad D. Huff, Jong Y. Park, Wei Zheng, Elaine A. Ostrander, Shiv Srivastava, John Carpten, Thomas A. Sellers, Maureen Sanderson, Dana C. Crawford, Olivier Cussenot, Jennifer Cullen, Rick A. Kittles, Jianfeng Xu, Zsofia Kote-Jarai, Luc Multigner, Marie-Elise Parent, Florence Menegaux, Geraldine Cancel-Tassin, Adam S. Kibel, Eric A. Klein, Phyllis J. Goodman, Jennifer J. Hu, Graham Casey, Anselm J. Hennis, Ian M. Thompson, Robin Leach, James L. Mohler, Elizabeth T. Fontham, Gary J. Smith, Jack A. Taylor, Rosalind A. Eeles, Laurent Brureau, Stephen J. Chanock, Stephen Watya, Janet L. Stanford, Diptasri Mandal, William B. Isaacs, Kathleen A. Cooney, William J. Blot, David V. Conti, Christopher A. Haiman. A germline variant at 8q24 contributes to familial clustering of prostate cancer in men of African ancestry [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 3517.

Integrated case-control and somatic-germline interaction analyses of soft-tissue sarcoma

PURPOSE

The contribution of rare genetic variation in the development of soft-tissue sarcoma (STS) remains underexplored. To address this gap, we conducted a whole-exome case-control and somatic-germline interaction study to identify and characterise STS susceptible genes.

METHODS

The study involved 219 STS cases from The Cancer Genome Atlas and 3507 controls. All cases and controls were matched genetically onEuropean ancestry based on the 1000 Genomes project. Cross-platform technological stratification was performed with XPAT and gene-based association tests with VAAST 2.

RESULTS

NF1 exhibited the strongest genome-wide signal across the six subtypes, with p=1×10^-5. We also observed nominally significant association signals for three additional genes of interest, TP53 (p=0.0025), RB1 (p=0.0281), and MSH2 (p=0.0085). BAG1, which has not previously been implicated in STS, exhibited the strongest genome-wide signal after NF1, with p=6×10^-5. The association signals for NF1 and MSH2 were driven primarily by truncating variants, with ORs of 39 (95% CI: 7.1 to 220) for NF1 and 33 (95% CI: 2.4 to 460) for MSH2. In contrast, the association signals for RB1 and BAG1 were driven primarily by predicted damaging missense variants, with estimated ORs of 12 (95% CI: 2.4 to 59) for RB1 and 20 (95% CI: 1.4 to 300) for BAG1.

CONCLUSIONS

Our results confirm that pathogenic variants in NF1, RB1 and TP53 confer large increases in the risk of developing multiple STS subtypes, provide support for the role of MSH2 in STS susceptibility and identify BAG1 as a novel candidate STS risk gene.

Large-scale analysis of acquired chromosomal alterations in non-tumor samples from patients with cancer

Mosaicism, the presence of subpopulations of cells bearing somatic mutations, is associated with disease and aging and has been detected in diverse tissues, including apparently normal cells adjacent to tumors. To analyze mosaicism on a large scale, we surveyed haplotype-specific somatic copy number alterations (sCNAs) in 1,708 normal-appearing adjacent-to-tumor (NAT) tissue samples from 27 cancer sites and in 7,149 blood samples from The Cancer Genome Atlas. We find substantial variation across tissues in the rate, burden and types of sCNAs, including those spanning entire chromosome arms. We document matching sCNAs in the NAT tissue and the adjacent tumor, suggesting a shared clonal origin, as well as instances in which both NAT tissue and tumor tissue harbor a gain of the same oncogene arising in parallel from distinct parental haplotypes. These results shed light on pan-tissue mutations characteristic of field cancerization, the presence of oncogenic processes adjacent to cancer cells.

Assessing the performance of in silico methods for predicting the pathogenicity of variants in the gene CHEK2, among Hispanic females with breast cancer

The availability of disease-specific genomic data is critical for developing new computational methods that predict the pathogenicity of human variants and advance the field of precision medicine. However, the lack of gold standards to properly train and benchmark such methods is one of the greatest challenges in the field. In response to this challenge, the scientific community is invited to participate in the Critical Assessment for Genome Interpretation (CAGI), where unpublished disease variants are available for classification by in silico methods. As part of the CAGI-5 challenge, we evaluated the performance of 18 submissions and three additional methods in predicting the pathogenicity of single nucleotide variants (SNVs) in checkpoint kinase 2 (CHEK2) for cases of breast cancer in Hispanic females. As part of the assessment, the efficacy of the analysis method and the setup of the challenge were also considered. The results indicated that though the challenge could benefit from additional participant data, the combined generalized linear model analysis and odds of pathogenicity analysis provided a framework to evaluate the methods submitted for SNV pathogenicity identification and for comparison to other available methods. The outcome of this challenge and the approaches used can help guide further advancements in identifying SNV-disease relationships.

Association Between Birth Defects and Cancer Risk Among Children and Adolescents in a Population-Based Assessment of 10 Million Live Births

Importance

Birth defects affect approximately 1 in 33 children. Some birth defects are known to be strongly associated with childhood cancer (eg, trisomy 21 and acute leukemia). However, comprehensive evaluations of childhood cancer risk in those with birth defects have been limited in previous studies by insufficient sample sizes.

Objectives

To identify specific birth defect-childhood cancer (BD-CC) associations and characterize cancer risk in children by increasing number of nonchromosomal birth defects.

Design, Setting, and Participants

This multistate, population-based registry linkage study pooled statewide data on births, birth defects, and cancer from Texas, Arkansas, Michigan, and North Carolina on 10 181 074 children born from January 1, 1992, to December 31, 2013. Children were followed up to 18 years of age for a diagnosis of cancer. Data were retrieved between September 26, 2016, and September 21, 2017, and data analysis was performed from September 2, 2017, to March 21, 2019.

Exposures

Birth defects diagnoses (chromosomal anomalies and nonchromosomal birth defects) recorded by statewide, population-based birth defects registries.

Main Outcomes and Measures

Cancer diagnosis before age 18 years, as recorded in state cancer registries. Cox regression models were used to generate hazard ratios (HRs) and 95% CIs to evaluate BD-CC associations and the association between number of nonchromosomal defects and cancer risk.

Results

Compared with children without any birth defects, children with chromosomal anomalies were 11.6 (95% CI, 10.4-12.9) times more likely to be diagnosed with cancer, whereas children with nonchromosomal birth defects were 2.5 (95% CI, 2.4-2.6) times more likely to be diagnosed with cancer before 18 years of age. An increasing number of nonchromosomal birth defects was associated with a corresponding increase in the risk of cancer. Children with 4 or more major birth defects were 5.9 (95% CI, 5.3-6.4) times more likely to be diagnosed with cancer compared with those without a birth defect. In the analysis of 72 specific BD-CC patterns, 40 HRs were statistically significant (adjusted P < .05) after accounting for multiple comparisons. Cancers most frequently associated with nonchromosomal defects were hepatoblastoma and neuroblastoma.

Conclusions and Relevance

Several significant and novel associations were observed between specific birth defects and cancers. Among children with nonchromosomal birth defects, the number of major birth defects diagnosed was significantly and directly associated with cancer risk. These findings could inform clinical treatment for children with birth defects and may elucidate mechanisms that lead to these complex outcomes.

Functionally oriented analysis of cardiometabolic traits in a trans-ethnic sample

Interpretation of genetic association results is difficult because signals often lack biological context. To generate hypotheses of the functional genetic etiology of complex cardiometabolic traits, we estimated the genetically determined component of gene expression from common variants using PrediXcan (1) and determined genes with differential predicted expression by trait. PrediXcan imputes tissue-specific expression levels from genetic variation using variant-level effect on gene expression in transcriptome data. To explore the value of imputed genetically regulated gene expression (GReX) models across different ancestral populations, we evaluated imputed expression levels for predictive accuracy genome-wide in RNA sequence data in samples drawn from European-ancestry and African-ancestry populations and identified substantial predictive power using European-derived models in a non-European target population. We then tested the association of GReX on 15 cardiometabolic traits including blood lipid levels, body mass index, height, blood pressure, fasting glucose and insulin, RR interval, fibrinogen level, factor VII level and white blood cell and platelet counts in 15 755 individuals across three ancestry groups, resulting in 20 novel gene-phenotype associations reaching experiment-wide significance across ancestries. In addition, we identified 18 significant novel gene-phenotype associations in our ancestry-specific analyses. Top associations were assessed for additional support via query of S-PrediXcan (2) results derived from publicly available genome-wide association studies summary data. Collectively, these findings illustrate the utility of transcriptome-based imputation models for discovery of cardiometabolic effect genes in a diverse dataset.

ALPK1 missense pathogenic variant in five families leads to ROSAH syndrome, an ocular multisystem autosomal dominant disorder

Purpose

To identify the molecular cause in five unrelated families with a distinct autosomal dominant ocular systemic disorder we called ROSAH syndrome due to clinical features of retinal dystrophy, optic nerve edema, splenomegaly, anhidrosis, and migraine headache.

Methods

Independent discovery exome and genome sequencing in families 1, 2, and 3, and confirmation in families 4 and 5. Expression of wild-type messenger RNA and protein in human and mouse tissues and cell lines. Ciliary assays in fibroblasts from affected and unaffected family members.

Results

We found the heterozygous missense variant in the ɑ-kinase gene, ALPK1, (c.710C>T, [p.Thr237Met]), segregated with disease in all five families. All patients shared the ROSAH phenotype with additional low-grade ocular inflammation, pancytopenia, recurrent infections, and mild renal impairment in some. ALPK1 was notably expressed in retina, retinal pigment epithelium, and optic nerve, with immunofluorescence indicating localization to the basal body of the connecting cilium of the photoreceptors, and presence in the sweat glands. Immunocytofluorescence revealed expression at the centrioles and spindle poles during metaphase, and at the base of the primary cilium. Affected family member fibroblasts demonstrated defective ciliogenesis.

Conclusion

Heterozygosity for ALPK1, p.Thr237Met leads to ROSAH syndrome, an autosomal dominant ocular systemic disorder.

Introgression of regulatory alleles and a missense coding mutation drive plumage pattern diversity in the rock pigeon

Birds and other vertebrates display stunning variation in pigmentation patterning, yet the genes controlling this diversity remain largely unknown. Rock pigeons (Columba livia) are fundamentally one of four color pattern phenotypes, in decreasing order of melanism: T-check, checker, bar (ancestral), or barless. Using whole-genome scans, we identified NDP as a candidate gene for this variation. Allele-specific expression differences in NDP indicate cis-regulatory divergence between ancestral and melanistic alleles. Sequence comparisons suggest that derived alleles originated in the speckled pigeon (Columba guinea), providing a striking example of introgression. In contrast, barless rock pigeons have an increased incidence of vision defects and, like human families with hereditary blindness, carry start-codon mutations in NDP. In summary, we find that both coding and regulatory variation in the same gene drives wing pattern diversity, and post-domestication introgression supplied potentially advantageous melanistic alleles to feral populations of this ubiquitous urban bird.

XPAT: a toolkit to conduct cross-platform association studies with heterogeneous sequencing datasets

High-throughput sequencing data are increasingly being made available to the research community for secondary analyses, providing new opportunities for large-scale association studies. However, heterogeneity in target capture and sequencing technologies often introduce strong technological stratification biases that overwhelm subtle signals of association in studies of complex traits. Here, we introduce the Cross-Platform Association Toolkit, XPAT, which provides a suite of tools designed to support and conduct large-scale association studies with heterogeneous sequencing datasets. XPAT includes tools to support cross-platform aware variant calling, quality control filtering, gene-based association testing and rare variant effect size estimation. To evaluate the performance of XPAT, we conducted case-control association studies for three diseases, including 783 breast cancer cases, 272 ovarian cancer cases, 205 Crohn disease cases and 3507 shared controls (including 1722 females) using sequencing data from multiple sources. XPAT greatly reduced Type I error inflation in the case-control analyses, while replicating many previously identified disease-gene associations. We also show that association tests conducted with XPAT using cross-platform data have comparable performance to tests using matched platform data. XPAT enables new association studies that combine existing sequencing datasets to identify genetic loci associated with common diseases and other complex traits.

Integrated case-control and somatic-germline interaction analyses of melanoma susceptibility genes

While a number of genes have been implicated in melanoma susceptibility, the role of protein-coding variation in melanoma development and progression remains underexplored. To better characterize the role of germline coding variation in melanoma, we conducted a whole-exome case-control and somatic-germline interaction study involving 322 skin cutaneous melanoma cases from The Cancer Genome Atlas and 3607 controls of European ancestry. We controlled for cross-platform technological stratification using XPAT and conducted gene-based association tests using VAAST 2. Four established melanoma susceptibility genes achieved nominal statistical significance, MC1R (p = .0014), MITF (p = .0165) BRCA2 (p = .0206), and MTAP (p = .0393). We also observed a suggestive association for FANCA (p = .002), a gene previously implicated in melanoma survival. The association signal for BRCA2 was driven primarily by likely gene disrupting (LGD) variants, with an Odds Ratio (OR) of 5.62 (95% Confidence Interval (CI) 1.03-30.1). In contrast, the association signals for MC1R and MITF were driven primarily by predicted pathogenic missense variants, with estimated ORs of 1.4 to 3.0 for MC1R and 4.1 for MITF. MTAP exhibited an excess of both LGD and predicted damaging missense variants among cases, with ORs of 5.62 and 3.72, respectively, although neither category was significant. For individuals with known or predicted damaging variants, age of disease onset was significantly lower for two of the four genes, MC1R (p = .005) and MTAP (p = .035). In an analysis of germline carrier status and overlapping copy number alterations, we observed no evidence to support a two-hit model of carcinogenesis in any of the four genes. Although MC1R carriers were represented proportionally among the four molecular tumor subtypes, these individuals accounted for 69% of ultraviolet (UV) radiation mutational signatures among triple-wild type tumors (p = .040), highlighting the increased sensitivity to UV exposure among individuals with loss-of-function variants in MC1R.

Evolutionary selected Tibetan variants of HIF pathway and risk of lung cancer

Tibetans existed in high altitude for ~25 thousand years and have evolutionary selected unique haplotypes assumed to be beneficial to hypoxic adaptation. EGLN1/PHD2 and EPAS1/HIF-2α, both crucial components of hypoxia sensing, are the two best-established loci contributing to high altitude adaptation. The co-adapted Tibetan-specific haplotype encoding for PHD2:p.[D4E/C127S] promotes increased HIF degradation under hypoxic conditions. The Tibetan-specific 200 kb EPAS1 haplotype introgressed from an archaic human population related to Denisovans which underwent evolutionary decay; however, the functional variant(s) responsible for high-altitude adaptation at EPAS1/HIF-2α have not yet been identified. Since HIF modulates the behavior of cancer cells, we hypothesized that these Tibetan selected genomic variants may modify cancer risk predisposition. Here, we ascertained the frequencies of EGLN1^D4E/C127S and EGLN1^C127S variants and ten EPAS1/HIF-2α variants in lung cancer patients and controls in Nepal, whose population consists of people with Indo-Aryan origin and Tibetan-related Mongoloid origin. We observed a significant association between the selected Tibetan EGLN1/PHD2 haplotype and lung cancer (p=0.0012 for D4E, p=0.0002 for C127S), corresponding to a two-fold increase in lung cancer risk. We also observed a two-fold or greater increased risk for two of the ten EPAS1/HIF-2α variants, although the association was not significant after correcting for multiple comparisons (p=0.12). Although these data cannot address the role of these genetic variants on lung cancer initiation or progression, we conclude that some selected Tibetan variants are strongly associated with a modified risk of lung cancer.

Gain-of-function EGLN1 prolyl hydroxylase (PHD2 D4E:C127S) in combination with EPAS1 (HIF-2α) polymorphism lowers hemoglobin concentration in Tibetan highlanders

Tibetans have lived at high altitude for generations and are thought to be genetically adapted to hypoxic environments. Most are protected from hypoxia-induced polycythemia, and a haplotype of EPAS1, encoding hypoxia-inducible factor (HIF-2α), has been associated with lower hemoglobin levels. We earlier reported a Tibetan-specific EGLN1 haplotype encoding PHD2 which abrogates HIF augmentation in hypoxia. We genotyped 347 Tibetan individuals from varying altitudes for both the Tibetan-specific EGLN1 haplotype and 10 candidate SNPs in the EPAS1 haplotype and correlated their association with hemoglobin levels. The effect of the EGLN1 haplotype on hemoglobin exhibited age dependency at low altitude, while at higher altitudes, it showed a trend to lower hemoglobin levels in the presence of the Tibetan-selected EPAS1 rs142764723 C/C allele. The observed gene-environment and gene-gene interactions and the moderate effect of the EGLN1 and EPAS1 haplotypes on hemoglobin indicate that other modifiers exist. It remains to be determined whether a blunting of erythropoiesis or other physiological consequences of HIF downregulation are the primary drivers of these genetic adaptations among Tibetans.

KEY MESSAGE

Most Tibetans are protected from polycythemia while living in high altitude. An EGLN1 co-adapted haplotype, EGLN1 c.12C>G, c.380G>C is uniquely Tibetan. The Tibetan EPAS1 haplotype has introgressed from the Denisovan genome. While EGLN1 and EPAS1 genotypes lower Hb, this study indicates additional Hb modifiers.

Evolutionary history of Tibetans inferred from whole-genome sequencing

The indigenous people of the Tibetan Plateau have been the subject of much recent interest because of their unique genetic adaptations to high altitude. Recent studies have demonstrated that the Tibetan EPAS1 haplotype is involved in high altitude-adaptation and originated in an archaic Denisovan-related population. We sequenced the whole-genomes of 27 Tibetans and conducted analyses to infer a detailed history of demography and natural selection of this population. We detected evidence of population structure between the ancestral Han and Tibetan subpopulations as early as 44 to 58 thousand years ago, but with high rates of gene flow until approximately 9 thousand years ago. The CMS test ranked EPAS1 and EGLN1 as the top two positive selection candidates, and in addition identified PTGIS, VDR, and KCTD12 as new candidate genes. The advantageous Tibetan EPAS1 haplotype shared many variants with the Denisovan genome, with an ancient gene tree divergence between the Tibetan and Denisovan haplotypes of about 1 million years ago. With the exception of EPAS1, we observed no evidence of positive selection on Denisovan-like haplotypes.

Discovery and fine-mapping of adiposity loci using high density imputation of genome-wide association studies in individuals of African ancestry: African Ancestry Anthropometry Genetics Consortium

Genome-wide association studies (GWAS) have identified >300 loci associated with measures of adiposity including body mass index (BMI) and waist-to-hip ratio (adjusted for BMI, WHR_adjBMI), but few have been identified through screening of the African ancestry genomes. We performed large scale meta-analyses and replications in up to 52,895 individuals for BMI and up to 23,095 individuals for WHR_adjBMI from the African Ancestry Anthropometry Genetics Consortium (AAAGC) using 1000 Genomes phase 1 imputed GWAS to improve coverage of both common and low frequency variants in the low linkage disequilibrium African ancestry genomes. In the sex-combined analyses, we identified one novel locus (TCF7L2/HABP2) for WHR_adjBMI and eight previously established loci at P < 5×10⁻⁸: seven for BMI, and one for WHR_adjBMI in African ancestry individuals. An additional novel locus (SPRYD7/DLEU2) was identified for WHR_adjBMI when combined with European GWAS. In the sex-stratified analyses, we identified three novel loci for BMI (INTS10/LPL and MLC1 in men, IRX4/IRX2 in women) and four for WHR_adjBMI (SSX2IP, CASC8, PDE3B and ZDHHC1/HSD11B2 in women) in individuals of African ancestry or both African and European ancestry. For four of the novel variants, the minor allele frequency was low (<5%). In the trans-ethnic fine mapping of 47 BMI loci and 27 WHR_adjBMI loci that were locus-wide significant (P < 0.05 adjusted for effective number of variants per locus) from the African ancestry sex-combined and sex-stratified analyses, 26 BMI loci and 17 WHR_adjBMI loci contained ≤ 20 variants in the credible sets that jointly account for 99% posterior probability of driving the associations. The lead variants in 13 of these loci had a high probability of being causal. As compared to our previous HapMap imputed GWAS for BMI and WHR_adjBMI including up to 71,412 and 27,350 African ancestry individuals, respectively, our results suggest that 1000 Genomes imputation showed modest improvement in identifying GWAS loci including low frequency variants. Trans-ethnic meta-analyses further improved fine mapping of putative causal variants in loci shared between the African and European ancestry populations.

Genome-wide association studies (GWAS) have identified >300 genetic regions that influence body size and shape as measured by body mass index (BMI) and waist-to-hip ratio (WHR), respectively, but few have been identified in populations of African ancestry. We conducted large scale high coverage GWAS and replication of these traits in 52,895 and 23,095 individuals of African ancestry, respectively, followed by additional replication in European populations. We identified 10 genome-wide significant loci in all individuals, and an additional seven loci by analyzing men and women separately. We combined African and European ancestry GWAS and were able to narrow down 43 out of 74 African ancestry associated genetic regions to contain small number of putative causal variants. Our results highlight the improvement of applying high density genome coverage and combining multiple ancestries in the identification and refinement of location of genetic regions associated with adiposity traits.

VARPRISM: incorporating variant prioritization in tests of de novo mutation association

Background

Patients with certain genetic diseases, such as autism spectrum disorder, have increased rates of de novo mutations within some protein-coding genes.

Results

We introduce the VARiant PRIoritization SuM (VARPRISM), a software package which incorporates functional variant prioritization information to improve the power to detect de novo mutations influencing disease risk. VARPRISM evaluates the consequence of any given exonic mutation on the protein sequence to estimate the likelihood that the mutation is benign or damaging and conducts a likelihood ratio test on the gene level. We analyzed the Simons Simplex Collection of 2508 parent-offspring autism trios using VARPRISM, replicating 44 genes previously implicated in autism susceptibility and identifying 20 additional candidate genes, including MYO1E, KCND3, PDCD1, DLX3, and TSPAN4 (false discovery rate < 0.3).

Conclusion

By incorporating functional predictions, VARPRISM improved the statistical power to identify de novo mutations increasing disease risks. VARPRISM is available at http://www.hufflab.org/software/VARPRISM.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-016-0341-9) contains supplementary material, which is available to authorized users.

Adaptive genetic changes related to haemoglobin concentration in native high-altitude Tibetans

NEW FINDINGS

What is the topic of this review? Tibetans have genetic adaptations that are hypothesized to underlie the distinct set of traits they exhibit at altitude. What advances does it highlight? Several adaptive signatures in the same genomic regions have been identified among Tibetan populations resident throughout the Qinghai-Tibetan Plateau. Many highland Tibetans exhibit a haemoglobin concentration within the range expected at sea level, and this trait is associated with putatively adaptive regions harbouring the hypoxia-inducible factor pathway genes EGLN1, EPAS1 and PPARA. Precise functional variants at adaptive loci and relationships to physiological traits, beyond haemoglobin concentration, are currently being examined in this population. Some native Tibetan, Andean and Ethiopian populations have lived at altitudes ranging from 3000 to >4000 m above sea level for hundreds of generations and exhibit distinct combinations of traits at altitude. It was long hypothesized that genetic factors contribute to adaptive differences in these populations, and recent advances in genomics provide evidence that some of the strongest signatures of positive selection in humans are those identified in Tibetans. Many of the top adaptive genomic regions highlighted thus far harbour genes related to hypoxia sensing and response. Putatively adaptive copies of three hypoxia-inducible factor pathway genes, EPAS1, EGLN1 and PPARA, are associated with sea-level range, rather than elevated, haemoglobin concentration observed in many Tibetans at high altitude, and recent studies provide insight into some of the precise adaptive variants, timing of adaptive events and functional roles. While several studies in highland Tibetans have converged on a few hypoxia-inducible factor pathway genes, additional candidates have been reported in independent studies of Tibetans located throughout the Qinghai-Tibetan Plateau. Various aspects of adaptive significance have yet to be identified, integrated, and fully explored. Given the rapid technological advances and interdisciplinary efforts in genomics, physiology and molecular biology, careful examination of Tibetans and comparisons with other distinctively adapted highland populations will provide valuable insight into evolutionary processes and models for both basic and clinical research.

PADRE: Pedigree-Aware Distant-Relationship Estimation

Accurate estimation of shared ancestry is an important component of many genetic studies; current prediction tools accurately estimate pairwise genetic relationships up to the ninth degree. Pedigree-aware distant-relationship estimation (PADRE) combines relationship likelihoods generated by estimation of recent shared ancestry (ERSA) with likelihoods from family networks reconstructed by pedigree reconstruction and identification of a maximum unrelated set (PRIMUS), improving the power to detect distant relationships between pedigrees. Using PADRE, we estimated relationships from simulated pedigrees and three extended pedigrees, correctly predicting 20% more fourth- through ninth-degree simulated relationships than when using ERSA alone. By leveraging pedigree information, PADRE can even identify genealogical relationships between individuals who are genetically unrelated. For example, although 95% of 13(th)-degree relatives are genetically unrelated, in simulations, PADRE correctly predicted 50% of 13(th)-degree relationships to within one degree of relatedness. The improvement in prediction accuracy was consistent between simulated and actual pedigrees. We also applied PADRE to the HapMap3 CEU samples and report new cryptic relationships and validation of previously described relationships between families. PADRE greatly expands the range of relationships that can be estimated by using genetic data in pedigrees.

The Relationship between Native American Ancestry, Body Mass Index and Diabetes Risk among Mexican-Americans

Higher body mass index (BMI) is a well-established risk factor for type 2 diabetes, and rates of obesity and type 2 diabetes are substantially higher among Mexican-Americans relative to non-Hispanic European Americans. Mexican-Americans are genetically diverse, with a highly variable distribution of Native American, European, and African ancestries. Here, we evaluate the role of Native American ancestry on BMI and diabetes risk in a well-defined Mexican-American population. Participants were randomly selected among individuals residing in the Houston area who are enrolled in the Mexican-American Cohort study. Using a custom Illumina GoldenGate Panel, we genotyped DNA from 4,662 cohort participants for 87 Ancestry-Informative Markers. On average, the participants were of 50.2% Native American ancestry, 42.7% European ancestry and 7.1% African ancestry. Using multivariate linear regression, we found BMI and Native American ancestry were inversely correlated; individuals with <20% Native American ancestry were 2.5 times more likely to be severely obese compared to those with >80% Native American ancestry. Furthermore, we demonstrated an interaction between BMI and Native American ancestry in diabetes risk among women; Native American ancestry was a strong risk factor for diabetes only among overweight and obese women (OR = 1.190 for each 10% increase in Native American ancestry). This study offers new insight into the complex relationship between obesity, genetic ancestry, and their respective effects on diabetes risk. Findings from this study may improve the diabetes risk prediction among Mexican-American individuals thereby facilitating targeted prevention strategies.

Rare mutations in RINT1 predispose carriers to breast and Lynch syndrome-spectrum cancers

Approximately half of the familial aggregation of breast cancer remains unexplained. A multiple-case breast cancer family exome-sequencing study identified three likely pathogenic mutations in RINT1 (NM_021930.4) not present in public sequencing databases: RINT1 c.343C>T (p.Q115X), c.1132_1134del (p.M378del), and c.1207G>T (p.D403Y). On the basis of this finding, a population-based case-control mutation-screening study was conducted that identified 29 carriers of rare (minor allele frequency < 0.5%), likely pathogenic variants: 23 in 1,313 early-onset breast cancer cases and six in 1,123 frequency-matched controls [OR, 3.24; 95% confidence interval (CI), 1.29-8.17; P = 0.013]. RINT1 mutation screening of probands from 798 multiple-case breast cancer families identified four additional carriers of rare genetic variants. Analysis of the incidence of first primary cancers in families of women carrying RINT1 mutations estimated that carriers were at increased risk of Lynch syndrome-spectrum cancers [standardized incidence ratio (SIR), 3.35; 95% CI, 1.7-6.0; P = 0.005], particularly for relatives diagnosed with cancer under the age of 60 years (SIR, 10.9; 95% CI, 4.7-21; P = 0.0003).

SIGNIFICANCE

The work described in this study adds RINT1 to the growing list of genes in which rare sequence variants are associated with intermediate levels of breast cancer risk. Given that RINT1 is also associated with a spectrum of cancers with mismatch repair defects, these findings have clinical applications and raise interesting biological questions.

A unified test of linkage analysis and rare-variant association for analysis of pedigree sequence data

High-throughput sequencing of related individuals has become an important tool for studying human disease. However, owing to technical complexity and lack of available tools, most pedigree-based sequencing studies rely on an ad hoc combination of suboptimal analyses. Here we present pedigree-VAAST (pVAAST), a disease-gene identification tool designed for high-throughput sequence data in pedigrees. pVAAST uses a sequence-based model to perform variant and gene-based linkage analysis. Linkage information is then combined with functional prediction and rare variant case-control association information in a unified statistical framework. pVAAST outperformed linkage and rare-variant association tests in simulations and identified disease-causing genes from whole-genome sequence data in three human pedigrees with dominant, recessive and de novo inheritance patterns. The approach is robust to incomplete penetrance and locus heterogeneity and is applicable to a wide variety of genetic traits. pVAAST maintains high power across studies of monogenic, high-penetrance phenotypes in a single pedigree to highly polygenic, common phenotypes involving hundreds of pedigrees.

Shared and unique signals of high-altitude adaptation in geographically distinct Tibetan populations

Recent studies have used a variety of analytical methods to identify genes targeted by selection in high-altitude populations located throughout the Tibetan Plateau. Despite differences in analytic strategies and sample location, hypoxia-related genes, including EPAS1 and EGLN1, were identified in multiple studies. By applying the same analytic methods to genome-wide SNP information used in our previous study of a Tibetan population (n = 31) from the township of Maduo, located in the northeastern corner of the Qinghai-Tibetan Plateau (4200 m), we have identified common targets of natural selection in a second geographically and linguistically distinct Tibetan population (n = 46) in the Tuo Tuo River township (4500 m). Our analyses provide evidence for natural selection based on iHS and XP-EHH signals in both populations at the p<0.02 significance level for EPAS1, EGLN1, HMOX2, and CYP17A1 and for PKLR, HFE, and HBB and HBG2, which have also been reported in other studies. We highlight differences (i.e., stratification and admixture) in the two distinct Tibetan groups examined here and report selection candidate genes common to both groups. These findings should be considered in the prioritization of selection candidate genes in future genetic studies in Tibet.

Relationship estimation from whole-genome sequence data

The determination of the relationship between a pair of individuals is a fundamental application of genetics. Previously, we and others have demonstrated that identity-by-descent (IBD) information generated from high-density single-nucleotide polymorphism (SNP) data can greatly improve the power and accuracy of genetic relationship detection. Whole-genome sequencing (WGS) marks the final step in increasing genetic marker density by assaying all single-nucleotide variants (SNVs), and thus has the potential to further improve relationship detection by enabling more accurate detection of IBD segments and more precise resolution of IBD segment boundaries. However, WGS introduces new complexities that must be addressed in order to achieve these improvements in relationship detection. To evaluate these complexities, we estimated genetic relationships from WGS data for 1490 known pairwise relationships among 258 individuals in 30 families along with 46 population samples as controls. We identified several genomic regions with excess pairwise IBD in both the pedigree and control datasets using three established IBD methods: GERMLINE, fastIBD, and ISCA. These spurious IBD segments produced a 10-fold increase in the rate of detected false-positive relationships among controls compared to high-density microarray datasets. To address this issue, we developed a new method to identify and mask genomic regions with excess IBD. This method, implemented in ERSA 2.0, fully resolved the inflated cryptic relationship detection rates while improving relationship estimation accuracy. ERSA 2.0 detected all 1(st) through 6(th) degree relationships, and 55% of 9(th) through 11(th) degree relationships in the 30 families. We estimate that WGS data provides a 5% to 15% increase in relationship detection power relative to high-density microarray data for distant relationships. Our results identify regions of the genome that are highly problematic for IBD mapping and introduce new software to accurately detect 1(st) through 9(th) degree relationships from whole-genome sequence data.

Detecting statistical interaction between somatic mutational events and germline variation from next-generation sequence data

The two-hit model of carcinogenesis provides a valuable framework for understanding the role of DNA repair and tumor suppressor genes in cancer development and progression. Under this model, tumor development can initiate from a single somatic mutation in individuals that inherit an inactivating germline variant. Although the two-hit model can be an overgeneralization, the tendency for the pattern of somatic mutations to differ in cancer patients that inherit predisposition alleles is a signal that can be used to identify and validate germline susceptibility variants. Here, we present the Somatic-Germline Interaction (SGI) tool, which is designed to identify statistical interaction between germline variants and somatic mutational events from next-generation sequence data. SGI interfaces with rare-variant association tests and variant classifiers to identify candidate germline susceptibility variants from case-control sequencing data. SGI then analyzes tumor-normal pair next-generation sequence data to evaluate evidence for somatic-germline interaction in each gene or pathway using two tests: the Allelic Imbalance Rank Sum (AIRS) test and the Somatic Mutation Interaction Test (SMIT). AIRS tests for preferential allelic imbalance to evaluate whether somatic mutational events tend to amplify candidate germline variants. SMIT evaluates whether somatic point mutations and small indels occur more or less frequently than expected in the presence of candidate germline variants. Both AIRS and SMIT control for heterogeneity in the mutational process resulting from regional variation in mutation rates and inter-sample variation in background mutation rates. The SGI test combines AIRS and SMIT to provide a single, unified measure of statistical interaction between somatic mutational events and germline variation. We show that the tests implemented in SGI have high power with relatively modest sample sizes in a wide variety of scenarios. We demonstrate the utility of SGI to increase the power of rare variant association studies in cancer and to validate the potential role in cancer causation of germline susceptibility variants.

VAAST 2.0: improved variant classification and disease-gene identification using a conservation-controlled amino acid substitution matrix

The need for improved algorithmic support for variant prioritization and disease-gene identification in personal genomes data is widely acknowledged. We previously presented the Variant Annotation, Analysis, and Search Tool (VAAST), which employs an aggregative variant association test that combines both amino acid substitution (AAS) and allele frequencies. Here we describe and benchmark VAAST 2.0, which uses a novel conservation-controlled AAS matrix (CASM), to incorporate information about phylogenetic conservation. We show that the CASM approach improves VAAST’s variant prioritization accuracy compared to its previous implementation, and compared to SIFT, PolyPhen-2, and MutationTaster. We also show that VAAST 2.0 outperforms KBAC, WSS, SKAT, and variable threshold (VT) using published case-control datasets for Crohn disease (NOD2), hypertriglyceridemia (LPL), and breast cancer (CHEK2). VAAST 2.0 also improves search accuracy on simulated datasets across a wide range of allele frequencies, population-attributable disease risks, and allelic heterogeneity, factors that compromise the accuracies of other aggregative variant association tests. We also demonstrate that, although most aggregative variant association tests are designed for common genetic diseases, these tests can be easily adopted as rare Mendelian disease-gene finders with a simple ranking-by-statistical-significance protocol, and the performance compares very favorably to state-of-art filtering approaches. The latter, despite their popularity, have suboptimal performance especially with the increasing case sample size.

Genomic diversity and evolution of the head crest in the rock pigeon

The geographic origins of breeds and the genetic basis of variation within the widely distributed and phenotypically diverse domestic rock pigeon (Columba livia) remain largely unknown. We generated a rock pigeon reference genome and additional genome sequences representing domestic and feral populations. We found evidence for the origins of major breed groups in the Middle East and contributions from a racing breed to North American feral populations. We identified the gene EphB2 as a strong candidate for the derived head crest phenotype shared by numerous breeds, an important trait in mate selection in many avian species. We also found evidence that this trait evolved just once and spread throughout the species, and that the crest originates early in development by the localized molecular reversal of feather bud polarity.

Metabolic insight into mechanisms of high-altitude adaptation in Tibetans

Recent studies have identified genes involved in high-altitude adaptation in Tibetans. Genetic variants/haplotypes within regions containing three of these genes (EPAS1, EGLN1, and PPARA) are associated with relatively decreased hemoglobin levels observed in Tibetans at high altitude, providing corroborative evidence for genetic adaptation to this extreme environment. The mechanisms that afford adaptation to high-altitude hypoxia, however, remain unclear. Considering the strong metabolic demands imposed by hypoxia, we hypothesized that a shift in fuel preference to glucose oxidation and glycolysis at the expense of fatty acid oxidation would improve adaptation to decreased oxygen availability. Correlations between serum free fatty acid and lactate concentrations in Tibetan groups living at high altitude and putatively selected haplotypes provide insight into this hypothesis. An EPAS1 haplotype that exhibits a signal of positive selection is significantly associated with increased lactate concentration, the product of anaerobic glycolysis. Furthermore, the putatively advantageous PPARA haplotype is correlated with serum free fatty acid concentrations, suggesting a possible decrease in the activity of fatty acid oxidation. Although further studies are required to assess the molecular mechanisms underlying these patterns, these associations suggest that genetic adaptation to high altitude involves alteration in energy utilization pathways.

Crohn's disease and genetic hitchhiking at IBD5

Inflammatory bowel disease 5 (IBD5) is a 250 kb haplotype on chromosome 5 that is associated with an increased risk of Crohn's disease in Europeans. The OCTN1 gene is centrally located on IBD5 and encodes a transporter of the antioxidant ergothioneine (ET). The 503F variant of OCTN1 is strongly associated with IBD5 and is a gain-of-function mutation that increases absorption of ET. Although 503F has been implicated as the variant potentially responsible for Crohn's disease susceptibility at IBD5, there is little evidence beyond statistical association to support its role in disease causation. We hypothesize that 503F is a recent adaptation in Europeans that swept to relatively high frequency and that disease association at IBD5 results not from 503F itself, but from one or more nearby hitchhiking variants, in the genes IRF1 or IL5. To test for evidence of recent positive selection on the 503F allele, we employed the iHS statistic, which was significant in the European CEU HapMap population (P=0.0007) and European Human Genome Diversity Panel populations (P≤0.01). To evaluate the hypothesis of disease-variant hitchhiking, we performed haplotype association tests on high-density microarray data in a sample of 1,868 Crohn's disease cases and 5,550 controls. We found that 503F haplotypes with recombination breakpoints between OCTN1 and IRF1 or IL5 were not associated with disease (odds ratio [OR]: 1.05, P=0.21). In contrast, we observed strong disease association for 503F haplotypes with no recombination between these three genes (OR: 1.24, P=2.6×10(-8)), as expected if the sweeping haplotype harbored one or more disease-causing mutations in IRF1 or IL5. To further evaluate these disease-gene candidates, we obtained expression data from lower gastrointestinal biopsies of healthy individuals and Crohn's disease patients. We observed a 72% increase in gene expression of IRF1 among Crohn's disease patients (P=0.0006) and no significant difference in expression of OCTN1. Collectively, these data indicate that the 503F variant has increased in frequency due to recent positive selection and that disease-causing variants in linkage disequilibrium with 503F have hitchhiked to relatively high frequency, thus forming the IBD5 risk haplotype. Finally, our association results and expression data support IRF1 as a strong candidate for Crohn's disease causation.

Maximum-likelihood estimation of recent shared ancestry (ERSA)

Accurate estimation of recent shared ancestry is important for genetics, evolution, medicine, conservation biology, and forensics. Established methods estimate kinship accurately for first-degree through third-degree relatives. We demonstrate that chromosomal segments shared by two individuals due to identity by descent (IBD) provide much additional information about shared ancestry. We developed a maximum-likelihood method for the estimation of recent shared ancestry (ERSA) from the number and lengths of IBD segments derived from high-density SNP or whole-genome sequence data. We used ERSA to estimate relationships from SNP genotypes in 169 individuals from three large, well-defined human pedigrees. ERSA is accurate to within one degree of relationship for 97% of first-degree through fifth-degree relatives and 80% of sixth-degree and seventh-degree relatives. We demonstrate that ERSA's statistical power approaches the maximum theoretical limit imposed by the fact that distant relatives frequently share no DNA through a common ancestor. ERSA greatly expands the range of relationships that can be estimated from genetic data and is implemented in a freely available software package.

Ancestry of the Iban is predominantly Southeast Asian: genetic evidence from autosomal, mitochondrial, and Y chromosomes

Humans reached present-day Island Southeast Asia (ISEA) in one of the first major human migrations out of Africa. Population movements in the millennia following this initial settlement are thought to have greatly influenced the genetic makeup of current inhabitants, yet the extent attributed to different events is not clear. Recent studies suggest that south-to-north gene flow largely influenced present-day patterns of genetic variation in Southeast Asian populations and that late Pleistocene and early Holocene migrations from Southeast Asia are responsible for a substantial proportion of ISEA ancestry. Archaeological and linguistic evidence suggests that the ancestors of present-day inhabitants came mainly from north-to-south migrations from Taiwan and throughout ISEA approximately 4,000 years ago. We report a large-scale genetic analysis of human variation in the Iban population from the Malaysian state of Sarawak in northwestern Borneo, located in the center of ISEA. Genome-wide single-nucleotide polymorphism (SNP) markers analyzed here suggest that the Iban exhibit greatest genetic similarity to Indonesian and mainland Southeast Asian populations. The most common non-recombining Y (NRY) and mitochondrial (mt) DNA haplogroups present in the Iban are associated with populations of Southeast Asia. We conclude that migrations from Southeast Asia made a large contribution to Iban ancestry, although evidence of potential gene flow from Taiwan is also seen in uniparentally inherited marker data.