Application of Genome-Wide Single Nucleotide Polymorphism Typing: Simple Association and Beyond

The International HapMap Project and the arrival of technologies that type more than 100,000 SNPs in a single experiment have made genome-wide single nucleotide polymorphism (GW-SNP) assay a realistic endeavor. This has sparked considerable debate regarding the promise of GW-SNP typing to identify genetic association in disease. As has already been shown, this approach has the potential to localize common genetic variation underlying disease risk. The data provided from this technology also lends itself to several other lines of investigation; autozygosity mapping in consanguineous families and outbred populations, direct detection of structural variation, admixture analysis, and other population genetic approaches. In this review we will discuss the potential uses and practical application of GW-SNP typing including those above and beyond simple association testing

Chelator free gallium-68 radiolabelling of silica coated iron oxide nanorods via surface interactions

The commercial availability of combined magnetic resonance imaging (MRI)/positron emission tomography (PET) scanners for clinical use has increased demand for easily prepared agents which offer signal or contrast in both modalities. Herein we describe a new class of silica coated iron–oxide nanorods (NRs) coated with polyethylene glycol (PEG) and/or a tetraazamacrocyclic chelator (DO3A). Studies of the coated NRs validate their composition and confirm their properties as in vivo T₂ MRI contrast agents. Radiolabelling studies with the positron emitting radioisotope gallium-68 (t1/2 = 68 min) demonstrate that, in the presence of the silica coating, the macrocyclic chelator was not required for preparation of highly stable radiometal-NR constructs. In vivo PET-CT and MR imaging studies show the expected high liver uptake of gallium-68 radiolabelled nanorods with no significant release of gallium-68 metal ions, validating our innovation to provide a novel simple method for labelling of iron oxide NRs with a radiometal in the absence of a chelating unit that can be used for high sensitivity liver imaging

Exploring dementia and neuronal ceroid lipofuscinosis genes in 100 FTD-like patients from 6 towns and rural villages on the Adriatic Sea cost of Apulia

Frontotemporal dementia (FTD) refers to a complex spectrum of clinically and genetically heterogeneous disorders. Although fully penetrant mutations in several genes have been identified and can explain the pathogenic mechanisms underlying a great portion of the Mendelian forms of the disease, still a significant number of families and sporadic cases remains genetically unsolved. We performed whole exome sequencing in 100 patients with a late-onset and heterogeneous FTD-like clinical phenotype from Apulia and screened mendelian dementia and neuronal ceroid lipofuscinosis genes. We identified a nonsense mutation in SORL1 VPS domain (p.R744X), in 2 siblings displaying AD with severe language problems and primary progressive aphasia and a near splice-site mutation in CLCN6 (p.S116P) segregating with an heterogeneous phenotype, ranging from behavioural FTD to FTD with memory onset and to the logopenic variant of primary progressive aphasia in one family. Moreover 2 sporadic cases with behavioural FTD carried heterozygous mutations in the CSF1R Tyrosin kinase flanking regions (p.E573K and p.R549H). By contrast, only a minority of patients carried pathogenic C9orf72 repeat expansions (1%) and likely moderately pathogenic variants in GRN (p.C105Y, p.C389fs and p.C139R) (3%). In concert with recent studies, our findings support a common pathogenic mechanisms between FTD and neuronal ceroid lipofuscinosis and suggests that neuronal ceroid lipofuscinosis genes should be investigated also in dementia patients with predominant frontal symptoms and language impairments

Whole-genome sequencing to understand the genetic architecture of common gene expression and biomarker phenotypes.

