Fine-mapping analysis including over 254,000 East Asian and European descendants identifies 136 putative colorectal cancer susceptibility genes.

Genome-wide association studies (GWAS) have identified more than 200 common genetic variants independently associated with colorectal cancer (CRC) risk, but the causal variants and target genes are mostly unknown. We sought to fine-map all known CRC risk loci using GWAS data from 100,204 cases and 154,587 controls of East Asian and European ancestry. Our stepwise conditional analyses revealed 238 independent association signals of CRC risk, each with a set of credible causal variants (CCVs), of which 28 signals had a single CCV. Our cis-eQTL/mQTL and colocalization analyses using colorectal tissue-specific transcriptome and methylome data separately from 1299 and 321 individuals, along with functional genomic investigation, uncovered 136 putative CRC susceptibility genes, including 56 genes not previously reported. Analyses of single-cell RNA-seq data from colorectal tissues revealed 17 putative CRC susceptibility genes with distinct expression patterns in specific cell types. Analyses of whole exome sequencing data provided additional support for several target genes identified in this study as CRC susceptibility genes. Enrichment analyses of the 136 genes uncover pathways not previously linked to CRC risk. Our study substantially expanded association signals for CRC and provided additional insight into the biological mechanisms underlying CRC development.

Deciphering colorectal cancer genetics through multi-omic analysis of 100,204 cases and 154,587 controls of European and east Asian ancestries.

Colorectal cancer (CRC) is a leading cause of mortality worldwide. We conducted a genome-wide association study meta-analysis of 100,204 CRC cases and 154,587 controls of European and east Asian ancestry, identifying 205 independent risk associations, of which 50 were unreported. We performed integrative genomic, transcriptomic and methylomic analyses across large bowel mucosa and other tissues. Transcriptome- and methylome-wide association studies revealed an additional 53 risk associations. We identified 155 high-confidence effector genes functionally linked to CRC risk, many of which had no previously established role in CRC. These have multiple different functions and specifically indicate that variation in normal colorectal homeostasis, proliferation, cell adhesion, migration, immunity and microbial interactions determines CRC risk. Crosstissue analyses indicated that over a third of effector genes most probably act outside the colonic mucosa. Our findings provide insights into colorectal oncogenesis and highlight potential targets across tissues for new CRC treatment and chemoprevention strategies.

Vitamin D treatment induces in vitro and ex vivo transcriptomic changes indicating anti-tumor effects.

Vitamin D deficiency is associated with risk of several common cancers, including colorectal cancer (CRC). Here we have utilized patient derived epithelial organoids (ex vivo) and CRC cell lines (in vitro) to show that calcitriol (1,25OHD) increased the expression of the CRC tumor suppressor gene, CDH1, at both the transcript and protein level. Whole genome expression analysis demonstrated significant differential expression of a further six genes after 1,25OHD treatment, including genes with established links to carcinogenesis GADD45, EFTUD1 and KIAA1199. Furthermore, gene ontologies relevant to carcinogenesis were enriched by 1,25OHD treatment (e.g., 'regulation of Wnt signaling pathway', 'regulation of cell death'), with common enriched processes across in vitro and ex vivo cultures including 'negative regulation of cell proliferation', 'regulation of cell migration' and 'regulation of cell differentiation'. Our results identify genes and pathways that are modifiable by calcitriol that have links to CRC tumorigenesis. Hence the findings provide potential mechanism to the epidemiological and clinical trial data indicating a causal association between vitamin D and CRC. We suggest there is strong rationale for further well-designed trials of vitamin D supplementation as a novel CRC chemopreventive and chemotherapeutic agent.

Oral vitamin D supplementation induces transcriptomic changes in rectal mucosa that are linked to anti-tumour effects.

