Rare genetic variants underlie outlying levels of DNA methylation and gene-expression.

Testing the effect of rare variants on phenotypic variation is difficult due to the need for extremely large cohorts to identify associated variants given expected effect sizes. An alternative approach is to investigate the effect of rare genetic variants on DNA methylation (DNAm) as effect sizes are expected to be larger for molecular traits compared with complex traits. Here, we investigate DNAm in healthy ageing populations-the Lothian Birth Cohorts of 1921 and 1936-and identify both transient and stable outlying DNAm levels across the genome. We find an enrichment of rare genetic single nucleotide polymorphisms (SNPs) within 1 kb of DNAm sites in individuals with stable outlying DNAm, implying genetic control of this extreme variation. Using a family-based cohort, the Brisbane Systems Genetics Study, we observed increased sharing of DNAm outliers among more closely related individuals, consistent with these outliers being driven by rare genetic variation. We demonstrated that outlying DNAm levels have a functional consequence on gene expression levels, with extreme levels of DNAm being associated with gene expression levels toward the tails of the population distribution. This study demonstrates the role of rare SNPs in the phenotypic variation of DNAm and the effect of extreme levels of DNAm on gene expression.

A locus conferring tolerance to Theileria infection in African cattle.

East Coast fever, a tick-borne cattle disease caused by the Theileria parva parasite, is among the biggest natural killers of cattle in East Africa, leading to over 1 million deaths annually. Here we report on the genetic analysis of a cohort of Bos indicus (Boran) cattle demonstrating heritable tolerance to infection with T. parva (h2 = 0.65, s.e. 0.57). Through a linkage analysis we identify a 6 Mb genomic region on bovine chromosome 15 that is significantly associated with survival outcome following T. parva exposure. Testing this locus in an independent cohort of animals replicates this association with survival following T. parva infection. A stop gained variant in a paralogue of the FAF1 gene in this region was found to be highly associated with survival across both related and unrelated animals, with only one of the 20 homozygote carriers (T/T) of this change succumbing to the disease in contrast to 44 out of 97 animals homozygote for the reference allele (C/C). Consequently, we present a genetic locus linked to tolerance of one of Africa's most important cattle diseases, raising the promise of marker-assisted selection for cattle that are less susceptible to infection by T. parva.

Investigating the Role of Brain Natriuretic Peptide (BNP) and N-Terminal-proBNP in Thrombosis and Acute Ischemic Stroke Etiology.