Initial results from sequencing studies suggest that there are relatively few low-frequency (<5%) variants associated with large effects on common phenotypes. We performed low-pass whole-genome sequencing in 680 individuals from the InCHIANTI study to test two primary hypotheses: (i) that sequencing would detect single low-frequency-large effect variants that explained similar amounts of phenotypic variance as single common variants, and (ii) that some common variant associations could be explained by low-frequency variants. We tested two sets of disease-related common phenotypes for which we had statistical power to detect large numbers of common variant-common phenotype associations-11 132 cis-gene expression traits in 450 individuals and 93 circulating biomarkers in all 680 individuals. From a total of 11 657 229 high-quality variants of which 6 129 221 and 5 528 008 were common and low frequency (<5%), respectively, low frequency-large effect associations comprised 7% of detectable cis-gene expression traits [89 of 1314 cis-eQTLs at P < 1 × 10(-06) (false discovery rate ∼5%)] and one of eight biomarker associations at P < 8 × 10(-10). Very few (30 of 1232; 2%) common variant associations were fully explained by low-frequency variants. Our data show that whole-genome sequencing can identify low-frequency variants undetected by genotyping based approaches when sample sizes are sufficiently large to detect substantial numbers of common variant associations, and that common variant associations are rarely explained by single low-frequency variants of large effect

Genetic Variability in CLU and Its Association with Alzheimer's Disease

Background: Recently, two large genome wide association studies in Alzheimer disease (AD) have identified variants in three different genes (CLU, PICALM and CR1) as being associated with the risk of developing AD. The strongest association was reported for an intronic single nucleotide polymorphism (SNP) in CLU.Methodology/Principal Findings: To further characterize this association we have sequenced the coding region of this gene in a total of 495 AD cases and 330 healthy controls. A total of twenty-four variants were found in both cases and controls. For the changes found in more than one individual, the genotypic frequencies were compared between cases and controls. Coding variants were found in both groups (including a nonsense mutation in a healthy subject), indicating that the pathogenicity of variants found in this gene must be carefully evaluated. We found no common coding variant associated with disease. In order to determine if common variants at the CLU locus effect expression of nearby (cis) mRNA transcripts, an expression quantitative loci (eQTL) analysis was performed. No significant eQTL associations were observed for the SNPs previously associated with AD.Conclusions/Significance: We conclude that common coding variability at this locus does not explain the association, and that there is no large effect of common genetic variability on expression in brain tissue. We surmise that the most likely mechanism underpinning the association is either small effects of genetic variability on resting gene expression, or effects on damage induced expression of the protein

Measures of Autozygosity in Decline: Globalization, Urbanization, and Its Implications for Medical Genetics

This research investigates the influence of demographic factors on human genetic sub-structure. In our discovery cohort, we show significant demographic trends for decreasing autozygosity associated with population variation in chronological age. Autozygosity, the genomic signature of consanguinity, is identifiable on a genome-wide level as extended tracts of homozygosity. We identified an average of 28.6 tracts of extended homozygosity greater than 1 Mb in length in a representative population of 809 unrelated North Americans of European descent ranging in chronological age from 19–99 years old. These homozygous tracts made up a population average of 42 Mb of the genome corresponding to 1.6% of the entire genome, with each homozygous tract an average of 1.5 Mb in length. Runs of homozygosity are steadily decreasing in size and frequency as time progresses (linear regression, p<0.05). We also calculated inbreeding coefficients and showed a significant trend for population-wide increasing heterozygosity outside of linkage disequilibrium. We successfully replicated these associations in a demographically similar cohort comprised of a subgroup of 477 Baltimore Longitudinal Study of Aging participants. We also constructed statistical models showing predicted declining rates of autozygosity spanning the 20th century. These predictive models suggest a 14.0% decrease in the frequency of these runs of homozygosity and a 24.3% decrease in the percent of the genome in runs of homozygosity, as well as a 30.5% decrease in excess homozygosity based on the linkage pruned inbreeding coefficients. The trend for decreasing autozygosity due to panmixia and larger effective population sizes will likely affect the frequency of rare recessive genetic diseases in the future. Autozygosity has declined, and it seems it will continue doing so

Allelic heterogeneity and more detailed analyses of known loci explain additional phenotypic variation and reveal complex patterns of association