BACKGROUND: The risk for several common cancers is influenced by the transcriptomic landscape of the respective tissue-of-origin. Vitamin D influences in vitro gene expression and cancer cell growth. We sought to determine whether oral vitamin D induces beneficial gene expression effects in human rectal epithelium and identify biomarkers of response. METHODS: Blood and rectal mucosa was sampled from 191 human subjects and mucosa gene expression (HT12) correlated with plasma vitamin D (25-OHD) to identify differentially expressed genes. Fifty subjects were then administered 3200IU/day oral vitamin D3 and matched blood/mucosa resampled after 12 weeks. Transcriptomic changes (HT12/RNAseq) after supplementation were tested against the prioritised genes for gene-set and GO-process enrichment. To identify blood biomarkers of mucosal response, we derived receiver-operator curves and C-statistic (AUC) and tested biomarker reproducibility in an independent Supplementation Trial (BEST-D). RESULTS: Six hundred twenty-nine genes were associated with 25-OHD level (P < 0.01), highlighting 453 GO-term processes (FDR<0.05). In the whole intervention cohort, vitamin D supplementation enriched the prioritised mucosal gene-set (upregulated gene-set P < 1.0E-07; downregulated gene-set P < 2.6E-05) and corresponding GO terms (P = 2.90E-02), highlighting gene expression patterns consistent with anti-tumour effects. However, only 9 individual participants (18%) showed a significant response (NM gene-set enrichment P < 0.001) to supplementation. Expression changes in HIPK2 and PPP1CC expression served as blood biomarkers of mucosal transcriptomic response (AUC=0.84 [95%CI 0.66-1.00]) and replicated in BEST-D trial subjects (HIPK2 AUC=0.83 [95%CI 0.77-0.89]; PPP1CC AUC=0.91 [95%CI 0.86-0.95]). CONCLUSIONS: Higher plasma 25-OHD correlates with rectal mucosa gene expression patterns consistent with anti-tumour effects, and this beneficial signature is induced by short-term vitamin D supplementation. Heterogenous gene expression responses to vitamin D may limit the ability of randomised trials to identify beneficial effects of supplementation on CRC risk. However, in the current study blood expression changes in HIPK2 and PPP1CC identify those participants with significant anti-tumour transcriptomic responses to supplementation in the rectum. These data provide compelling rationale for a trial of vitamin D and CRC prevention using easily assayed blood gene expression signatures as intermediate biomarkers of response.

Differential genetic influences over colorectal cancer risk and gene expression in large bowel mucosa.

Site-specific variation in colorectal cancer (CRC) incidence, biology and prognosis are poorly understood. We sought to determine whether common genetic variants influencing CRC risk might exhibit topographical differences on CRC risk through regional differences in effects on gene expression in the large bowel mucosa. We conducted a site-specific genetic association study (10 630 cases, 31 331 controls) to identify whether established risk variants exert differential effects on risk of proximal, compared to distal CRC. We collected normal colorectal mucosa and blood from 481 subjects and assessed mucosal gene expression using Illumina HumanHT-12v4 arrays in relation to germline genotype. Expression quantitative trait loci (eQTLs) were explored by anatomical location of sampling. The rs3087967 genotype (chr11q23.1 risk variant) exhibited significant site-specific effects-risk of distal CRC (odds ratio [OR] = 1.20, P = 8.20 × 10-20 ) with negligible effects on proximal CRC risk (OR = 1.05, P = .10). Expression of 1261 genes differed between proximal and distal colonic mucosa (top hit PRAC gene, fold-difference = 10, P = 3.48 × 10-57 ). In eQTL studies, rs3087967 genotype was associated with expression of 8 cis- and 21 trans-genes. Four of these (AKAP14, ADH5P4, ASGR2, RP11-342M1.7) showed differential effects by site, with strongest trans-eQTL signals in proximal colonic mucosa (eg, AKAP14, beta = 0.61, P = 5.02 × 10-5 ) and opposite signals in distal mucosa (AKAP14, beta = -0.17, P = .04). In summary, genetic variation at the chr11q23.1 risk locus imparts greater risk of distal rather than proximal CRC and exhibits site-specific differences in eQTL effects in normal mucosa. Topographical differences in genomic control over gene expression relevant to CRC risk may underlie site-specific variation in CRC. Results may inform individualised CRC screening programmes.