The need for biomarkers for acute ischemic stroke (AIS) to understand the mechanisms implicated in pathological clot formation is critical. The levels of the brain natriuretic peptides known as brain natriuretic peptide (BNP) and NT-proBNP have been shown to be increased in patients suffering from heart failure and other heart conditions. We measured their expression in AIS clots of cardioembolic (CE) and large artery atherosclerosis (LAA) etiology, evaluating their location inside the clots, aiming to uncover their possible role in thrombosis. We analyzed 80 thrombi from 80 AIS patients in the RESTORE registry of AIS clots, 40 of which were of CE and 40 of LAA etiology. The localization of BNP and NT-BNP, quantified using immunohistochemistry and immunofluorescence, in AIS-associated white blood cell subtypes was also investigated. We found a statistically significant positive correlation between BNP and NT-proBNP expression levels (Spearman's rho = 0.668 p < 0.0001 *). We did not observe any statistically significant difference between LAA and CE clots in BNP expression (0.66 [0.13-3.54]% vs. 0.53 [0.14-3.07]%, p = 0.923) or in NT-proBNP expression (0.29 [0.11-0.58]% vs. 0.18 [0.05-0.51]%, p = 0.119), although there was a trend of higher NT-proBNP expression in the LAA clots. It was noticeable that BNP was distributed throughout the thrombus and especially within platelet-rich regions. However, NT-proBNP colocalized with neutrophils, macrophages, and T-lymphocytes, suggesting its association with the thrombo-inflammatory process.

Graphia: A platform for the graph-based visualisation and analysis of high dimensional data.

Graphia is an open-source platform created for the graph-based analysis of the huge amounts of quantitative and qualitative data currently being generated from the study of genomes, genes, proteins metabolites and cells. Core to Graphia's functionality is support for the calculation of correlation matrices from any tabular matrix of continuous or discrete values, whereupon the software is designed to rapidly visualise the often very large graphs that result in 2D or 3D space. Following graph construction, an extensive range of measurement algorithms, routines for graph transformation, and options for the visualisation of node and edge attributes are available, for graph exploration and analysis. Combined, these provide a powerful solution for the interpretation of high-dimensional data from many sources, or data already in the form of a network or equivalent adjacency matrix. Several use cases of Graphia are described, to showcase its wide range of applications in the analysis biological data. Graphia runs on all major desktop operating systems, is extensible through the deployment of plugins and is freely available to download from https://graphia.app/.

Mapping restricted introgression across the genomes of admixed indigenous African cattle breeds.

BACKGROUND: The genomes of indigenous African cattle are composed of components with Middle Eastern (taurine) and South Asian (indicine) origins, providing a valuable model to study hybridization and to identify genetic barriers to gene flow. In this study, we analysed indigenous African cattle breeds as models of hybrid zones, considering taurine and indicine samples as ancestors. In a genomic cline analysis of whole-genome sequence data, we considered over 8 million variants from 144 animals, which allows for fine-mapping of potential genomic incompatibilities at high resolution across the genome. RESULTS: We identified several thousand variants that had significantly steep clines ('SCV') across the whole genome, indicating restricted introgression. Some of the SCV were clustered into extended regions, with the longest on chromosome 7, spanning 725 kb and including 27 genes. We found that variants with a high phenotypic impact (e.g. indels, intra-genic and missense variants) likely represent greater genetic barriers to gene flow. Furthermore, our findings provide evidence that a large proportion of breed differentiation in African cattle could be linked to genomic incompatibilities and reproductive isolation. Functional evaluation of genes with SCV suggest that mitonuclear incompatibilities and genes associated with fitness (e.g. resistance to paratuberculosis) could account for restricted gene flow in indigenous African cattle. CONCLUSIONS: To our knowledge, this is the first time genomic cline analysis has been applied to identify restricted introgression in the genomes of indigenous African cattle and the results provide extended insights into mechanisms (e.g. genomic incompatibilities) contributing to hybrid differentiation. These results have important implications for our understanding of genetic incompatibilities and reproductive isolation and provide important insights into the impact of cross-breeding cattle with the aim of producing offspring that are both hardy and productive.

Using machine learning to detect the differential usage of novel gene isoforms.

BACKGROUND: Differential isoform usage is an important driver of inter-individual phenotypic diversity and is linked to various diseases and traits. However, accurately detecting the differential usage of different gene transcripts between groups can be difficult, in particular in less well annotated genomes where the spectrum of transcript isoforms is largely unknown. RESULTS: We investigated whether machine learning approaches can detect differential isoform usage based purely on the distribution of reads across a gene region. We illustrate that gradient boosting and elastic net approaches can successfully identify large numbers of genes showing potential differential isoform usage between Europeans and Africans, that are enriched among relevant biological pathways and significantly overlap those identified by previous approaches. We demonstrate that diversity at the 3' and 5' ends of genes are primary drivers of these differences between populations. CONCLUSION: Machine learning methods can effectively detect differential isoform usage from read fraction data, and can provide novel insights into the biological differences between groups.

Assessment of genotyping array performance for genome-wide association studies and imputation in African cattle.

BACKGROUND: In cattle, genome-wide association studies (GWAS) have largely focused on European or Asian breeds, using genotyping arrays that were primarily designed for European cattle. Because there is growing interest in performing GWAS in African breeds, we have assessed the performance of 23 commercial bovine genotyping arrays for capturing the diversity across African breeds and performing imputation. We used 409 whole-genome sequences (WGS) spanning global cattle breeds, and a real cohort of 2481 individuals (including African breeds) that were genotyped with the Illumina high-density (HD) array and the GeneSeek bovine 50 k array. RESULTS: We found that commercially available arrays were not effective in capturing variants that segregate among African indicine animals. Only 6% of these variants in high linkage disequilibrium (LD) (r2 > 0.8) were on the best performing arrays, which contrasts with the 17% and 25% in African and European taurine cattle, respectively. However, imputation from available HD arrays can successfully capture most variants (accuracies up to 0.93), mainly when using a global, not continent-specific, reference panel, which partially reflects the unusually high levels of admixture on the continent. When considering functional variants, the GGPF250 array performed best for tagging WGS variants and imputation. Finally, we show that imputation from low-density arrays can perform almost as well as HD arrays, if a two-stage imputation approach is adopted, i.e. first imputing to HD and then to WGS, which can potentially reduce the costs of GWAS. CONCLUSIONS: Our results show that the choice of an array should be based on a balance between the objective of the study and the breed/population considered, with the HD and BOS1 arrays being the best choice for both taurine and indicine breeds when performing GWAS, and the GGPF250 being preferable for fine-mapping studies. Moreover, our results suggest that there is no advantage to using the indicus-specific arrays for indicus breeds, regardless of the objective. Finally, we show that using a reference panel that better represents global bovine diversity improves imputation accuracy, particularly for non-European taurine populations.

Increased ultra-rare variant load in an isolated Scottish population impacts exonic and regulatory regions.

Human population isolates provide a snapshot of the impact of historical demographic processes on population genetics. Such data facilitate studies of the functional impact of rare sequence variants on biomedical phenotypes, as strong genetic drift can result in higher frequencies of variants that are otherwise rare. We present the first whole genome sequencing (WGS) study of the VIKING cohort, a representative collection of samples from the isolated Shetland population in northern Scotland, and explore how its genetic characteristics compare to a mainland Scottish population. Our analyses reveal the strong contributions played by the founder effect and genetic drift in shaping genomic variation in the VIKING cohort. About one tenth of all high-quality variants discovered are unique to the VIKING cohort or are seen at frequencies at least ten fold higher than in more cosmopolitan control populations. Multiple lines of evidence also suggest relaxation of purifying selection during the evolutionary history of the Shetland isolate. We demonstrate enrichment of ultra-rare VIKING variants in exonic regions and for the first time we also show that ultra-rare variants are enriched within regulatory regions, particularly promoters, suggesting that gene expression patterns may diverge relatively rapidly in human isolates.

Comparison of V̇o2peak and V̇o2max at Different Sampling Intervals in Collegiate Wrestlers.

ABSTRACT: Sell, KM, Ghigiarelli, JJ, Prendergast, JM, Ciani, GJ, Martin, J, and Gonzalez, AM. Comparison of V̇o2peak and V̇o2max at different sampling intervals in collegiate wrestlers. J Strength Cond Res 35(10): 2915-2917, 2021-The purpose of this study was to determine the difference in the highest oxygen uptake (V̇o2peak) achieved during a maximal effort graded exercise test (GXT) in 20 NCAA Division I male wrestlers using breath-by-breath (BbB) values to the maximal uptake averaged across different time- and breath-based oxygen consumption sampling intervals (V̇o2max). Given the need for aerobic fitness and anaerobic power in wrestling, the accurate determination of V̇o2max is imperative if it is to be used to identify current aerobic fitness and consequently guide sport-specific training programs to address weaknesses in this area. Each subject completed a cycle ergometer GXT during which BbB data were collected via indirect calorimetry and V̇o2peak determined as the highest value. V̇o2max was considered as the average value of 3-s, 5-s, 10-s, 20-s, and 30-s sampling, and 3-b, 7-b, and 11-b sampling during the GXT. Results show that the BbB V̇o2peak was significantly higher than the 5-s, 10-s, 20-s, 30-s, and 11-b (p < 0.05). The 3-b V̇o2max was significantly higher than the 20-s and 30-s V̇o2max values (p < 0.05). The underestimation of V̇o2peak for each time-based interval sampling approach compared with BbB V̇o2peak is consistent with previous research, but the comparison of BbB data to breath-based interval sampling has not been widely addressed in prior research. The use of a 7-b sampling interval for the determination of V̇o2max may be a promising approach to minimize the systematic errors associated with BbB or less frequent sampling intervals, but future research is needed to further support its application with elite athletic populations such as those in the current study.

Shared regulatory sites are abundant in the human genome and shed light on genome evolution and disease pleiotropy.

Large-scale gene expression datasets are providing an increasing understanding of the location of cis-eQTLs in the human genome and their role in disease. However, little is currently known regarding the extent of regulatory site-sharing between genes. This is despite it having potentially wide-ranging implications, from the determination of the way in which genetic variants may shape multiple phenotypes to the understanding of the evolution of human gene order. By first identifying the location of non-redundant cis-eQTLs, we show that regulatory site-sharing is a relatively common phenomenon in the human genome, with over 10% of non-redundant regulatory variants linked to the expression of multiple nearby genes. We show that these shared, local regulatory sites are linked to high levels of chromatin looping between the regulatory sites and their associated genes. In addition, these co-regulated gene modules are found to be strongly conserved across mammalian species, suggesting that shared regulatory sites have played an important role in shaping human gene order. The association of these shared cis-eQTLs with multiple genes means they also appear to be unusually important in understanding the genetics of human phenotypes and pleiotropy, with shared regulatory sites more often linked to multiple human phenotypes than other regulatory variants. This study shows that regulatory site-sharing is likely an underappreciated aspect of gene regulation and has important implications for the understanding of various biological phenomena, including how the two and three dimensional structures of the genome have been shaped and the potential causes of disease pleiotropy outside coding regions.

The host ubiquitin-dependent segregase VCP/p97 is required for the onset of human cytomegalovirus replication.

The human cytomegalovirus major immediate early proteins IE1 and IE2 are critical drivers of virus replication and are considered pivotal in determining the balance between productive and latent infection. IE1 and IE2 are derived from the same primary transcript by alternative splicing and regulation of their expression likely involves a complex interplay between cellular and viral factors. Here we show that knockdown of the host ubiquitin-dependent segregase VCP/p97, results in loss of IE2 expression, subsequent suppression of early and late gene expression and, ultimately, failure in virus replication. RNAseq analysis showed increased levels of IE1 splicing, with a corresponding decrease in IE2 splicing following VCP knockdown. Global analysis of viral transcription showed the expression of a subset of viral genes is not reduced despite the loss of IE2 expression, including UL112/113. Furthermore, Immunofluorescence studies demonstrated that VCP strongly colocalised with the viral replication compartments in the nucleus. Finally, we show that NMS-873, a small molecule inhibitor of VCP, is a potent HCMV antiviral with potential as a novel host targeting therapeutic for HCMV infection.

Comparison of Physical Fitness Parameters for Starters vs. Nonstarters in an NCAA Division I Men's Lacrosse Team.

Sell, KM, Prendergast, JM, Ghigiarelli, JJ, Gonzalez, AM, Biscardi, LM, Jajtner, AR, and Rothstein, AS. Comparison of physical fitness parameters for starters vs. nonstarters in an NCAA Division I men's lacrosse team. J Strength Cond Res 32(11): 3160-3168, 2018-The purpose of this study was to present a fitness profile of Division I male lacrosse players and compare the fitness attributes across different positions and starting status. Forty-one Division I men's lacrosse players (19.6 ± 1.6 years, 82.5 ± 9.5 kg, 182.0 ± 5.4 cm) volunteered to participate in the study. Fitness attributes assessed included aerobic fitness (1.5-mile run), muscular strength (1 repetition maximum bench press, squat, and hang clean), grip strength (hand dynamometer), explosive power (vertical jump), agility (3-cone drill, pro-agility), body composition (7-site skinfold), and speed (20- and 40-yard sprint). All testing was conducted by a certified strength and conditioning coach and occurred at the conclusion of pre-season training. The only significant difference across positions was for body mass, whereby defensemen were significantly heavier than attacking players (p < 0.05). Starters were significantly faster on the 3-cone drill, 20- and 40-yard sprint, and jumped significantly higher on the vertical jump compared with nonstarters (p < 0.05). Attributes pertaining to anaerobic fitness (speed, agility, and explosive power) may be better predictors of starting status than aerobic fitness in men's NCAA Division I lacrosse players. This differs from previous research on men's club lacrosse players where a difference in aerobic fitness and body composition was shown between starters and nonstarters. The normative data presented in this study may assist strength and conditioning coaches in the development of sport-specific training programs and motivate athletes toward achieving sport-specific fitness goals by helping identify areas of weakness before the start of the season.

A recurrent de novo mutation in ACTG1 causes isolated ocular coloboma.

Ocular coloboma (OC) is a defect in optic fissure closure and is a common cause of severe congenital visual impairment. Bilateral OC is primarily genetically determined and shows marked locus heterogeneity. Whole-exome sequencing (WES) was used to analyze 12 trios (child affected with OC and both unaffected parents). This identified de novo mutations in 10 different genes in eight probands. Three of these genes encoded proteins associated with actin cytoskeleton dynamics: ACTG1, TWF1, and LCP1. Proband-only WES identified a second unrelated individual with isolated OC carrying the same ACTG1 allele, encoding p.(Pro70Leu). Both individuals have normal neurodevelopment with no extra-ocular signs of Baraitser-Winter syndrome. We found this mutant protein to be incapable of incorporation into F-actin. The LCP1 and TWF1 variants each resulted in only minor disturbance of actin interactions, and no further plausibly causative variants were identified in these genes on resequencing 380 unrelated individuals with OC.

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.

hapbin: An Efficient Program for Performing Haplotype-Based Scans for Positive Selection in Large Genomic Datasets.

Understanding how the genome is shaped by selective processes forms an integral part of modern biology. However, as genomic datasets continue to grow larger it is becoming increasingly difficult to apply traditional statistics for detecting signatures of selection to these cohorts. There is therefore a pressing need for the development of the next generation of computational and analytical tools for detecting signatures of selection in large genomic datasets. Here, we present hapbin, an efficient multithreaded implementation of extended haplotype homzygosity-based statistics for detecting selection, which is up to 3,400 times faster than the current fastest implementations of these algorithms.

Genome-wide association scan identifies a colorectal cancer susceptibility locus on chromosome 8q24.

Using a multistage genetic association approach comprising 7,480 affected individuals and 7,779 controls, we identified markers in chromosomal region 8q24 associated with colorectal cancer. In stage 1, we genotyped 99,632 SNPs in 1,257 affected individuals and 1,336 controls from Ontario. In stages 2-4, we performed serial replication studies using 4,024 affected individuals and 4,042 controls from Seattle, Newfoundland and Scotland. We identified one locus on chromosome 8q24 and another on 9p24 having combined odds ratios (OR) for stages 1-4 of 1.18 (trend; P = 1.41 x 10(-8)) and 1.14 (trend; P = 1.32 x 10(-5)), respectively. Additional analyses in 2,199 affected individuals and 2,401 controls from France and Europe supported the association at the 8q24 locus (OR = 1.16, trend; 95% confidence interval (c.i.): 1.07-1.26; P = 5.05 x 10(-4)). A summary across all seven studies at the 8q24 locus was highly significant (OR = 1.17, c.i.: 1.12-1.23; P = 3.16 x 10(-11)). This locus has also been implicated in prostate cancer.

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.

Abundant pleiotropy in human complex diseases and traits.

We present a systematic review of pleiotropy among SNPs and genes reported to show genome-wide association with common complex diseases and traits. We find abundant evidence of pleiotropy; 233 (16.9%) genes and 77 (4.6%) SNPs show pleiotropic effects. SNP pleiotropic status was associated with gene location (p = 0.024; pleiotropic SNPs more often exonic [14.5% versus 4.9% for nonpleiotropic, trait-associated SNPs] and less often intergenic [15.8% versus 23.6%]), "predicted transcript consequence" (p = 0.001; pleiotropic SNPs more often predicted to be structurally deleterious [5% versus 0.4%] but not more often in regulatory sequences), and certain disease classes. We develop a method to calculate the likelihood that pleiotropic links between traits occurred more often than expected and demonstrate that this approach can identify etiological links that are already known (such as between fetal hemoglobin and malaria risk) and those that are not yet established (e.g., between plasma campesterol levels and gallstones risk; and between immunoglobulin A and juvenile idiopathic arthritis). Examples of pleiotropy will accumulate over time, but it is already clear that pleiotropy is a common property of genes and SNPs associated with disease traits, and this will have implications for identification of molecular targets for drug development, future genetic risk-profiling, and classification of diseases.

Widespread signatures of recent selection linked to nucleosome positioning in the human lineage.

In this study we investigated the strengths and modes of selection associated with nucleosome positioning in the human lineage through the comparison of interspecies and intraspecies rates of divergence. We identify significant evidence for both positive and negative selection linked to human nucleosome positioning for the first time, implicating a widespread and important role for DNA sequence in the location of well-positioned nucleosomes. Selection appears to be acting on particular base substitutions to maintain optimum GC compositions in core and linker regions, with, e.g., unexpectedly elevated rates of C→T substitutions during recent human evolution at linker regions 60-90 bp from the nucleosome dyad but significant depletion of the same substitutions within nucleosome core regions. These patterns are strikingly consistent with the known relationships between genomic sequence composition and nucleosome assembly. By stratifying nucleosomes according to the GC content of their genomic neighborhood, we also show that the strength and direction of selection detected is dictated by local GC content. Intriguingly these signatures of selection are not restricted to nucleosomes in close proximity to exons, suggesting the correct positioning of nucleosomes is not only important in and around coding regions. This analysis provides strong evidence that the genomic sequences associated with nucleosomes are not evolving neutrally, and suggests that underlying DNA sequence is an important factor in nucleosome positioning. Recent signatures of selection linked to genomic features as ubiquitous as the nucleosome have important implications for human genome evolution and disease.

Multiple common susceptibility variants near BMP pathway loci GREM1, BMP4, and BMP2 explain part of the missing heritability of colorectal cancer.

Genome-wide association studies (GWAS) have identified 14 tagging single nucleotide polymorphisms (tagSNPs) that are associated with the risk of colorectal cancer (CRC), and several of these tagSNPs are near bone morphogenetic protein (BMP) pathway loci. The penalty of multiple testing implicit in GWAS increases the attraction of complementary approaches for disease gene discovery, including candidate gene- or pathway-based analyses. The strongest candidate loci for additional predisposition SNPs are arguably those already known both to have functional relevance and to be involved in disease risk. To investigate this proposition, we searched for novel CRC susceptibility variants close to the BMP pathway genes GREM1 (15q13.3), BMP4 (14q22.2), and BMP2 (20p12.3) using sample sets totalling 24,910 CRC cases and 26,275 controls. We identified new, independent CRC predisposition SNPs close to BMP4 (rs1957636, P = 3.93×10(-10)) and BMP2 (rs4813802, P = 4.65×10(-11)). Near GREM1, we found using fine-mapping that the previously-identified association between tagSNP rs4779584 and CRC actually resulted from two independent signals represented by rs16969681 (P = 5.33×10(-8)) and rs11632715 (P = 2.30×10(-10)). As low-penetrance predisposition variants become harder to identify-owing to small effect sizes and/or low risk allele frequencies-approaches based on informed candidate gene selection may become increasingly attractive. Our data emphasise that genetic fine-mapping studies can deconvolute associations that have arisen owing to independent correlation of a tagSNP with more than one functional SNP, thus explaining some of the apparently missing heritability of common diseases.