The identification of multiple signals at individual loci could explain additional phenotypic variance (‘missing heritability’) of common traits, and help identify causal genes. We examined gene expression levels as a model trait because of the large number of strong genetic effects acting in cis. Using expression profiles from 613 individuals, we performed genome-wide single nucleotide polymorphism (SNP) analyses to identify cis-expression quantitative trait loci (eQTLs), and conditional analysis to identify second signals. We examined patterns of association when accounting for multiple SNPs at a locus and when including additional SNPs from the 1000 Genomes Project. We identified 1298 cis-eQTLs at an approximate false discovery rate 0.01, of which 118 (9%) showed evidence of a second independent signal. For this subset of 118 traits, accounting for two signals resulted in an average 31% increase in phenotypic variance explained (Wilcoxon P< 0.0001). The association of SNPs with cis gene expression could increase, stay similar or decrease in significance when accounting for linkage disequilibrium with second signals at the same locus. Pairs of SNPs increasing in significance tended to have gene expression increasing alleles on opposite haplotypes, whereas pairs of SNPs decreasing in significance tended to have gene expression increasing alleles on the same haplotypes. Adding data from the 1000 Genomes Project showed that apparently independent signals could be potentially explained by a single association signal. Our results show that accounting for multiple variants at a locus will increase the variance explained in a substantial fraction of loci, but that allelic heterogeneity will be difficult to define without resequencing loci and functional work

Additional rare variant analysis in Parkinson's Disease cases with and without known pathogenic mutations: evidence for oligogenic inheritance

Oligogenic inheritance implies a role for several genetic factors in disease etiology. We studied oligogenic inheritance in Parkinson’s (PD) by assessing the potential burden of additional rare variants in established Mendelian genes and/or GBA, in individuals with and without a primary pathogenic genetic cause in two large independent cohorts totaling 7,900 PD cases and 6,166 controls. An excess (≥30%) of cases with a recognized primary genetic cause had ≥1 additional rare variants in Mendelian PD genes, as compared with no known mutation PD cases (17%) and unaffected controls (16%), supporting our hypothesis. Carriers of additional Mendelian gene variants have younger ages at onset (AAO). The effect of additional Mendelian variants in LRRK2 G2019S mutation carriers, of which ATP13A2 variation is particularly common, may account for some of the variation in penetrance. About 10% of no known mutation PD cases harbor a rare GBA variant compared to known pathogenic mutation PD cases (8%) and controls (5%), with carriers having earlier AAOs. Together, the data suggest that the oligogenic inheritance of rare Mendelian variants may be important in patient with a primary pathogenic cause, whereas GBA increases risk across all forms of PD. This study highlights potential genetic complexity of Mendelian PD. The identification of potential modifying variants provides new insights into disease mechanisms by potentially separating relevant from benign variants and by the interaction between genes in specific pathways. In the future this may be relevant to genetic testing and counselling of PD patients and their families

Imputation of variants from the 1000 Genomes Project modestly improves known associations and can identify low-frequency variant-phenotype associations undetected by HapMap based imputation

notes: PMCID: PMC3655956This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.Genome-wide association (GWA) studies have been limited by the reliance on common variants present on microarrays or imputable from the HapMap Project data. More recently, the completion of the 1000 Genomes Project has provided variant and haplotype information for several million variants derived from sequencing over 1,000 individuals. To help understand the extent to which more variants (including low frequency (1% ≤ MAF <5%) and rare variants (<1%)) can enhance previously identified associations and identify novel loci, we selected 93 quantitative circulating factors where data was available from the InCHIANTI population study. These phenotypes included cytokines, binding proteins, hormones, vitamins and ions. We selected these phenotypes because many have known strong genetic associations and are potentially important to help understand disease processes. We performed a genome-wide scan for these 93 phenotypes in InCHIANTI. We identified 21 signals and 33 signals that reached P<5×10(-8) based on HapMap and 1000 Genomes imputation, respectively, and 9 and 11 that reached a stricter, likely conservative, threshold of P<5×10(-11) respectively. Imputation of 1000 Genomes genotype data modestly improved the strength of known associations. Of 20 associations detected at P<5×10(-8) in both analyses (17 of which represent well replicated signals in the NHGRI catalogue), six were captured by the same index SNP, five were nominally more strongly associated in 1000 Genomes imputed data and one was nominally more strongly associated in HapMap imputed data. We also detected an association between a low frequency variant and phenotype that was previously missed by HapMap based imputation approaches. An association between rs112635299 and alpha-1 globulin near the SERPINA gene represented the known association between rs28929474 (MAF = 0.007) and alpha1-antitrypsin that predisposes to emphysema (P = 2.5×10(-12)). Our data provide important proof of principle that 1000 Genomes imputation will detect novel, low frequency-large effect associations