Colorectal cancer risk variants rs10161980 and rs7495132 are associated with cancer survival outcome by a recessive mode of inheritance.

Previous studies using additive genetic models failed to identify robust evidence of associations between colorectal cancer (CRC) risk variants and survival outcomes. However, additive models can be prone to false negative detection if the underlying inheritance mode is recessive. Here, we tested all currently known CRC-risk variants (n = 129) in a discovery analysis of 5675 patients from a Scottish cohort. Significant associations were then validated in 2474 CRC cases from UK Biobank. We found that the TT genotype of the intron variant rs7495132 in the CRTC3 gene was associated with clinically relevant poorer CRC-specific survival in both the discovery (hazard ratio [HR] = 1.97, 95% confidence interval [CI] = 1.41-2.74, P = 6.1 × 10-5 ) and validation analysis (HR = 1.69, 95% CI = 1.03-2.79, P = .038). In addition, the GG genotype of rs10161980 (intronic variant of AL139383.1 lncRNA) was associated with worse overall survival in the discovery cohort (HR = 1.24, 95% CI = 1.10-1.39, P = 3.4 × 10-4 ) and CRC-specific survival in the validation cohort (HR = 1.26, 95% CI = 1.01-1.56, P = .040). Our findings show that common genetic risk factors can also influence CRC survival outcome.

Gene Co-Expression Network Analysis Identifies Vitamin D-Associated Gene Modules in Adult Normal Rectal Epithelium Following Supplementation.

Colorectal cancer (CRC) is a common, multifactorial disease. While observational studies have identified an association between lower vitamin D and higher CRC risk, supplementation trials have been inconclusive and the mechanisms by which vitamin D may modulate CRC risk are not well understood. We sought to perform a weighted gene co-expression network analysis (WGCNA) to identify modules present after vitamin D supplementation (when plasma vitamin D level was sufficient) which were absent before supplementation, and then to identify influential genes in those modules. The transcriptome from normal rectal mucosa biopsies of 49 individuals free from CRC were assessed before and after 12 weeks of 3200IU/day vitamin D (Fultium-D3) supplementation using paired-end total RNAseq. While the effects on expression patterns following vitamin D supplementation were subtle, WGCNA identified highly correlated genes forming gene modules. Four of the 17 modules identified in the post-vitamin D network were not preserved in the pre-vitamin D network, shedding new light on the biochemical impact of supplementation. These modules were enriched for GO terms related to the immune system, hormone metabolism, cell growth and RNA metabolism. Across the four treatment-associated modules, 51 hub genes were identified, with enrichment of 40 different transcription factor motifs in promoter regions of those genes, including VDR:RXR. Six of the hub genes were nominally differentially expressed in studies of vitamin D effects on adult normal mucosa organoids: LCN2, HLA-C, AIF1L, PTPRU, PDE4B and IFI6. By taking a gene-correlation network approach, we have described vitamin D induced changes to gene modules in normal human rectal epithelium in vivo, the target tissue from which CRC develops.

Prediction of colorectal cancer risk based on profiling with common genetic variants.

Increasing numbers of common genetic variants associated with colorectal cancer (CRC) have been identified. Our study aimed to determine whether risk prediction based on common genetic variants might enable stratification for CRC risk. Meta-analysis of 11 genome-wide association studies comprising 16 871 cases and 26 328 controls was performed to capture CRC susceptibility variants. Genetic prediction models with several candidate polygenic risk scores (PRSs) were generated from Scottish CRC case-control studies (6478 cases and 11 043 controls) and the score with the best performance was then tested in UK Biobank (UKBB) (4800 cases and 20 287 controls). A weighted PRS of 116 CRC single nucleotide polymorphisms (wPRS116 ) was found with the best predictive performance, reporting a c-statistics of 0.60 and an odds ratio (OR) of 1.46 (95% confidence interval [CI] = 1.41-1.50, per SD increase) in Scottish data set. The predictive performance of this wPRS116 was consistently validated in UKBB data set with c-statistics of 0.61 and an OR of 1.49 (95% CI = 1.44-1.54, per SD increase). Modeling the levels of PRS with age and sex in the general UK population shows that employing genetic risk profiling can achieve a moderate degree of risk discrimination that could be helpful to identify a subpopulation with higher CRC risk due to genetic susceptibility.

Systematic meta-analyses, field synopsis and global assessment of the evidence of genetic association studies in colorectal cancer.

OBJECTIVE: To provide an understanding of the role of common genetic variations in colorectal cancer (CRC) risk, we report an updated field synopsis and comprehensive assessment of evidence to catalogue all genetic markers for CRC (CRCgene2). DESIGN: We included 869 publications after parallel literature review and extracted data for 1063 polymorphisms in 303 different genes. Meta-analyses were performed for 308 single nucleotide polymorphisms (SNPs) in 158 different genes with at least three independent studies available for analysis. Scottish, Canadian and Spanish data from genome-wide association studies (GWASs) were incorporated for the meta-analyses of 132 SNPs. To assess and classify the credibility of the associations, we applied the Venice criteria and Bayesian False-Discovery Probability (BFDP). Genetic associations classified as 'positive' and 'less-credible positive' were further validated in three large GWAS consortia conducted in populations of European origin. RESULTS: We initially identified 18 independent variants at 16 loci that were classified as 'positive' polymorphisms for their highly credible associations with CRC risk and 59 variants at 49 loci that were classified as 'less-credible positive' SNPs; 72.2% of the 'positive' SNPs were successfully replicated in three large GWASs and the ones that were not replicated were downgraded to 'less-credible' positive (reducing the 'positive' variants to 14 at 11 loci). For the remaining 231 variants, which were previously reported, our meta-analyses found no evidence to support their associations with CRC risk. CONCLUSION: The CRCgene2 database provides an updated list of genetic variants related to CRC risk by using harmonised methods to assess their credibility.

A Comprehensive Study of the Effect on Colorectal Cancer Survival of Common Germline Genetic Variation Previously Linked with Cancer Prognosis.

BACKGROUND: Germline genetic variants may influence pathways of tumor progression common to multiple cancer types. Here, we investigated the association between survival after colorectal cancer diagnosis and 128 common genetic variants previously associated with prognosis in genome-wide association studies in different cancer types. METHODS: We studied survival outcomes in a large well-documented, prospective, population-based cohort (5,675 patients with colorectal cancer) with up to 20 years' follow-up. RESULTS: None of the 128 variants were significantly associated with overall or colorectal cancer-specific survival (P < 5 × 10-4, Bonferroni-corrected threshold). We observed suggestive evidence (P < 0.05) for eight variants (rs17026425, rs17057166, rs6854845, rs1728400, rs17693104, rs202280, rs6797464, and rs823920) in all colorectal cancer and two variants (rs17026425 and rs6854845) in rectal cancer that were concordant with previous reports. CONCLUSIONS: Given good statistical power (>0.80 for 75% of variants), this study indicates that most previously reported variants associated with cancer survival have limited influence on colorectal cancer prognosis. IMPACT: Although small effects cannot be excluded, clinically meaningful germline influences on patients with colorectal cancer as a group are unlikely.

Effects of common genetic variants associated with colorectal cancer risk on survival outcomes after diagnosis: A large population-based cohort study.