Integration of GWAS SNPs and tissue specific expression profiling reveal discrete eQTLs for human traits in blood and brain

Our knowledge of the transcriptome has become much more complex since the days of the central dogma of molecular biology. We now know that splicing takes place to create potentially thousands of isoforms from a single gene, and we know that RNA does not always faithfully recapitulate DNA if RNA editing occurs. Collectively, these observations show that the transcriptome is amazingly rich with intricate regulatory mechanisms for overall gene expression, splicing, and RNA editing. Genetic variability can play a role in controlling gene expression, which can be identified by examining expression quantitative trait loci (eQTLs). eQTLs are genomic regions where genetic variants, including single nucleotide polymorphisms (SNPs) show a statistical association with expression of mRNA transcripts. In humans, many SNPs are also associated with disease, and have been identified using genome wide association studies (GWAS) but the biological effects of those SNPs are usually not known. If SNPs found in GWAS are also found in eQTLs, then one could hypothesize that expression levels may contribute to disease risk. Performing eQTL analysis with GWAS SNPs in both blood and brain, specifically the frontal cortex and the cerebellum, we found both shared and tissue unique eQTLS. The identification of tissue-unique eQTLs supports the argument that choice of tissue type is important in eQTL studies (Paper I). Aging is a complex process with the mechanisms underlying aging still being poorly defined. There is evidence that the transcriptome changes with age, and hence we used the brain dataset from our first paper as a discovery set, with an additional replication dataset, to investigate any aging-gene expression associations. We found evidence that many genes were associated with aging. We further found that there were more statically significant expression changes in the frontal cortex versus the cerebellum, indicating that brain regions may age at different rates. As the brain is a heterogeneous tissue including both neurons and non-neuronal cells, we used LCM to capture Purkinje cells as a representative neuronal type and repeated the age analysis. Looking at the discovery, replication and Purkinje cell datasets we found five genes with strong, replicated evidence of age-expression associations (Paper II). Being able to capture and quantify the depth of the transcriptome has been a lengthy process starting with methods that could only measure a single gene to genome-wide techniques such as microarray. A recently developed technology, RNA-Seq, shows promise in its ability to capture expression, splicing, and editing and with its broad dynamic range quantification is accurate and reliable. RNA-Seq is, however, data intensive and a great deal of computational expertise is required to fully utilize the strengths of this method. We aimed to create a small, well-controlled, experiment in order to test the performance of this relatively new technology in the brain. We chose embryonic versus adult cerebral cortex, as mice are genetically homogenous and there are many known differences in gene expression related to brain development that we could use as benchmarks for analysis testing. We found a large number of differences in total gene expression between embryonic and adult brain. Rigorous technical and biological validation illustrated the accuracy and dynamic range of RNA-Seq. We were also able to interrogate differences in exon usage in the same dataset. Finally we were able to identify and quantify both well-known and novel A-to-I edit sites. Overall this project helped us develop the tools needed to build usable pipelines for RNA-Seq data processing (Paper III). Our studies in the developing brain (Paper III) illustrated that RNA-Seq was a useful unbiased method for investigating RNA editing. To extend this further, we utilized a genetically modified mouse model to study the transcriptomic role of the RNA editing enzyme ADAR2. We found that ADAR2 was important for editing of the coding region of mRNA as a large proportion of RNA editing sites in coding regions had a statistically significant decrease in editing percentages in Adar2 -/-Gria2 R/R mice versus controls. However, despite indications in the literature that ADAR2 may also be involved in splicing and expression regulatory machinery we found no changes in gene expression or exon utilization in Adar2 -/-Gria2 R/R mice as compared to their littermate controls (Paper IV). In our final study, based on the methods developed in Papers III and IV, we revisited the idea of age related gene expression associations from Paper II. We used a subset of human frontal cortices for RNA sequencing. Interestingly we found more gene expression changes with aging compared to the previous data using microarrays in Paper II. When the significant gene lists were analysed for gene ontology enrichment, we found that there was a large number of downregulated genes involved in synaptic function while those that were upregulated had enrichment in immune function. This dataset illustrates that the aging brain may be predisposed to the processes found in neurodegenerative diseases (Paper V)