Genome-wide association studies have thus far identified 130 genetic variants linked to colorectal cancer (CRC) risk (r2 < 0.2). Given their implication in disease causation, and thus plausible biologically effects on cancer-relevant biological pathways, we investigated whether these variants are associated with CRC prognosis and also whether they might provide predictive value for survival outcome. We conducted the analysis in a well-characterized population-based study of 5,675 patients after CRC diagnosis in Scotland. None of the genetic risk variants were associated with either overall survival (OS) or CRC-specific survival. Next, we combined the variants in a polygenic risk score, but again we observed no association between survival outcome and overall genetic susceptibility to CRC risk-as defined by common genetic variants (OS: hazard ratio = 1.00, 95% confidence interval = 0.96-1.05). Furthermore, we found no incremental increase in the discriminative performance when adding these genetic variants to the baseline CRC-survival predictive model of age, sex and stage at diagnosis. Given that our study is well-powered (>0.88) to detect effects on survival for 74% of the variants, we conclude that effects of common variants associated with CRC risk which have been identified to date are unlikely to have clinically relevant effect on survival outcomes for patients diagnosed with CRC.

Flexible and scalable diagnostic filtering of genomic variants using G2P with Ensembl VEP.

We aimed to develop an efficient, flexible and scalable approach to diagnostic genome-wide sequence analysis of genetically heterogeneous clinical presentations. Here we present G2P ( www.ebi.ac.uk/gene2phenotype ) as an online system to establish, curate and distribute datasets for diagnostic variant filtering via association of allelic requirement and mutational consequence at a defined locus with phenotypic terms, confidence level and evidence links. An extension to Ensembl Variant Effect Predictor (VEP), VEP-G2P was used to filter both disease-associated and control whole exome sequence (WES) with Developmental Disorders G2P (G2PDD; 2044 entries). VEP-G2PDD shows a sensitivity/precision of 97.3%/33% for de novo and 81.6%/22.7% for inherited pathogenic genotypes respectively. Many of the missing genotypes are likely false-positive pathogenic assignments. The expected number and discriminative features of background genotypes are defined using control WES. Using only human genetic data VEP-G2P performs well compared to other freely-available diagnostic systems and future phenotypic matching capabilities should further enhance performance.

Association analyses identify 31 new risk loci for colorectal cancer susceptibility.

Colorectal cancer (CRC) is a leading cause of cancer-related death worldwide, and has a strong heritable basis. We report a genome-wide association analysis of 34,627 CRC cases and 71,379 controls of European ancestry that identifies SNPs at 31 new CRC risk loci. We also identify eight independent risk SNPs at the new and previously reported European CRC loci, and a further nine CRC SNPs at loci previously only identified in Asian populations. We use in situ promoter capture Hi-C (CHi-C), gene expression, and in silico annotation methods to identify likely target genes of CRC SNPs. Whilst these new SNP associations implicate target genes that are enriched for known CRC pathways such as Wnt and BMP, they also highlight novel pathways with no prior links to colorectal tumourigenesis. These findings provide further insight into CRC susceptibility and enhance the prospects of applying genetic risk scores to personalised screening and prevention.

Genome-wide scan of the effect of common nsSNPs on colorectal cancer survival outcome.

BACKGROUND: We conducted a genome-wide scan to identify non-synonymous SNPs (nsSNPs) that might influence survival after a diagnosis of colorectal cancer (CRC). METHODS: We genotyped 7679 nsSNPs in 1939 Scottish patients from the Scottish Colorectal Cancer Study recruited soon after a CRC diagnosis and prospectively followed for survival outcomes. All-cause and CRC-specific survival analyses were conducted using Cox proportional hazard models adjusted for stage, age and sex for all cancer cases, after cancer type stratification and assuming additive and recessive models of inheritance. For all the SNPs that had a p-value < 0.10 a meta-analysis was performed combining the results of the discovery set and a replication set of 899 Scottish CRC patients. The p-value threshold of significance was set as at p < 10-8. RESULTS: 897 and 894 nsSNPs were associated with all-cause and CRC-specific mortality, respectively, at a p-value level < 0.10 in the discovery set. Meta-analysis of the results from the discovery and replication sets was performed overall and for cancers of colon and rectum separately and none of the variants reached a p-value < 10-8. CONCLUSIONS: This large scale well-powered analysis demonstrates that common nsSNPs are not associated with CRC prognosis overall.

Exploring causality in the association between circulating 25-hydroxyvitamin D and colorectal cancer risk: a large Mendelian randomisation study.

BACKGROUND: Whilst observational studies establish that lower plasma 25-hydroxyvitamin D (25-OHD) levels are associated with higher risk of colorectal cancer (CRC), establishing causality has proven challenging. Since vitamin D is modifiable, these observations have substantial clinical and public health implications. Indeed, many health agencies already recommend supplemental vitamin D. Here, we explore causality in a large Mendelian randomisation (MR) study using an improved genetic instrument for circulating 25-OHD. METHODS: We developed a weighted genetic score for circulating 25-OHD using six genetic variants that we recently reported to be associated with circulating 25-OHD in a large genome-wide association study (GWAS) meta-analysis. Using this score as instrumental variable in MR analyses, we sought to determine whether circulating 25-OHD is causally linked with CRC risk. We conducted MR analysis using individual-level data from 10,725 CRC cases and 30,794 controls (Scotland, UK Biobank and Croatia). We then applied estimates from meta-analysis of 11 GWAS of CRC risk (18,967 cases; 48,168 controls) in a summary statistics MR approach. RESULTS: The new genetic score for 25-OHD was strongly associated with measured plasma 25-OHD levels in 2821 healthy Scottish controls (P = 1.47 × 10- 11), improving upon previous genetic instruments (F-statistic 46.0 vs. 13.0). However, individual-level MR revealed no association between 25-OHD score and CRC risk (OR 1.03/unit log-transformed circulating 25-OHD, 95% CI 0.51-2.07, P = 0.93). Similarly, we found no evidence for a causal relationship between 25-OHD and CRC risk using summary statistics MR analysis (OR 0.91, 95% CI 0.69-1.19, P = 0.48). CONCLUSIONS: Despite the scale of this study and employing an improved score capturing more of the genetic contribution to circulating 25-OHD, we found no evidence for a causal relationship between circulating 25-OHD and CRC risk. Although the magnitude of effect for vitamin D suggested by observational studies can confidently be excluded, smaller effects sizes and non-linear relationships remain plausible. Circulating vitamin D may be a CRC biomarker, but a causal effect on CRC risk remains unproven.

Systematic meta-analyses and field synopsis of genetic association studies in colorectal adenomas.

BACKGROUND: Low penetrance genetic variants, primarily single nucleotide polymorphisms, have substantial influence on colorectal cancer (CRC) susceptibility. Most CRCs develop from colorectal adenomas (CRA). Here we report the first comprehensive field synopsis that catalogues all genetic association studies on CRA, with a parallel online database [http://www.chs.med.ed.ac.uk/CRAgene/]. METHODS: We performed a systematic review, reviewing 9750 titles, and then extracted data from 130 publications reporting on 181 polymorphisms in 74 genes. We conducted meta-analyses to derive summary effect estimates for 37 polymorphisms in 26 genes. We applied the Venice criteria and Bayesian False Discovery Probability (BFDP) to assess the levels of the credibility of associations. RESULTS: We considered the association with the rs6983267 variant at 8q24 as 'highly credible', reaching genome-wide statistical significance in at least one meta-analysis model. We identified 'less credible' associations (higher heterogeneity, lower statistical power, BFDP > 0.02) with a further four variants of four independent genes: MTHFR c.677C>T p.A222V (rs1801133), TP53 c.215C>G p.R72P (rs1042522), NQO1 c.559C>T p.P187S (rs1800566), and NAT1 alleles imputed as fast acetylator genotypes. For the remaining 32 variants of 22 genes for which positive associations with CRA risk have been previously reported, the meta-analyses revealed no credible evidence to support these as true associations. CONCLUSIONS: The limited number of credible associations between low penetrance genetic variants and CRA reflects the lower volume of evidence and associated lack of statistical power to detect associations of the magnitude typically observed for genetic variants and chronic diseases. The CRA gene database provides context for CRA genetic association data and will help inform future research directions.

Recurrent Coding Sequence Variation Explains Only A Small Fraction of the Genetic Architecture of Colorectal Cancer.

Whilst common genetic variation in many non-coding genomic regulatory regions are known to impart risk of colorectal cancer (CRC), much of the heritability of CRC remains unexplained. To examine the role of recurrent coding sequence variation in CRC aetiology, we genotyped 12,638 CRCs cases and 29,045 controls from six European populations. Single-variant analysis identified a coding variant (rs3184504) in SH2B3 (12q24) associated with CRC risk (OR = 1.08, P = 3.9 × 10(-7)), and novel damaging coding variants in 3 genes previously tagged by GWAS efforts; rs16888728 (8q24) in UTP23 (OR = 1.15, P = 1.4 × 10(-7)); rs6580742 and rs12303082 (12q13) in FAM186A (OR = 1.11, P = 1.2 × 10(-7) and OR = 1.09, P = 7.4 × 10(-8)); rs1129406 (12q13) in ATF1 (OR = 1.11, P = 8.3 × 10(-9)), all reaching exome-wide significance levels. Gene based tests identified associations between CRC and PCDHGA genes (P < 2.90 × 10(-6)). We found an excess of rare, damaging variants in base-excision (P = 2.4 × 10(-4)) and DNA mismatch repair genes (P = 6.1 × 10(-4)) consistent with a recessive mode of inheritance. This study comprehensively explores the contribution of coding sequence variation to CRC risk, identifying associations with coding variation in 4 genes and PCDHG gene cluster and several candidate recessive alleles. However, these findings suggest that recurrent, low-frequency coding variants account for a minority of the unexplained heritability of CRC.

Homozygous loss-of-function variants in European cosmopolitan and isolate populations.

Homozygous loss of function (HLOF) variants provide a valuable window on gene function in humans, as well as an inventory of the human genes that are not essential for survival and reproduction. All humans carry at least a few HLOF variants, but the exact number of inactivated genes that can be tolerated is currently unknown­as are the phenotypic effects of losing function for most human genes. Here, we make use of 1432 whole exome sequences from five European populations to expand the catalogue of known human HLOF mutations; after stringent filtering of variants in our dataset, we identify a total of 173 HLOF mutations, 76 (44%) of which have not been observed previously. We find that population isolates are particularly well suited to surveys of novel HLOF genes because individuals in such populations carry extensive runs of homozygosity, which we show are enriched for novel, rare HLOF variants. Further, we make use of extensive phenotypic data to show that most HLOFs, ascertained in population-based samples, appear to have little detectable effect on the phenotype. On the contrary, we document several genes directly implicated in disease that seem to tolerate HLOF variants. Overall HLOF genes are enriched for olfactory receptor function and are expressed in testes more often than expected, consistent with reduced purifying selection and incipient pseudogenisation.

Exome sequencing to detect rare variants associated with general cognitive ability: a pilot study.

Variation in human cognitive ability is of consequence to a large number of health and social outcomes and is substantially heritable. Genetic linkage, genome-wide association, and copy number variant studies have investigated the contribution of genetic variation to individual differences in normal cognitive ability, but little research has considered the role of rare genetic variants. Exome sequencing studies have already met with success in discovering novel trait-gene associations for other complex traits. Here, we use exome sequencing to investigate the effects of rare variants on general cognitive ability. Unrelated Scottish individuals were selected for high scores on a general component of intelligence (g). The frequency of rare genetic variants (in n = 146) was compared with those from Scottish controls (total n = 486) who scored in the lower to middle range of the g distribution or on a proxy measure of g. Biological pathway analysis highlighted enrichment of the mitochondrial inner membrane component and apical part of cell gene ontology terms. Global burden analysis showed a greater total number of rare variants carried by high g cases versus controls, which is inconsistent with a mutation load hypothesis whereby mutations negatively affect g. The general finding of greater non-synonymous (vs. synonymous) variant effects is in line with evolutionary hypotheses for g. Given that this first sequencing study of high g was small, promising results were found, suggesting that the study of rare variants in larger samples would be worthwhile.